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SYNOPSIS 

Mechanical behaviour of a granular material is 
controlled by i to constituent particle characteristics . The 
particle size, gradation, shape and angularity together with 
the relative density play an important role in the 
deformation characteristics of sand®. Although 

classification of a sand for engineering purposes is based on 
the size attributes, the presence of mineral constituents 
with compositional differences such as quartz, feldspar, mica 
and carbonates results in significant changes in the 
engineering behaviour of the system. While the Importance of 
all these parameters in understanding the behaviour of 
granular material has been recognized earlier, adequate 
attention has not been bestowed to study in detail the 
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influence of all these parameters on the system behaviour. 
Further , not enough attempts have been made to quantify some 
of these aspects. There is also a need to study the particle 
modification at low stress levels prior to particle crushing 
especially with reference to the liquefaction potential of 
sands. Available literature on this aspect is scanty. In 
the present study, an attempt has been made to investigate 
all these aspects. 

Three natural sands of differing mineral composition 
and an artificially prepared quartz sand were chosen . In 
addition, six more samples have also been prepared using 
these four sands through size fractionation, mineral 
separation and selective mixing . Seven samples of sands as 
obtained from RP I New York wore also used in the study . All 
these samples have been characterized In terms of their 
sphericity, shape factor, and angularity. The quantification 
of particle shape has been achieved using these three 
parameters . Compressibility tests wore carried out for the 
ten samples under differing relative density conditions and 
varying stress levels. A scheme for classification of sands 
on the basis of compressibility was evolved using the results 
from the present study and test data available from different 
sources in the literature. For the different classes of 
compressibility as proposed in the present work, studies on 
constrained modulus (N) have enabled the estimation of 
coefficients in the constrained modulus expression. In 
addition, the modifications in particle size and shape have 
been systematically monitored for specimens subjected to 
compression. The stress levels within which the particle 
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modification takes place for different sands , prior to 
crushing, have been establ i shed . Role of particle size , 
angularity, shape , relative density and mineral composition 
on compressibility, constrained modulus, modification and 
crushing has been established. Predictions of cone 
penetration resistance (q c ) for sands of different 
compressibility on the basis of results in the present study 
and available theoretical models ( Janbu 1974 ) are shown to be 
in complete agreement with the trends reported . 

The work contained in the thesis has been organized 
into six chapters. In chapter one , the topic has been 
introduced and the scope of the present study has been 
outlined. The available literature on the characterization 
of shape of particles, its influence on, stress-strain- 
dilatancy behaviour, compressibility characteristics, and 
mechanism of grain modification as well as crushing under 
different stress conditions has been reviewed in Chapter Two. 
The materials and methods involved are described in detail in 
Chapter Three . The four sands (Ganga , Calcareous , Kalpi and 
Standard samples) have been characterized for their 
mineralogical composition using X-ray diffraction. Particle 
shape analysis has been carried out using an image analyzer 
system. Details of this equipment and procedures of 
measurements are presented in this chapter. The various 
expressions used for the quantification of particle shape 
such as sphericity, shape factor and angularity are 
explained. Sample preparation and the test procedure for 
compression test has also been outlined. Chapter four deals 
with the experimental data concerning mineralogical 
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composition , physical character-istics , morphology, 

compressibility and crushing. While the Gangs Sand is 
composed of quartz (60-65%), feldspar (20-25%) and mica (8 to 
10%) with minor amount (upto about 3%) of chlorite and 
kaolonite, Kalpi sand has calcite upto around 18% with quartz 
and feldspar in equal proportion (40% each) together with 
minor quantity of mica ( upto 2%). The cal carceous sand has 
carbonate upto 95 % with subordinate quartz (upto 5%), 
Aragonite and calcite are the .main minerals . The standard 
sand is a pure quartz sand. Grain size distribution and 
morphological characteristics for all these sands have been 
presented . Data on compressibility for dry sands with loose 

and dense relative densities were obtained upto a confining 

2 

pressure of 176.32 Kg/cm and deformation characteristic® are 
presented. Grain size distributions for sand samples on 
completion of compression tests are contrasted with the 
original gradation. In Chapter Five, the test data obtained 
are interpreted in terms of various operative parameters and 
the details are presented on morphological clanai f ication and 
role of grain modification and crushing on compressibility 
and constrained modulus. The role of raineralogical 
composition in the behaviour of sand under stress has also 
been brought out. Interrelationships between various 
parameters used for the description of surface 
characteristics are established. It has been brought out 
that for the classification of sands, in addition to the 
usual tests (e^x. e , e m j; n » grain size distribution etc.), 
shape and angularity characteristics of sand grains should 
be evaluated and this has been accomplished using a 


also 
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technique proposed in the present study. I t is shown that 

based on grain size, shape, angularity, mineralogy, gradation 

and relative density natural sands experience particle 

degradation to a varying degree under stress. On the basis of 

particle gradation characteristics before and after 

compression, an expression for total potential for particle 

degradation as related to d^g has been established. It has 

also been brought out that there exists a critical stress 

level at which the grain modification of an alluvial sand is 

initiated. This stress is much less than the stress level 

required for actual grain crushing . I t has been indicated 

2 

that critical stress value is of the order of 8-9 kg/cm for 

7 . 

Ganga and calcareous sands whereas it is around 1 6 kg/cm f or 

Kalpi sand ; the standard ©and showing no modification at all 

2 

even upto stress level of 176 Kg/cm . It has been indicated 
that the compressibility of sands is strongly controlled 
among other factors by mineralogy and angularity. At a low 
stress level , the stress-strain behaviour of sands with 
significant amount of mica and or carbonates is mainly 
governed by grain inter-locking. At high stresses , the 
stress-strain behaviour is significantly different from that 
at low stress due to grain modification and crushing. 
Depending on mineral composition, particle size, shape, 
angularity and gradation, grouping of sands into three broad 
categories - sands of low, moderate and high compressibility, 
has been proposed on the basis of their comprens ibi 1 i ty 
characteristics. The power laws governing variation of 
constrained modulus with effective confining stress and 
relative density as reported in the literature, tacitly 
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assume that the sand remains unchanged during loading to high 
stress. It has been established in the present study that 
particle modification and crushing at high stress level 
(greater than the critical stress) influence the stress- 
strain behaviour of sand . Hence, these power laws are 
strictly applicable for stress levels less than the critical 
stress. The power laws of the type proposed by Belotti et 
al . (1985) between q c and effective stress as well as PI and 
effective stress, for a given relative density, has been 
observed to be strongly controlled by the compressibility of 
sands. The coefficients in these power laws have been 
quantitatively correlated with grain mineralogy and 
angularity . The predicted trends for q c are in agreement 
with those reported by Robertson et al . (1983). Based on 
these relationships for q c and H, the value of OC =q c /PI is 
shown to vary between 3 and 11 as generally reported. 

The conclusions arrived at on the basis of the present 
study are highl ight ed in Chapter Six. 



CHAPTER I 


INTRODUCTION 

1.1 G«neral 

It ia known that the mechanical behaviour of granular 
material ia controlled by the char act er iatca of ita conatl tuent 
particlea. Considerable work has been done on the material 
character ization in terms of particle size, gradation, and 
relative density as well as on the effect of these parameters on 
mechanical character i st i cs such as liquefaction potential, 

penetration resistance, shearing resistance, di latancy , 

compressibility and crushing. While the importance of grain 
shape, angularity and mineral composition in under standing the 
mechanical behaviour of granular material has been recognized 
(Hoi ube c and D'Appolonia, 1973? National Research Council, 19831, 
enough attempts have not been made to quantify these aspects and 
also to study their influence on the mechanical behaviour of 
sands . A detailed review of available literature as outlined in 
Chapter II , has clear ly brought out several aspects of granular 
materials that deserve an immediate and a greater attention. 
Particle shape and angularity, which have hitherto been specified 
in a descriptive manner , need to be quantified. No systemic 
details are available on grain morphology and mineralogy as 
related to their mechanical behaviour. No attempt seems to have 
been made by earlier workers to recognize the particle 

modification, in contrast to particle crushing, and also to 
systematically monitor the same at different stress levels. For 
sands of dominant carbonate composition, possibility of 
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modification and crushing at relatively lower stress levels 
exists compared to quartz sands and the available information is 
scarce. In fact, there is a need to classify sands in regard to 
their compressibility taking into consideration the size, shape 
angularity and mineral composition of constituent particles. 

In the present study, an affort has been made to develop 
quantitative procedures for describing particle shape and 
angularity. Experimental programmes were designed on a systematic 
basis to examine the influence of grain morphology and mineralogy 
nn the compressibility, modification, crushing and constrained 
modulus for a variety of sands . 

1.2 ORGANISATION OF THE THESIS 

In Chapter I, an introduction indicating the need of the 
present study is given with a discussion on the and the 
presentation of the material in the thesis. 

Available literature on the characterization of particle 
shape and its influence on mechanical behaviour of sands, stress- 
strain dilatancy behaviour of sands under varied confining 
pressures, compressibility of sands, and the mechanisms of grain 
modification and crushing under different stress conditions and 
their influence on engineering properties of sands has been 
reviewed in Chapter II. 

The characteristics of the starting materials and the test 
procedures adopted together with the details of the equipment 
involved are dealt in Chapter III, 



Results of experimental investigations carried out to 
characterize all the sand samples involved in the present study 
in terms of their mineralogy, physical and morphological 
characteristics form the theme of Chapter IV. 

In Chapter V, test data in terms of particle sphericity, 
shape factor and angularity have been interpreted. The 
classification of sands in terms of grain angularity has been 
presented. Results of compressibility, grain modification and 
crushing have been correlated with grain size, angularity, 
mineralogy and stress levels. The influence of grain angularity 
and mineralogy on cone penetration resistance ( q f > as well as 
constrained modulus <M) have been reported for different stress 
levels . 

The conclusions arrived at in the present study are 
highlighted in Chapter VI. Relevant references pertaining to the 
present study are given at the end of the thesis. 



CHAPTER 


I I 

LITERATURE REVIEW 


2.1 General 

The significant effect of confining pressure on behaviour of 
granular material has been recognized since the earlier days of 
Soil Mechanics . Yet, until recent years, the entire knowledge of 
mechanical behaviour of such material has been baaed on observa- 
tions limited to stress levels that could be achieved in 

conventional triaxial apparatus (mean normal stresses not excee- 
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ding 10 Kg/cm ). Even from a practical point of view, this 

situation has frequently led to unjustified extrapolations . For 

example, the stresses under deep foundations as in the case of 

very high earth and rockf ill dams may be of the order of 100 
2 

Kg/cm . In recent years intensive use of the principles of Soil 

and Rock Mechanics in geological and geophysical analyses, and 

Interest in evaluating response of geologic materials under 

blast loads , have extended the range of pressures of practical 

interest to well above 10,000 Kg/cm 2 .' 

The term high pressure obviously has only a relative 

meaning. For most physicists it would mean pressures exceeding 

1000 Kg/cra , whereas for practising geotechnical engineers it 

may mean any pressure that exceeds the capacity of ordinary 

2 

triaxial apparatus (about 10 Kg/ cm ) . The terminology for range 
of pressure used by different investigators, has been varied, 
depending on the nature of problems dealt with- Bishop (1965) and 
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Leps (1970) used the following terminology for pressure ranges 

dealing with response of soils under high dams: 

2 

Low Pressure 0.0 to 0.7 Kg/cm 

Medium Pressure 0.7 to 7.0 Kg/cm 2 

High Pressure 7.0 to 70.0 Kg/cm 2 

Very High Pressure > 70.0 Kg/cm 

whereas the pressure range terminology used by Vesic and Clough 
(1968) was: 

? 

Low Pressure 0.0 to 10 Kg/cm 

2 

Elevated Pressure 10 to 100 Kg/cm 

High Pressure 100 to 1000 Kg/cm 2 

Very High Pressure 1000 to 10000 Kg/cm 2 

Ultra High Pressure beyond 10000 Kg/cm 2 

In civil engineering, in some situations the range of 

7 

pressure can be as high as several hundred Kg/cm . Several 
investigations on the material properties of soils subjected to 
high pressures have been made by earlier workers using elevated 
cell pressure apparatus (Bridgman, -1918; Robert and De Souza, 
1958; Lee and Seed, 1967; Vesic and Clough, 1968; Agar et al , 
1987). Available data clearly point out to a need for 
differentiating the mechanical' properties of soil under high 
pressures from those determined at low confining pressures. 

At low pressure there is hardly any significant amount of 
particle modification or crushing and the sand particles being 
relatively free to move with respect to each other, the dilatancy 
effect can be quite significant. As the mean normal stress 
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increases , crushing becomes more pronounced and the dilatancy 

effect gradually disappears. Crushing appears to be most intense 

2 

in the elevated pressure range (10 - 100 Kg/cm ) until the so- 
called "break down stress” (Vesic and Clough 1969) is reached 
(this stress is defined as the stress needed to eliminate the 
entire effect of initial void ratio of sand). The main effects 
of grain crushing on properties of granular materials (Vesic and 
Clough 1969) are: 

1. Principal effective stress ratio at failure decreases with 

increasing confining pressure due to the suppression of 
dilatancy, but it does not tend towards a lower limiting 
constant value. 

2. Large volume reduction occurs. 

3 . Particle breakage and crushing produces fine material 
(powder) with consequent reduction of the permeability and 
dilatancy and increase the density , and resulting in further 
additional settlement due to reorientation of the grains. 

4. Much of the grain crushing occurs among the coarse and 
medium size particles in a specimen. 

5. The strength and stress-strain behaviour of an element of 

soil is affected greatly by the extent of crushing or 

particle breakage during loading and deformation. 

While most of the studies at high pressures have concentrated on 
total grain-crushing, the minor modification in the grain bound- 
ary characteristics have not received adequate attention. 
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In this study, the available literature has been reviewed 
under the following heads: 

1. Morphology of sands (physical characteristics of soil). 

2. Behaviour of sand at low and high confining pressures. 

3. Compressibility. 

4. Grain modification and crushing. 

2.2 MORPHOLOGY OF SANDS 
General 

Size classification of particles is commonly based on 
either the arithmetic or the geometric mean of the greatest, 
intermediate and least diameters. The values obtained in both 
cases are dependent upon two factors -the shape and the volume of 
the particle. Two particles, one irregular having varying 
diameters and the other a sphere with the three diameters of the 
same length, may both have the same "size" value as a result of 
the arithmetic or geometric mean of their diameters, although 
they may differ distinctly in respect to both shape and volume . 

Engineering classification of sands is generally based on 
particle size distribution and relative density without much 
regard to the particle morphology. For example, two sands at the 
same relative density and with similar gradation curves may have 
widely different maximum and minimum void ratio values depending 
upon the grain shape. Thus the maximum volume change potential 
(liquefaction under undrained cyclic loading) of a sand deposit 
will be significantly governed by the shape of the grains 



8 


(Burmister, 1938, 1962; Ki lbuzewski , 1963; Iehihara and Uatanabe, 

1976; Castro, 1969). 

2.2.1 Grain Shape 

Shape is a complex property of a grain and it is difficult 
to describe shape precisely except for grains that approximate 
regular geometrical shapes. The following aspects of shape are 
generally adopted by sedimetologists to characterize the shape of 
particle (Blatt et al , 1971): 

i Surface texture 
ii Roundnees 
iii Sphericity 

2. 2. 1.1 Surface Texture 

Surface texture is used to describe the type of surfaces of 
particle that are so small to affect the overal shape. Termino- 
logy commonly adopted for distinguishing surface textures in- 
cludes terras such as polished surface, greasy surface and 
frosting of the grain. 

2. 2. 1.2 Roundness 

Roundness refers to the sharpness of the corners and edges 
of the grain. It therefore refers to those aspects of the grain 
surface that are on a larger scale than those classed as surface 
texture, but still smaller than the overall dimension of the 
grain. Uentworth (1919) was the first to quantify roundness. He 
defined roundness as the ratio of radius of sharpest corner to 
the radius of largest inscribed circle. 
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Pent land (1922) was the first to express roundness of 
individual sand grains as a ratio of percentage area of the grain 
projection to that of a circle with equal diameter of the grain. 
The area of the grain projection was determined from camera 
lucida drawings. 

Cox (1927), on the basis of measurements made on the image 
of grains as projected on a screen, calculated the roundness or 
circularity as 


4AH 


K « 


( 2 . 1 ) 


where A - area measured by a planimeter 

p - perimeter measured by a map measurer 
K - roundness value. 

Tickell (1931) used ratio of the area of the projected grain to 
the area of the smallest circumscribed circle to express ’’round- 
ness” . 


Uadel 1 (1932) was the first to show that the commonly used 
terms "shape” and "roundness” were not synonymous, but really 
included two geomet er i cal ly distinct concepts. He defined round- 
ness as the ratio of the average of radii of corners of the grain 
image and the maximum radius of the inscribed circle. This 
definition is a modification of the Uentworth roundness. Russel 
and Taylor (1937) placed particles into classes based on 
comparison with photographs of type grains using five grade terras 
(Table 2.1). The class limits for the grade terms are based on 
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Table 2.1 Roundness Grades 


Grade Russel and Taylor Petti John 

terms Class limits Arithmatic Class limits Geometric 

Mid-Point Mid-Point 


Angular 

0.00 

to 

0.15 

0.075 

o . no 

to 

0.15 

0.125 

Subanftul ar 

0.15 

to 

0 . 30 

0.225 

0.15 

to 

0 .25 

0.200 

Subrounded 

0 . 30 

to 

0.50 

0 . 400 

0.25 

to 

0 .40 

0.315 

Rounded 

0 . 50 

to 

0 .70 

0.600 

0 . 40 

to 

0 . 60 

0 . 500 

Well Rounded 

0.70 

to 

1.00 

0.850 

0.60 

to 

1 .00 

0.800 
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the roundneos estimates using Uadel's method. These roundness 
grades are subsequently modified by Petti John ( 1949 ) using 
"geometric scale" and rounding off the class limits (Table 2 . 1 ). 
Pettijohn used Silhouette of the grains rather than their 
photographs (as used earlier) for the grouping of particles into 
classes . 

Powers (1953) has further modified the roundness scale by 
adding another descriptive term, "very angular" to Petti John's 
and Russel and Taylor’s set of roundness terms, which would 
provide an additional division of the lower classes. The inter- 
val from 0.12 to 1.00 has been subdivided into six intervals in 
such a way that the ratio of the upper to the lower limit of any 
interval is 0.70 (Table 2,2). The characteristic point of each 
interval is the geometric mean of the class limits. 

Folk (1955) suggested a logarithmic transformation of round- 
ness analogous to the <p scale for sire. In this scale, 
designated as the rho( p) scale, the p values extend from 0 . 0 
(perfectly angular) to 6.0 (perfectly rounded) with each unit 
corresponding to one class on the Powers' scale of roundneos. 

Pobklns and Folk (1970) suggested use of the radius of the 
largest corner divided by the radius of the largest inscribed 
circle as a measure of roundness. For precise measurement of 
small particles, the projected image of the grain must be enlar- 
ged to a standard size and the radii measured using a set of 
circles engraved on a plastic template. 



Table 2.2 Roundness Grades Powers (1953) 


Grade T or' ms 

Very AnfMi 1 ai 
Atif»u la? 

fjnb Anr u 1 ar 
,Gob Pounded 
R ou ml r*<j 


Class Interval Geometric flean 


n . 

i ?, 

n . 

1 7 

CJ . 

1 4 

n . 

1 7 

o . 

2 5 

f 1 . 

2 1 

o . 

2 5 

n . 

15 

n . 

in 

0 . 

,15 - 

n . 

, 4 9 

n , 

. 4 I 

n , 

<19 

n . 

7 0 

n . 

5 9 

0. 

.70 

* i , 

, 00 

o . 

. 04 


Well Rounded 



13 


Youd (1973) adopted the procedure ae suggested by Powers 
but with the modification to consider particles from each size 
fraction rather than the specimen as a unit. The procedure used 
for his study consisted of examining under a microscope a repre- 
sentative number of sand particles (at least 50) from each sieve 
fraction and visually assigning each particle to a category 
listed by Powers (Table 2.2), The average roundness value for the 
sieve fraction (Rj) was then calculated from the equation 

n 

R j « (I R x )/n (2.2) 

i = 1 

where R j is the roundness value assigned to the particle i, and 
n is the number of particle® examined. The average value of 
roundness for the sand specimen was then calculated from 

R = ( I P j R j )7100 (2.3) 

where P j is the percent (by weight) of particles retained in each 
sieve fraction. 

2. 2. 1.3 Sphericity 

The shapes of object® may be classified in a number of ways. 
The geometrician has defined a series of regular shapes such 
as cube, prism, sphere, cylinder, cone. The crystal lographer 
likewise uses a classification of solids bounded by plane 
surface. Neither system is adequate for the natural sediments. 
At best, pebble shapes only approximate the regular solids of the 
geometrician. Terms such as pr ismoldal , bipyramidal, pyramidal, 
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wedge-shaped and parallel-tabular may be used to describe 
particle shape (Uentvorth , 1936). Such a classification however, 
is at best a qualitative description only. Instead a single 
number index of shape would be desirable which is amenable to 
mathematical or graphic analysis and by means of which a shape 
distribution or frequency curve can be constructed. 

A quantitative or numerical index of shape requires some 
standard of reference. Uadel 1 (1932) was the first to choose 
sphere as a standard. A sphere is not only the limiting shape 
eventually assumed by many rock and mineral fragments upon pro- 
longed abrasion during transportation, but it has certain 
unique properties which make it an acceptable standard of refe- 
rence. 

Ideally, the property of sphericity might be defined as s/S 
(Uadell 1932) where "s” is the surface area of a sphere of the 
same volume as the fragment in question and S is the actual 
surface area of the object. For practical purposes this ratio is 
difficult to measure and the actual measurement is expressed in 
terms of a ratio of the volume of particle to the volume of its 
circumscribed sphere. The cube root of this ratio is called the 
sphericity of the particle. 

Uadell expressed the volume of the particle In terms of a 
sphere having the same volume, the diameter of the corresponding 
sphere is the "nominal particle diameter” (d). On this basis 
volume of particle is (JT/6)d 3 . The volume of the circumscribed 



sphere is in general based on the longest diameter (a) of the 
particle, so that the volume of this sphere is IT/ 6 a . Therefore, 
the Uadell's sphericity reduces to the ratio of the nominal 
diameter of the particle to its longest diameter (d/a). 

As a further improvement,, Krumbein (1941) used the concept 
of intercept sphericity, wherein the volume of particle is expre- 
ssed in terms of a triaxial ellipsoid having the three diameters, 
a, b, c with a>b>c. The volume of such an ellipsoid is JT abc/6 
and sphericity is given as 

3 p bc~ 

= J -J— (2-4) 

& 

Sneed and Folk (1958) have objected to the use of sphericity 
as defined by Uadell on the ground that although it accurately 
measures the degree to which a particle approaches a sphere, it 
does not correctly express the dynamic behaviour of the particle 
in a fluid. This is because particles tend to orient themselves 
with their maximum projection area normal to the flow and this is 
true both for particles settling in a fluid and for particles 
resting on the bottom. In either case, the fluid drag is almost 
the same as the force on a sphere with the same projection area 
as the particle, whereas the gravity force is the same as on a 
sphere with the same volume of the particle. Defining the "maxi- 
mum projection sphericity" as a ratio of the diameters of these 
two spheres, Sneed and Folk (1958) expressed sphericity as: 



16 


Based on a measure o £ flatness and elongation ratios (p and 
q respectively), Aschenbr enner (1956) and Lees (1964), proposed 
the following relationship for particle sphericity: 

12.8(p 2 q)^ 3 

V = o (2.6) 

1 +p( 1 + q ) + 6 CCl + p 4 Cl+q^)] 1/Z 

2.2.2 Classification of Shape Type 

"Definition of roundness and sphericity of a particle does 
not define the shape uniquely, even for the case of ideal 
ellipsoids” (Blatt et al 1971). Several other indices have been 
proposed and used by sedlmentologists to define other aspects of 
shape . 

Zingg (1935) proposed the use of two shape indices, dj/d^ 
and d s /dj, where d L ,dj and d s represent the longest, intermediate 
and the shortest particle dimensions. These indices are used to 
construct a diagram, commonly called the Zingg diagram (Fig. 
2.1), which is divided into four quadrant® defining four main 
shape classes - tabular (oblate spheroid or disk shape), equant , 
bladed and prolate spheroid (rod or roller shaped). The 
boundaries for each region are indicated in Table 2.3. 

Equant particles have highest sphericity but oblate, bladed 
and prolate particles may all have comparable sphericity values. 

Krumbein (1941) superimposed the Uadell sphericity fields 
on the Zingg diagram (Fig. 2.1). 

Heywood (1947) defined the two indices B/T and L/B as 
flakiness and elongation respectively, where L, B and T are 
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Table 2.3 

Quadrant 

I 

I 1 
I I I 
IV 



Zingg Classification of Particle Shape 


d I /d L 

d s /d j 

S h a p e 

2/3 

< 2/1 

Ob 1 ad o (dinks) 

2/3 

> 2/3 ' 

* 

Eq u an t. 

2/3 

> 2/3 

R 1 ad od 

2/3 

< 2/3 

Prolate (roller. 



Rod - I ike) 
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analogous to d^, dj and dg) . He further pointed out that if 
approximately B=L=T the particle is said to be bulky. The 
particle is said to be flaky, platey, scale-like if L=B>T and rod 
or needle-like when L> B=T. 

Sneed and Folk (1958) presented triangular diagram combining 
the two indices dg/d^ and (d^-dj )/(d^-dg) . The diagram is 
divided into a number of shape fields and line of equal maximum 
projection sphericity may be drawn on it. 

Lees (1964) defined the flatness (p) , elongation (q) and 
shape factor respectively as 

c 

p = — (2.7) 

b 

to* 

q = — (2.8) 

a 

P c a 

F • = — (2.9) 

q b 2 

where a > b > c (analogous to dg , dj , d L ) - 

It may be noted that for a given value of F , the grain shape can 
vary considerably depending upon either flatness or elongation 
ratio. For example when F=l, c/b=b/e = k (say) the particles 
will be equldimensional only for variation of k as 2/3 ^ k < 
1.0. Particles with larger elongation ratio (compared to 

I 

flatness ratio) are disc or blade like in shape. Uhereas grains 
having flatness ratio greater than elongation ratio 


are 
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dominantly rod or roller like in shape. Accordingly to Blatt et 
al (1971) equation 2.5 for sphericity is similar to the expre- 
ssion given below for shape factor (SF) as used earlier by 
engineers to express the main effect of shape on settling 
velocity: 



2.2.3. Effect of Morphology of Sand on Its Engineering 
Properties 

Burmister (1938, 1962) observed that the range of particle 

sizes is the most important factor controlling maximum and mini- 
mum densities. The density increases with the increase in the 

¥ 

range of particle sizes. Other factors he considered important 
were the type of gradation curve and particle shape. 

Smith (1965) found a 3 percent decrease in n ro ax“ n min 
(maximum porosity range) as grains increased in roundness from 
very angular (crushed basalt) to rounded (Erith sand), which is 
consistent with the trend observed by El Sobby (1964) who 
reported a 5.7 percent difference between the values of porosity 
interval for crushed feldspar and glass ballot ini. 



Dickin (1973) studied the influence of grain shape and size 
upon the limiting porosities of quartz sands (5 different sands) 
and glass ballotini and concluded that both the maximum and 
minimum porosities decrease with increasing sphericity for the 
fine sands chosen for their study. In the case of natural quartz 
sands, limiting porosities decrease with increasing grain size 



for any one sand as a result of the inherent increase in 
sphericity with grain size. 

Ishihara and Uatanabe (1976) proposed use of the term volume 
decrease potential, ®niax " e min (maximum void ratio range) for 
assessing liquefaction potential of sands. They suggested prog- 
ressive decrease in ( © roax - e m j[ n ) values with increasing average 
particle size ( d$g , mm). It may be pointed out that in case of 
sedimentary sand deposits, larger particle size would be approa- 
ching high values of roundneao due to abrasion during transpor- 
tation as compared to finer sand grains which are likely to be 
angular in shape f or the same distance of transportation. Data 
for Fuji river sand given by Ishihara and Uatanabe supports this 
conclusion (Dicken 1973) . Thus it is evident that volume 
decrease potential is also governed by the particle shape. Youd 
(1973) observed that the particle size has very little effect on 

®max an< * ®min‘ 

The ( n max ” n min^ values for dg 0 = 0.11 mm is 15 % for Fuji 
River sand and 8% for glass beads , giving a decrease of 7% as 
the roundness of the particle increases . This trend is In 
conformity with the results of El Sohby (1964) and Smith (1965). 

Charles and Uatt (1980) have carried out large scale 
drained triaxial compression tests on rockfill and found that the 
maximum principal stress ratio at failure is related to the rate 
of dilation which is itself a function of particle shape. 


strength and relative density. 
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Usesugi and Kiahida (1986) examined the effect of signifi- 
cant factor on the frictional coefficient using three different 
oands °nd found that the particle shape, mineral ingredients and 
surface roughness of the particles have considerable influence 
on the frictional coefficient. 

2.3 BEHAVIOUR OF SANDS AT LOU AND HIGH CONFINING PRESSURES 

Considerable amount of research work has been published on 
the behaviour of granular materials , based on oedometer and 
triaxial shear tests . Many of such studies have been limited to 
low to elevated cell pressure range . Research work carried out 
under high confining pressures was restricted mainly to small 
size specimens (typical volume being equal to a few cubic centi- 
meters or a fraction of a cubic inch). Available information on 
large size particulate materials tested at high pressures is very- 
limited. In this review, the behaviour of sands under cyclic 
loading is not considered. 

On* of the earliest investigations concerning the 
compressibility of sand was by Terzaghi (1925). On the basis of a 
number of confining pressure tests on sand and clay soils, 
Terzaghi observed that the compressibility of both types of soil 
was similar. However, the volume changes in sand were not as 
great as those which occurred In clay soil when tested under the 
same magnitude of load. Nevertheless, the sand was found compre- 
ssible even at pressures of the order of 3 Kg/cm . An extension 
of this work Is quoted by Terzaghi and Peck (1948). Several sands 
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and sand mica mixtures were nub j ect ed to one dimensional 
compr onelon, and the particle crushing that was observed was 
suggested to be responsible for the compressibility of the sand 
at high pressures. 

Bridgman <1918) was able to compress an angular sand 

2 

(grounded quartz ) to 30000 Kg/cm , reducing its porosity to about 
2 per cent . 

Robert and De-Souza (1958) tested different sands and 
observed that at moderately low pressure, crushing and 
consequent compression are more in angular sands compared to the 
rounded ones . However , at very high pressure little difference 
wan observed between the behaviour of angular and rounded sands . 

Hardin, Borg, Fredman and Higgs (1960) and Griggs, Turner 

and Heard (1960) conducted triaxial tests on small samples of 

2 

sand at very high pressure ( upto 5000 Kg/cm ) combined with high 
temperature and reported the ductile failure at large strains in 
most samples with no visible effect of temperature upto 600°C. 

Hirschfeld and Foul as ( 1 9 6 3) on the basis of tests on silt 
and sand at different confining pressures concluded that while 
there is a net increase in volume during shear at low confining 
pressure, a net decrease in volume results at high confining 
pressure . 

Vesic and Barksdale (1963) investigated the mechanical 

behaviour of sand in the high pressure range on samples of dense 

2 

sand (2,2 cm diameter) at all pressures upto 630 Kg/cm and found 
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that the curvature of the flohr envelope of the sand tested was 
restricted to an initial pressure range. They showed by a limited 
number of volume change measurement in shear at constant mean 
normal stress that the effect of dilatancy of the material and 
the curvature of its strength envelope disappeared at high 
pressure . 

Lee and Farhoomand (1967) tested crushed granitic material 

o 

(angular to subangu 1 ar ) upto confining pressures of 140 Kg/ cm 

under isotropic and anisotropic conditions and reported increase 

in compr eonib 1 1 i ty with increase in confining pressures. 

Lee and Seed (1967), In the case of an intially loose 

specimen of fine uniform Sacramento river sand subjected to 

2 

conf ining pressure greater than 100 Kg/cm , were able to obtain a 
final void ratio that was less than the minimum void ratio 
obtained under normal tests conditions in the laboratory. 

Veoic and Clough (1968) studied the behaviour of granular 

material at high pressures extending beyond the mean normal 

2 

stresses of 1000 Kg/cm . Data reported by them on Chattachoochee 
river sand revealed the existance of a mean normal stress beyond 
which the curvature of strength envelope for all Initial void 
ratio vanishes, and the shear strength of sand is not affected by 
its initial void ratio. This stress is designated as the crush- 
ing stress or breakdown stress as it represents the stress level 
at which all dilatancy effects disappear and particle crushing 
becomes the only mechanism, in addition to shearing displacement 
due to simple slip along the slip planes. It appears that this 
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stress is affected by gradation and shape of particles of the 
granular material as well as the mineral composition. 

Lee and Seed (1970) conduct ed drained and undrained tria- 


xial 

teat on 

fine, uniform clean graded Sacramento river 

sand 

^ e max 

= 1.03, 

e min 

0.61) with relative density 

varying 

from 

loose 

to dense 

(dense 

having relative density of 

78%), and 

con- 

eluded 

that , 

for the 

granular soils the volume 

changes 

which 

occur 

during 

consolidation appear to depend 

only on 

the 

magn 1 t 

ud e of the major 

principle stress and to be 

entirely 

inde- 


pend ent of the minor principal stress during consolidation. 

Miura and Yamanouchi (1973) conducted isotropic consolida- 
tion tests on "Toyoura Standard Sand” for two relative densities 

of 24 percent and 90 percent. Tests were performed upto isotropic 

2 2 

pressure of 500 Kg/cm . At 300 Kg/cm' isotropic pressure, two 
samples showed 21 and 10 percent volumetric compression respecti- 
vely. They also concluded that void ratio of the two samples 

2 

approaches a certain value at a pressure of 500 Kg/cm'. 

Agar et al (1987) investigated the strength and stress 
deformation behaviour of Athabasca oil sand at elevated 
temperature (25 - 200 °C) and confining pressure of 800 Kg/cm 2 
using triaxial compression test and concluded that the shear 
strength of dense sand increases non-linearly with effective 
confining pressure. 

2.4 COMPRESSIBILITY 

The compressibility of granular materials has been studied 
both in anisotropic (i.e. in oedometer) and in isotropic (I.e. 



in trl axial ) conditions. 
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From these studies it is concluded that 
the compressibility of granular media is a consequence of complex 
phenomena resulting from displacements between particles combined 
with particle breakage. The deformation of granular soils under 
loading is a result of three main mechanisms (Loung 1982). 

1. Consolidation mechanism - Change in shape and 
compressibility of the part icl e assemblies. 

2 . Distortion mechanism - Bending of the flat particles , 
sliding and rolling of rounded grains. 

3 . Attrition mechanism - Particle crushing and breakage, 
modification of the soil fabric and packing . 

Stress-strain behaviour of sands as studied in oedometer 
test ( H endron , 1963; Roberts, 1964) has been considered by 
Lambe and Whitman (1969) to be in three stages as detailed below: 

1. Locking behaviour during which sand gets stlffer and 

stiffer as the level of stress increases ”upto 150 
? 

Kg/cm', and finally, a stage is reached in which already 
dense arrays (of particles) are being squeezed more 
tightly together as contact points crush, thus allowing 
a little more sliding”. 

2. Yield behaviour as a "result of fracturing of Individual 
sand particles which permits large relative motions 
between particles” at stresses equal to and greater 
than 150 Kg/cm 2 . Considerable particle degradation 
takes place during this stage. 
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3 . Fracturing and crushing leading to tighter packing of 
particles and increase of number of particle contacts 
and "thus the sand once again becomes stiffer and 
stiffer as the stress increases still further”. 

These distinct stages of compression were observed for 
medium to coarse, well rounded quartz sand. However, actual 
stress level at which compression behaviour of sand changes from 
locking to yielding depends on many factors such as mineralogy, 
particle size, shape, roughness, angularity, gradation and rela- 
tive density, composition, cementation, average stress level , 
prentr ens ing , stress path and ageing. 

Effect of the different factors such as particle size, 
shape, roughness, angularity, relative density and mineralogy on 
the compressibility has been studied earlier. In the succeeding 
paragraphs, a brief reivew of the available literature pertinent 
to the present study Is made. 

One of the earliest investigations concerning the compress- 
ibility of sand was by Terzaghi (1925). He reported the increase 
in compression with confining pressure. Many such experiments 
have since been performed by investigators interested in 
compressibility of different granular materials. 

Burmister (1948) reported that the deformation increases 
with increasing void ratio and that different soils with 0% 
relative density are not equally compressible. 
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Schultz and Housa (1961) reported that compress! bl 1 ity o£ a 
sand is considerably increased when a small quantity of cohesive 
soil is mixed with the same. 

The study of compressibility on Leighton Buzzard sand by 
Kolbuszeweski and Fredrick ( 1963 ) has shown slight Increase in 
compressibility as roughness of particle increases when samples 
are dense, and a large increase in compressibility as roundness 
of particle decreases for loose samples. 

De Boer (1963) reported that the compressibility and crush- 
ing of particles are accelerated by the addit ion of water and 
increase with Increase In angularity and grain size of uniform 
soil, and decrease in strength of individual grains . 

K jaernsl i and Sande ( 1963 ) showed that for a given material, 
compressibility is low for well-graded material In comparison to 
uniformly graded material. For example, the compressibiity of 
crushed angular particles of syenite Increased from 4.75 to 16 
percent with a decrease in the coefficient of uniformity from 2.5 
to 1.4. They also observed that the material with weakest 
particles show the greatest amount of compression and particle 
breakage. 

Marshal (1967) performed drained triaxial test on coarse 
gravel and broken rock upto 20 cm in diameter, with confining 
pressures upto 25 Kg/cm and shoved that compressibility is 
related to the nature of the soil. The material with weakest 
particles showed the greatest amount of compressibility. 
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Koerner (1968) reported that deformation Increases with 
decreasing coefficient of uniformity and decreases with increas- 
ing roughness, 

El Sohby (1969) studied the effect of st r ess on resulting 
axial deformation and volume change for three different 
materials- Silver Sand ( Dg q = .25. Djq = .14), Glass Bollotini 
( D 50 = 0.1 mm), fine sand and found that the total deformation of 
a mass of sand subjected to a constant stress ratio can be 
divided into two main components i.e. sliding deformation and 
elastic deformation. He further reported that deformation 
increases with increasing void ratio. 

Koerner (1970) conducted tests on quartz, feldspar and 
calcite and reported more compressibility and degradation for 
uniform soils in comparison to the well graded soil having same 
maximum particle size. 

Charles Moore ( 1971 ) performed two series of testa on 
Chattahoochee river sand (with addition of mica upto 1, 5 , 10, 
50 and 100 percent) and residual silt from Atlanta. Results of 
testa on Chattahoochee river sand conducted with relative density 
values of 25,30 and 57 percent indicated that the compressibility 
Increases appreciably with Increasing mica content (5 percent of 
mica will double the relative compressibility). 

Bordon and Proctor (1971) performed a series of tests on 

river Ui 1 land sand using triaxial compression test with cell 

o 

pressure of 703 Kg/cm and showed that the dense sample is 
highly dilatant with insignificant crushing whereas loose sample 
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©hows no di latancy but has significant amount of particle 
crushing . 


El-Sohby and Andrews (1972) have reported results from an 

experimental study on variety of granular materials using confin- 

2 

ing pressure upto 6.33 Kg/cm only, and showed that the compress- 
ibility of granular material is mainly governed by the properties 
of its individual particles and geometry of its internal struc- 
ture, i . e . the deformation increases with increasing void ratio, 
decreasing coefficient of uniformity and decreasing mineral 


strength. 

Ramamurthy and Kanitkar (1972) conducted isotropic 

consolidation test on four fractions of Badarpur sand . 

Larger volume changes were bbserved for the coarsest 

material tested ( d^Q = 1.7 mm). The volumetric compression 

during isotropic consolidation in the pressure range of 8.5 
2 

to 84 Kg/cm varied from 1.00 to 4.2 percent for mean 
particle size of 0.096 mm, 1.9 to 7.00 percent for size of 
0.5 mm, 2.4 to 11.3 percent for size of 1.00 mm, and 3 to 
15.5 percent for size of 1.7 mm. 

y 


D Appolonia and Holubec (1973) showed that the particle 
shape has minor influence on the minimum void ratio, but has a 
major Influence on the maximum void ratio. Since the maximum 
void ratio increases considerably with angularity than in the 
case of the minimum void ratio, the void ratio difference 
increases markedly with increase in angularity. They suggested 
that the def orroabi 1 i ty of sands having the same gradation and 
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relative density increases with increasing angularity. Also, 
strength, as measured by friction angle, increases with increa- 
sing angularity. 

Cornforth (1974) reported that deformation of granular mate- 
rial increases with increasing void ratio. 

Haruyama (1977) studied the deformation characteristics of 
highly compressible angular, uniformly graded volcanic deposit 
(Shi rasu sand) with different relative densities, using triaxial 
test and the results were compared with uniformly graded Toyaura 
standard sand. It was revealed that Shirasu sand showed higher 
compressibility as much as three times compared to that for 
Toyoura standard sand . 

Oberol (1983) performed oedomet er tests on dune sand, Chao- 

2 

parya river sand and micaceous sand upto 30 Kg/cm and reported 
volumetric strains of 4.5, 7.8 and 34 percent respectively. 

A1 Hussiani (1983) conducted drained triaxial compression 

o 

and plane strain tests at 31 Kg/cm using crushed basalt and 
found that compressibility decreases as the uniformity coeffi- 
cient Increases. 

Clayton et al (1985) have reviewed factors affecting compre- 
ssibility of sands. Following Daramala (1978) they have grouped 
these factors Into two categories - the material-dependent 
factors and the stress-dependent factors. According to them all 
the factors In both categories, angularities of particles and 
void ratio have the maximum effect on compressibility of sands. 
It has been shown that the compressibility of an uniform fine 
sand is strongly dependent on its previous stress path. 
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Kapoor (1985) tested four different sands with different 
relative densities, using oedometer and triaxial compression test 
and showed that dense samples are highly dilatant with insigni- 
ficant amount of particle breakage whereas loose samples show no 
dilatancy but have significant amount of particle breakage. 

Beloti et al (1985) observed the variation of tangent 
modulus with relative density and normal stress and suggested the 
following equation 


H 


m o p o 



exp (m 2 x D R ) 


( 2 . 10 ) 


where 

2 

P c = reference stress (1 Kg/cm'), 

O' = effective normal stress, and 

m c , mj , m 2 are experimental constants. 

Yudhbir and Rahim (1987) have reviewed the factors controll- 
ing compressibility. Based on available data, it was shown that 
"in general the compressibility is high when the sand is loose, 
the grains are angular , the coefficient of uniformity, Cu , is 
low, the average particle size is small, particle surface is 
smooth, the grain mineral strength is smal 1 , interparticle cemen- 
tation is weak and the ageing effects are Insignificant. Stress 
history, stress level and stress path have significant influence 
on compressibility of sands". It was also concluded that the 
angular well-graded sands attain closer packing of particles 
compared to poorly graded and rounded sands for same height of 
fall in a raining technique of sample preparation. Also angular 
well graded sands are more compressible than poorly graded 
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Band. It was further pointed out that in case of angular sands, 
the grain boundaries are modified at a stress level less than the 
critical stress, and at stress levels greater than the critical 
stress, sands behave like normally consolidated clay, exhibi- 
ting compressibility. 

Phool Chand (1988) investigated the response of three sands, 
namely standard rounded quartz sand, angular micaceous Ganga sand 
and angular Kalpi sand (having quartz , feldspars and carbonates 
as major minerals ) , under one-dimensional compression and 
unloading. The modulus during unloading was shown to follow the 
relationships as proposed by Janbu (1985): 

M = H aw a' (2.11) 

d a' 

n = ( 2 . 12 ) 

de 

where 

It 0W is a modulus number, 


and 0 ' = effective vertical stress. 

He observed that for standard sand practically one value of tt BW 
was obtained for samples with two different relative densities 
(25 % and 55 %), whereas for angular Kalpi and Ganga sand 
distinct values of H 0U were obtained for each density. Also M sw 
was shown to be independent of the past maximum pressure in case 
of rounded standard sand whereas these values varied with past 
maximum pressure (higher value for lower past maximum pressure) 
in case of angular Ganga and Kalpi sands. 
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angular flaky micaceous sand. This ratio J1 sw /mo) appears to be 
independent of relative density (he studied only two relative 
densities) . 

2.5 PARTICLE CRUSHING 

When soil grains are hard , tough and well - rounded , they 
may take up reasonably high stresses without crushing . Angular 
particles of freshly quarried material undergo fragmentation 
under ordinary pressures due to breakdown of sharp angularities 
(grain modification). The decomposed particles split under 
ordinary pressures In addition to breakdown of angularities. 

The Hertz theory of contact stresses between elastic sphe- 
rical particle® can be used to estimate the radius of zone of 
particle contact and the magnitude and distribution of contact 
stress in this zone. The load per particle will obviously depend 

on the mode of particle packing and the particle shape (Harsal, 

* 

1963). Harsal s theory brings out the influence of void ratio 
(packing), particle shape factor (sphericity), particle diameter 
and the average number of contacts per particle. 

Grivas and Harr( 1980) have investigated particle contacts in 
discrete materials and suggested that, for all practical 
purposes, the product of number of contacts per particle, and the 
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porosity (n) is a constant equal to three. The response of sand 
particles at their contacts would be governed by the magnitude of 
both normal and tangential contact stress. 

Hindi in and Deresi ewks ( 1953) after introducing tangential 
stress component suggested that higher stress is developed in the 
anulus of the slip on the periphery of contact area. Smith and 
Liu (1953) have investigated the distribution of stresses due to 
normal and tangential forces on the contact between rolling 
bodies. They have shown that a rolling contact leaves behind it a 
tensile stress , while introducing a compressive stress ahead of 
i t . 

In addition to the existence of normal and tangential forces 
at the contacts of particles during reorientation, the contacts 
may also be subjected to twist due to the movement of 
neighbouring particle. Studies of Hentengi and Hc-Donald (1958) 
and others revealed that spheres, when subjected to combined 
effect of pressure and twist, experience severe stresses on the 
contact surface resulting In rapid destruction and failure of 
mat er ial . 

In the presence of sharp angularity and asperities, the 
stress concentrations are likely to be higher than that the 
points could possibly bear. The soil grains experience crushing 
under direct loading or shearing of angularities and asperities 
under shear stresses. Under Increasing stresses, the asperities 
are likely to undergo progressive crushing until the contact area 
Is sufficiently increased to withstand the stresses without any 
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further crushing, or the interference of the sheared material may 
restrict the development of resistance between the contacts. 

2.5.1 Types of Crushing 

The different types of crushing (Ramamurthy , 1966) occuring 

in granular soils during shear may be classified as follows: 

1. Crushing or shearing of asperities on the contact 
surface. 

2 . Splitting of soil grains. 

3 . Breakdown of sharp asperities. 

4. Breakdown of asperities due to the stress reversal. 

2. 5. 1.1 Shearing of Asperities 

The frictional resistance between soil grains depend® on 
their surface properties and mineral constituents. Locally 
stronger asperities act as indentors on the contact surface and 
develop resistance to plastic deformation of material. Uhen 
shear stresses are applied, these asperities plough through or 
scratch the surface indentation and increase the resistance. An 
asperity can scratch a surface if the surface hardness of an 
asperity is at least equal to 1.2 times the hardness of the 
surface to be scratched (Tabor, 1956 ). The resistance developed 
during scratching and indentation is essentially the plastic 
property of a material. Bridgeman ( 1953 ) using a high pressure 
equipment has demonstrated that under hydrostatic pressures, a 
material which is brittle under ordinary conditions can be 
prevented from fracture and made to deform plastically. 
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The proportion of the crushed material is very small, since 
it involves only the asperities In regions of high shear 
stresses. The degree of crushing varies according to the 
magnitude of the stress and movement of particle. Even though 
the magnitude of crushing may be of the order of a few per cent 
of the entire sample. Its effect on the strength will be dispro- 
portionately large. 


2. 5. 1.2 Splitting of Soil Grains 

Splitting of soil grains takes place under high compre- 
ss i v o stresses. Splitting is encountered during field/ laboratory 
compaction of granular soils, particularly in clayey sand and 
clayey gravel type of soils. Its extent depends upon the shape, 
size and strength of individual grains and also on contact pres- 
sures. Splitting may be limited to the contact region or may 
extend across whole grain. In the preparation of laboratory 
specimens of coarse grained soil, large local voids tend to 
develop during deposition. If these specimens are compressed 
statically in a mould to achieve the desired density (especially 
if the movement occurs at one end only), the particle will split 
at the contact. 

The fragments formed due to grain splitting are coarser than 
those formed due to the shearing or crushing of the surface 
asperities. The broken fragments are generally angular and may 
offer greater surface frictional resistance on the newly formed 
surfaces. The magnitude of crushing in static compaction varies 
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with the difference between the density achieved during deposi- 
tion and that after compaction. And, therefore, a very open 
initial soil structure experiences greater crushing while stati- 
cally compacting the specimen to a given density. 

2 . 5 . 1.3 Breakdown of Sharp Asperities 

In quarried materials and residual soils, sharp corners of 
particles break during compaction or shear, producing extra blunt 
corners, though reducing the overall angularities. The newly 
formed surface and corners are relatively stronger. The asperi- 
ties on the newly formed surfaces are stronger at low pressures 
and may not exhibit crushing, but they may behave differently 
under different stress ranges . 

2.5. 1.4 Breakdown of Asperities due to the Stress Reversals 
This category of crushing occurs when hard soil grains are 

subjected to stress release at their contact points. Under high 
contact stresses , the angularities and asperities of particles 
experience elastic and plastic deformations. When shear strains 
are induced, some of the initial contacts may be unloaded; when 
the load is transferred to new contacts such stress reversals 
results in fragmentation. Different kinds of crushing are 
indicated in Fig. 2.2. 

2.5.2 Quantitative Expressions for Degree of Crushing 
The degree of crushing has been expressed quantitatively by 
several Investigators. 



s 


I N 

(a) Crushing of asperities 
(introduces lubrication 
effect ) 



(c) Breaking of angularities 
(introduces more corners 
and stronger surfaces) 


(b) Shearing of asperities 
(introduces smoother 
surface) 



(d) Splitting of particle 
(introduces more frag- 
mentation and stronger 
surfaces) 


Fig, 2.2 Different types of crushing ( Ramamurthy , 1966) 
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Leslie (1963) used the percent passing the sieve on which 
the original material was 100 3 retained. 

flarsal (1965) has selected as a measure of grain breakage a 
parameter equal to the sum of the positive Ay^, where Auk * a 
the difference between the percentages of the total sample con- 
tained in each grain size fraction before and after the test. 
Hence B ft is the percentage by weight of the solid phase that has 
undergone breakage . The parameter B g ( expressed as a ratio) 
multiplied by the concentration of solid (q), is the volume of 
the broken particles per unit total volume 





1 + e 


(2.13) 


Lee and Farhoomand (1967) studied crushing of material under 
anisotropic compression and defined a parameter designated as 
relative crushing which is equal to where Dj ^ is the 
diameter through which 15 % of particles of the original material 
pass and Dj^ a is the diameter through which 15 % of the particle 
pass after being subjected to anisotropic compression. 

Ramamurthy and Lai (1970) have expressed quantitatively the 
degree of crushing as the area between the grading curves of a 
sand obtained before and after it has been subjected to shear. 

Lowrison (1974) suggested that the amount of crushing can be 
estimated by evaluating the increase in surface area of the 
particles after crushing. 
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Leslie ( 1975) modified his earlier method and used the 
increase in percent passing the sieve on which 90% of the 
original sample was retaind . 

Datta, Gulhatl and Rao (1980) expressed the magnitude of 

crushing in terms of a crushing coefficient, C c defined as: 

Percentage of particle of sand after being 
subjected to stress finer than D 10 of the 
original sand 

C c = (2.14) 

Percentage of particle of the original sand 
finer than Djq of the original sand. 

Hardin (1985) has critcal ly reviewed the factors affecting 
crushing of sands. He has proposed a numerical measure of degree 
of crushing in terms of "potential breakage”, "total breakage” 
and "relative breakage”, and developed equations that can be used 
to evaluate the above mentioned terms for a given soil subjected 
to a specified loading. The potential for breakage of a particle 
of a given size D , may be represented by : 


1 f 

f leant 

element 

divided 


D in mm 

= log 10 for D l 0.074 mm. (2.15) 

0.074 mm 

= 0 for D ^ 0.074 mm (because extremely large stresses 

will be required to produce crushing of particles 
with D < 0.074 mm). 

b p represents the potential for breakage that la signl 
to soil behaviour for a. given size fraction, df, in an 
of soil, where df is a differential of "percent passing” 
by 100, then the potential breakage is given by. 
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1 





B P = J 

o 

bp df 

( 2 . 16 ) 

The 

amount of 

crushing 

that is 

signi f i cant 

to soil behaviour 

given 

by 

1 

B t = / 
o 

(bp 0 ~ 

bpj) df 

(2.17) 


where 

b P 0 

the 

original value of bp 



bpj 

the 

values of 

bp after loading 


The relative breakage i o given by 

B t 

B r = — — (2.18) 

B P 

2.5.3 Particle Crushing and its Effect on Engineering Properties 
of Granular Soil 

The strength and stress-strain behaviour of an element of 
granular soil is affected greatly by the degree to which crush- 
ing or particle breakage takes place during, loading and defor- 
mation. For the type of deformation that primarily produces 
volume change, such as one dimensional strain or isotropic com- 
pression, particle breakage adds to the reduction in volume. For 
the type of deformation where particles are moving past or around 
one another, such as triaxial compression or simple shear, crush- 
ing at sliding contacts or breakage of particles decreases the 
rate of dilation corresponding to a given principal stress 
ratio. Sliding between particles is usually present at all 
stress levels. 

The stress level at which crushing and fracturing of 
particle becomes significant is called the "critical stress”. 
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Thin critical stress Is smallest when the particle size is 

large, the soil Is loose, the particles are angular, the strength 

of Individual mineral particles is low, and the soil has a 

uniform gradation (Lambe and Whitman , 1969). 

Terzaghi (1925) reported an increase of 4.6% fines (dust 

content) of a loose quartz sand (made by grinding grain size 0.7 

to 0.25 mm) loaded in consl idometer under a stress of 50 Kg/cm 2 . 

He cited this increase in dust cont ent as a proof that the 

compaction of structure (of sand grains under stress) was 

accompanied by breaking of corners and ' edges (grain modification). 

H 1 1 f ( 19 57 ) observed that the residual and coarse-grained 

soils exhibited field dry density values higher than laboratory 

7 

densities by about 32.04 to 48.05 Kg/ cm ‘ due to breakdown of 
soil grains. Robert and De Souza (1958), Rowel (1955), Housa 
(1961), Vesic and Borksdal e (1963), Doe Beer (1965) have des- 
cribed the significant amount of particle crushing and increase 
in compressibility of natural sands at pressures beyond about 150 
Kg/cm 2 . 

Hardin (1953), Borg et al (1960), and Griggs et al (1960) 

have carried out a certain number of triaxial tests on small 

2 

samples of sand at very high pressure (upto 5000 Kg/cm ) combined 
with high temperature. Elaborate data were presented on struc- 
tural breakdown of the material, as well as on the fracture 
phenomenon of Individual sand particles. 

Leslie (1963) reported on testing of large samples of 

several poorly-graded to well-graded gravelly materials at all 

2 

pressure upto 47.5 Kg/cm . Harsal (1967) tested large samples of 
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coarse gravel and rockf ill (45.72 cm diameter) at confining 
pressure upto 25 Kg/cm 2 . 

Bishop, Uebb and Skinner (1965) tested loose sand samples 
by applying all pressures upto 10 Kg/cm . Lee and Fahroomand 
(1967) performed triaxial tests on six different sand (fine and 
gravel size) samples under a confining pressure of 140 Kg/cm . 
All these investigators reported the structural breakdown of the 
tested material and also concluded that uniformly-graded soil 
undergoes more particle crushing compared to well-graded one. 

Hall and Gordon (1963) performed drained triaxial compre- 
ssion test on gravelly soil and observed considerable amount of 
crushing especially during the shearing stage. They too observed 
more crushing for uniformly-graded soil than for well-graded 
soil. 

Bishop (1966) showed that considerably more crushing occurs 
during the shearing stage of triaxial test than during the 
consolidation stage. 

Lee and Seed (1967) showed increasing degradation of 

Sacramento river sand (having 4% material finer than 75 microns) 

with increasing confining pressure. It Increased to 14%, 22% , 40% 

2 

and 52% under confining pressure of 10, 40, 80 and 140 Kg/cm . 
In case of well-rounded Ottowa sand, they reported that It was 
much less compressible than fine sub-angular to sub- rounded 

if 

sand, and the energy required for crushing of Ottowa sand was 
shown to make up to a significant, component of the total observed 
Shear strength of granular soil under high pressure. 
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Vesic and Clough (1968) performed triaxial teats on 50 
samples of fine Chattahoochee river sand (medium, uniform mica- 
ceous sand, composed of subangular quartz grains, d^p = 0.37 mm, 
*max ~ I-** e min = 0-^1) ancl reported that at very low pressure 
there is very little crushing . However , the sand particles were 
relatively free to move with respect to each other and the dila- 
tancy effect could be quite significant. As the mean normal 
stress Increased, the crushing became more pronounced and the 

di latancy effect gradually disappeared and the crushing was more 

2 

intense in the elevated pressure range (10 - 100 Kg/cm ) until a 

"breakdown stress" was reached beyond which the sand behaved 

essentially as a linear deformable sand . 

Ramamurthy (1970) studied the influence of fragmentation of 

particles of Badarpur sand with two different gradings, in tria- 

2 

xial compression test with confining pressures upto 70 Kg/cm and 
showed that the magnitude of. crushing Increased with mean 
principal stress during isotropic consolidation and shear , but 
major crushing occured during shear and in the coarser fraction 
of sands. This crushing produced major changes in the properties 
of sand like effective principal stress ratio, effective angle 
of shearing resistance, dilatancy ratio and coefficient of earth 
pressure at rest. 

Koerner (1970) conducted tests on quartz, feldspar and 
calcite and reported more compressibility and degradation for 
uniform soils in comparison with well-graded soiljs having 
maximum particle size. 


same 



Blllan (1971) shows that some aspects of high pressure 
behaviour are strongly influenced by the tensile strength of the 
grain substance. Grain strength may be easily estimated and used 
as a scale factor for comparing the behaviour of initially dense 
granular materials sheared at high pressures. 

Oda and Minami (1976) tested residual soil derived from 
weathering of granitic rock (called Masado ) using direct shear 
tests and found that, due to effect of particle breakage, the 
failure envelope derived from direct shear tests under constant 
volume and constant normal stress is never a straight line. 

Miura and Sukeo-o-Hara (1979) studied the effect of particle 
crushing on shear characteristics of a decomposed granite soil, 
using a series of triaxial compression tests under drained and 
constant confining pressure conditions and found that the 
particle crushing phenomena, under low stress, of a soil consis- 
ting of breakable particles Is substantially the same as the 
particle crushing phenomena under very high pressure of soil 
consisting of strong particles. Further they showed that the 
particle breakage rate has a close relationship with volume 
change characteristics of the sample tested. 

Datta, Gulhati and Rao (1980) studied the crushing of 

calcareous sands during drained shear, using triaxial 

2 

compression test and confining pressure of 64 Kg/cm , and found 
that calcareous sands with skeletal particles are significantly 
more susceptible to crushing than terrigeneous quartzitlc 
material. The extent of crushing is Influenced by particle 
characteristics and void ratio. They further pointed out that 
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there i a significant reduct i on in maximum principal effective 
stress ratio due to particle crushing. 

Dfttta, Gulhati and Rao (1982) studied the engineering 
behaviour of carbonate soils of India and reported that crushing 
of carbonate particles and cementation by carbonate material are 
the two most dominating factors which influence the engineering 
behaviour of carbonate soils. 

Vaid and Chern (1985) studied the effect of pressure and 

angularities on liquefaction of Ottawa and mine tailing sands 

and showed that there is some particle breakage in case of mine 

tailing sand (this implies that both consolidation and shearing 

of angular sands results in some breakage of sharp edges of 

particles, with no gross particle crushing under the level of 

confining pressure used) whereas no particle breakage could be 

2 

observed for the Ottawa sand even at a stress of the 40 Kg/cm . 

Phool Chand (1988) performed direct shear test on angular 
sands (Ganga and Kalpi) and well rounded quartz sand (standard 
sand) and observed significant amount of particle crushing for 
the angular sands with no crushing for uniformly-graded rounded 
quartz sand. 

2.6 SCOPE OF WORK 

The foregoing review covered the following aspects of 
response of granular materials subjected to both normal and shear 

stresses : 

1) Characterization of particle shape and its Influence on 
mechanical behaviour of sands. 
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2) Stress-strain dilatancy behaviour of sands under low and 
high confining pressures. 

3) Compressibility of sands and the significant controlling 
factors. 

4) Mechanisms of grain modification and crushing under normal 
and shear stress and its influence on engineering properties 
of sands. 

The review reveals that while geologists, dealing with large 
size particles, have successfully developed methods to quantify 
particle morphology, the same has not been done in case of sands 
in a standardized way. Highly descriptive and therefore subjec- 
tive descriptions of particle shape are usually given. In this 
study, an attempt is made to develop a quantitative procedure to 
describe particle morphology. This is considered necessary since 
it is well recognized that particle morphology has significant 
effect on compressibility, particle degradation, liquefaction 
potential etc. 

Literature review also reveals that data on natural sands 

demonstrating the effect of grain morphology and mineralogy is 

rather limited, and no systematic study involving a variety of 

natural sands is available. Most of the results reported relate 

usually to quartz sands with subrounded to subangular particles. 

In this study three natural sands with different mineralogy and 

mode of formation along with an artificial standard rounded 

quartz sand, have been used to investigate compressibility beha- 

2 

viour upto a normal stress of 176.32 Kg/cm . In addition. 
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different mixtures /eradat ions of sands were prepared by inter- 
mixing the four main sands and morphological and mineralogical 
characteristics and their compressibility were investigated. The 
relationship between relative density and height of free fall 
(using raining technique) were investigated for all the sands. 

The compression tests were also utilized to examine the 
nature and degree of grain modification/crushing. These results 
have been critically examined in terms of stress level, grain size 
and the nature of minerals present. The grain morphology after 
degradation was also quantified and compared with the unstressed 
natural sand particles. Nature of the laws governing variation of 
constrained modulus (PI) with effective confining stress and 
relative density has been Investigated for sands of varying 
compressibilities. Based on the Influence of grain angularity on 
the angle of shearing resistance and the effect of sand compre- 
ssibility on the mode of failure under the penetration cone tip, 
the laws governing the variation of ultimate cone tip resistance 
(q c ) with effective confining stress and relative density have 
been examined. 



CHAPTER III 


MATERIALS AND METHODS 

3-1 General 

Strength and compressibility are influenced significantly 
by the size and shape of the particles and the mineral 
composition of the soil constituents- While mechanical analysis 
and optical techniques are adopted for characterizing the size 
and shape, microscopic and X-ray difraction investigations have 
been resorted to for the mineral og i cal analysis. In the present 
Chapter, the materials chosen for the study are described in 
terms of their characteristics and the test procedures adopted 
together with the details of the equipment involved are 

presented. 

3.2 Materials Used 

In the present study four sands have been chosen of which 
three are natural sands while the fourth one is of an artificial 
type. In addition to the bulk samples of these four sands, six 
samples have also been prepared by size fractionation, mineral 
separation and also by selective mixing of the same. 

Few samples of sands as obtained from sources abroad 
(Rensselaer Polytechnic Institute, New York) have also been 
studied in terms of their size and morphological characteristics. 

The data obtained were used for evaluation in the present case. 

! 

3.2.1 Ganga Sand 

The sand from the flood plains of Ganga near Kanpur, 
chosen for the present investigation has 


uniform grain size 
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distribution. Examination with binocular microscope and image 
analyzer (Plate 3.1) revealed that this sand was predominantly 
angular with minor amounts of flaky as well as rounded grains. 
The mineral composition is dominantly of quartz and feldspar with 
subordinate amounts of mica and calcite (carbonate). Fragments 
of gneiss and quartzite also occur in minor quantity. The quartz 
grains are quite transparent and often noticed to be fractured, 
and some times coated with iron oxide. The feldspars are not as 
transparent as quartz and are partly altered. Presence of mica 
reflects the contribution of phyllltic material from the foot- 
hills of Himalayas as one of the possible sources. The carbonate 
constituent was possibly from associated carbonate formation in 
the Himalayan sequence. 

Ganga sand used here is similar in many respects to the 
Chattahoochee river sand (Vesic, 1968) which is also uniformly 
graded and consists of subangular quartz grains with subordinate 
mica . 

3.2.2 Calcareous Sand 

Calcareous sand used in the present study was obtained from 
the area in vicinity of Okhla light house (Gujarat). The sand is 
well-graded with subrounded to subangular grains mostly composed 
of sea shells (Plate 3.2). The X-ray diffraction pattern for the 
: same revealed dominant preserice of carbonate (aragonite and 

calcite) with subordinate amount of quartz. The source for this 

. 

n , 

carbonate material is from the oolites and also from the skeletal 
remains of marine organisms. The thin-walled skeletal calcareous 

P 

grains are characterized by the typical presence of intra- 
particle voids that increase the tendency of the material to 
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a. Typical morphology of shells. (x 36) 



b. Subrounded and Pitted Carbonate ^ 
grains. 


Plate 3,2 Micrograph of Calcareous Sand.(SEM) 
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crush under stress. 

3.2.3 Kalpi Sand 

The sand from Kalpi used in this study is from Yamuna river 
system. In the region west of Kalpi, ' two tributaries 
(Kanchamalanga and Non) of Yamuna have their confluence at 
Gorakalon. The sand from the flood plain of river is very 
characteristic in its features and is distinctly different from 
the Ganga sand. The particles are well-graded, relatively coarse 
and very angular in outline (Plate 3.3). Quartz and feldspar are 
present in equal proportion ( 40% each). The carbonate fraction, 
composed of calcite, is present in considerable quantity (upto 
about 187.) rendering the sand calcareous as compared to the Ganga 
sand in which the carbonate fraction is in minor amount only. 
The carbonate material is present largely as matrix agglomerated 
with quartz and feldspar grains. When treated with mild acid, 
this matrix dissolves leaving the silicate particles that 
constituted the nucleus for the growth of carbonate. In 
addition, carbonate is also present in minor amount in the farm 
of limestone grains. 

3.2.4 Standard Sand 

This sand, obtained from Ottawa (USA), is uniformly graded 
and is entirely of quartz. The grains are coarse, nearly 
spherical and highly rounded in outline (Plate 3.4). While the 
quartz grains in Ganga sand are very transparent, the grains in 
this sand are opaque with frosted surface possibly acquired 
during the process of its preparation. 






and quartz 
(x 150) 


Plat* 3.3 


Micrograph of Kalpi Sand 


(SEM) 



Plate 3.4 Micrograph o£ Standard Sand 

indicating rounded grains of 
quartz. ( x 46) 


CSEM) 
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3.3 Experimental 

Estimation of grain size distribution, void ratios and 
specific gravity have been carried out adopting standard 
procedures using the conventional test facilities in soil 
engineering practice. Using sieves of 4.75, 2.00, 1.00, 0.600, 

0.425, 0.212, 0.150, 0.075 and, 0.045 mm, sand specimens were 

analyzed for their grain size before and after loading. For the 
estimation of maximum and minimum void ratios as well as specific 
gravity, tests were conducted on. four specimens of each of the 
sand samples to obtain consistency of results. 

3.3.1 Preparation of Samples at Specific Relative Density 
Values 

For compressibility and crushing, test specimens with 
relative density values of 20, 25, 40 and 657. were used for 

i 

different sands. For this purpose calibration charts have been 
prepared for relative density as a function of the height of fall 
using the raining technique. The sand was poured from a 
container (funnel) at a fixed height and corresponding relative 
density was calculated. The relative density of sand is a unique 
function of height of free fall of sand as long as a constant 
rate of flow is maintained. Details of experimental programme are 
given in Table 3.1. 

3.3.2 Mineralogical Characterization 
Mineralogical composition of the sands used in the present 

study was determined using X-ray diffraction technique. This 
analysis » has been carried out on the General Electric Unit, 
operated at 30 KVA with 20 ma current, fitted with XRD-5 
diffractometer using nickel-filtered Copper radiation. Scanning 
was done between 5 and 50 degrees (2 0) at 2° (2 0) per cm. 




Table 3-1 Experimental Programme 


SI* Type of Testa No. of Sand Particles Total 
No. S amp I e s per s a mp 1 e 


1 . 


Morphology 

2 . 


Mineralogy 


a) 

without acid 
treatment 


b) 

with acid treatment 

3 . 


Gradat ion 


a) 

before loading 


b) 

after loading 

4 . 


Maximum void ratio 

5 . 


Minimum void ratio 

6 . 


Specific gravity 

7 . 


Density (Raining 
t eclmique) 

0 . 


Coiiiprcaa ion 


a ) 

relative density 

2 5 Z 


l>) 

relative density 


6 5 \ 


23 

100 

2 3 0 0 

10 

-- 

18 

0 



10 
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6 Z 

10 

-- 

10 

10 

-- 

10 

10 

-- 

I 0 

9 

-- 

9 

25 

~~ 

42 

17 

-- 
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Particles 


Tests 
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tt 


tf 
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if 


t* 


tf 



Identification was done by converting the "2 0" values of the 
peaks into the coresponding "d" values using Bragg’s equation and 
comparing the "d" values and intensities of the peaks with the 
data reported in the standard file. 

3.3.3 Compressibility 

The test specimen prepared with different relative density 
values in oedometer were subjected to varying stress levels upto 
176.32 Kg/cm^. For this purpose Material Testing System (MTS 
810) was used. This test facility has the following special 
features: 

i) Servo-hydraulic controls for smooth running 

ii) Electronic programming and read out. 

iii) Modular construction with the advantages of stroke, 
strain and lead controls. 

iv) Hydraulically operated grips for easy insertion and 
removal of test piece. 

v) X-ray recorders to suit various types of test output, 

vi) Large capacity of 10 tonnes with load cell having foil 
type strain guage bridge bonded to high 
strength metal-column. The bridge provides an out- 
put voltage proportional to applied load. 

The system operates on a close loop principle. A command 
signal in the form of analog programme voltage, representing the 
desired load, stroke or strain to be applied to the specimen is 
compared with a feedback signal that represents the actual load. 
Stroke or strain is measured by transducer. Any deviation 
between the command and feedback causes a correction of the 
control signal to be applied to a servo-valve which in response 
to the control signal causes the actuator to move in a direction 



2 up 

w Q *' 

3 2 S' 

2 2 si 


iSo 2 

u -* o 

^ U 


p* -1 

! i 1 


8 k T 

§£ 

°g 1 


xa 

uo 

< H 
OQ (J 

CONTROL 

cc 

a 

5 

qc 

1- 

z: 

b ui 


o 

o 

w -j 
U l uj 


o 

a. to 


> 



cr 

in 

V) 



— cr 
z o 

O 

;:< 


o cr 

Is 

£ s 

— . ^ 


, I o I 

cr ^ I 5 I 

Ik t » 'I* — . «M» m* 



I 1 

1 | o:l 

* JL I** jjjjj f 


— 

^ I ujl 
ol_ q J 

QC 

O 


t=r 

O £ 

j 1 

-J CL 

— \ -'i 

r\ 

< 1 /) 
fY _ 


£ a 

Q iu 

>" 

i s 


(L 


Fig. 3.1 Block diagram of Material Testing System -Model 810 













required to reduce the deviation to zero. Block diagram of the 
Material Testing System is given in Fig. 3.1. 

3.3.4 Particle Shape Analysis 

The particles in each of the sand specimens before and 
after the compressibility tests were analyzed for shape in terms 
of their sphericity, shape factor and angularity. This involved 
estimation of long, intermediate and short diameters as well as 
tangent counts and projected surface area for different grains in 
each of the specimens. Using these parameters, quantitative 
evaluation of shape has been carried out. 

An image Analyzer System of Bousch and Lamb make (Omnicon 
Alpha 500) has been used for this purpose. A brief description 
of the system is given in the succeeding paragraphs. The system 
is versatile and is capable of providing full field measurements 
or measurements of individual single feature in addition to the 
tangents and intercept counts. The system consists of two major 

assemblies the microscope with the scanner and the rack of 

modules - basic , display, control and measurement modules* 

3. 3. 4-1 Microscope and Scanner , 

A microscope is essentially an optical image forming device 
used in conjunction with the analyzer. The scanner resolves the 
optical image presented to it by scanning along parallel lines, 
producing an analog signal (as a function of position) which 
corresponds to the intensity profile of the image along each 
line. The overall function of scanner is therefore, essentially 
similar to that of television camera, but with a different design 
for obtaining excellent linearity, stability, shading and noise 
performance, as well as an ability to accommodate rapid changes 
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la the level of illumination depending on t;he specimens. 

3. 3. 4. 2 Basic Module 

The basic module supplies power to the Measurement Module 
and generates count tags, bright outline demarcating boundaries 
of the various features and the alpha-numeric characters that 
indicate the counts. After processing the video signal from the 
scanner, this module feeds the same into display. 

3.3. 4.3 Diplay Unit 

The Diplay Unit is an instrumentation type television 
monitor with a diagonal screen size of 229 mm. 

3. 3. 4. 4 Entire Field Count Unit 

The entire field count unit counts the selected features or 
their components in the image. It receives, its input in the form 
of binary video signal from the Basic Module and provides Binary 
Coded-Decimal (BCD) output back to the- Basic Module, (Plate 3.5). 

3. 3. 4. 5 Measurement Unit 

The Measurement Unit measures the desired dimension of 
detected features. Measurement of full field or any single 
feature Is possible with this unit. In the present study the 
latter has been adopted. Block diagram of the image analyzer is 
presented in Fig. 3.2. 

3.3.4. 6 Particle Measurement 

The measurement of diameters , projected surface area and 
tangent counts have been made using the image analyzer system. 
These data were used for the quantitative analysis of the 
particle shape. For determination of diameter and surface area 
(area excluding hole), single feature measurement mode of the 
image analyzer was used. These measurements are possible on 
features selected by the operator using either the light pen or 
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Plate 


3 . 5 


Image Analyzer Set Up 


Gray Level 

Detective Control 



Fig. 3.2 Block diagram of the Omnican Alpha Image Analyser with Microscope 
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madn to rest on their respective maximum pro j ec t ed 

a r e a . 

iii) The sample was placed on the movable stage under 
1. 1 1 e m i c r on cope . 

iv) An average of 25 to 30 particles from each sieve 
fraction were randomly chosen for the analysis. 

v) The shortest diameter (do) for each grain was deter- 
mined by focussing alternately on the glass plate and the 
top of grain (using the Vernier reading of the 
microscope) and measuring the differential deviation. 

vi) Using the movable frame facility, ” tangent: rectangular" 
for each individual grain was constructed. The long side 
of the tangent rectangular measured in vertical direction 
was taken as long diameter (djJ). 

vii) For determination of intermediate diameter (dj) the 
movable stage along with the glass slide was 

rotated through 90 degrees ' and the short side of 

already constructed tangent rectangular was measured in 
the vertical direction. 

The Count Mode of the image analyzer was used. to determine 
the tangent count. Four types of counts are possible - total 
feature (particle) count, lower positive tangent ,’ 1 over negative 
tangent and the intercept count. In this study, the lower posi- 
tive tangent (the sum of all tangent points formed at the lower 
convex boundaries of the particle with horizontal scan line) and 
lower negative tangent (the sum of all tangent points formed at 
the lower concave boundaries of the ' particle with the horizontal 
scan line) were counted. The negative tangent count acted as a 
check because its value is equal to positive tangent count minus 
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two. For tho dot grminat ion of total numb or of protrusions on tho 
upper portion of the grain can be counted by using either the 
dark mode of the image analyzer or by rotating the movable stage 
of microscope along with the glass slide by 90 degrees. Lower 
positive and negative tangents and features engaged in each other 
are shown in Fig. 3.4. 

3. 3. 4. 7 Quantification of Particle Shape 

Shape is a complex property of a grain and it is difficult 
to define shape properly except for grains that approximate 
regular geometrical shapes. An already indicated, roundness 
(angularity) and sphericity (shape factor) aid in quantifying the 
particle shape. In the present work, • detailed investigations 
were undertaken for the quantification' of shape. 

3. 3. 4. 7.1 Sphericity and Shape Factor 

Krumbein (1941) expres-sed sphericity iti terms of the three 
diameters ( d j , d j , d ^ ) of a triax.ial ell ipso id . A sphericity 
value of one implies a perfect sphere, and , the sphericity value 
decreases with the irregularity of .the ' par I. i c 3 e shape. based on 
measurements of d^> d j and dg particle sphericity may be computed 


as : , 

Volume of particle 

Y = 

Volume of circumscribed sphere 



(3.1) 


Shape factor Is another index used as a measure of 

shape. Lees (1964) proposed a shape factor 

c 

termg of flatness ratio P ( = d s /dj) and elongation 

b 


q(b/a = dj /d L ) . 


par tide 
(F) in 
r a t i o 



68 



Two features engaged Lower Positive Tangent 

in each other are 
counted as .one 



Fig. 3,4 Lower positive and negative tangents and two 
features engaged in each other. 
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In thin study, shape factor was calculated using the 
expression (Blatt et; al . 1972 ): 



^ 3 L d I (3-2) 

whore dg is the short diameter, 

dj is the intermediate diameter, and 
dj^ is the long diameter 

Procedure for the measurement of d^, dj and d^ has already been 
presented. 

3. 3. 4. 7. 2 Roundness (Angularity) 

Roundness refers to the sharpness of the corner's and edges 

* of the grain. It therefore refers to aspects of the grain 

t 

surface that are on larger scale than those covered under surface 

i ! 1 

texture but still smaller than the overall dimension of the 
grain. Uadoll (1932) defined roundness (angularity) as the ratio 
of fcho average radius of curvature of the corners to the radius 

fK 

of largest inscribed circle. Powers (1953) proposed a procedure 

for estimating particle roundness by visual comparison of grains 

, 

with standard images of grains of known roundness. The roundness 

- 

index assigned as per Powers’ classification is given in Table 

2 . 2 . 

■ 

The procedure (Youd, 1973) used for the determination of 
particle roundness index (R) is as follows: 

ifelf i 

1. Separate the sand ( < 2 mm ) into different sieve 

fractions . 

2. Examine under microscope a representative number of 
particles (at least 50) from each sieve fraction. 




Plat e 3.6 


Projected Image 
Indicating Edge 


of a Particle 
Protrusions 


3. Vi ou.illy compare shape of each particle with the shape 
classes and assign the value of roundness index (R) 
according to the classification scheme listed in Table 
2 . 2 . 

4. 'Die average value for each sieve fraction (Rj) was 
then calculated as 


Rj = C I Ri)/» 


C 3 - 3 ) 


where Rj is the R value assigned to par tide 
i , and 

u is the total 1 number of particles 
examined. 

S. After computing Rj for each sieve fraction, an 
average value of Rj for the sand under study was 
computed from the relationship 

R = ( £ R j P j)/10° ' ■ (3.4) 

where Pj is the percent (by weight) of the 

particles .in the applicable sieve 


fraction. 

The estimation of roundness index (R) by visual comparison 
of shape would appear to be some what subjective and therefore an 

; 

indedpendent approach was adopted to quantify roundness 
(angularity). Using the image analyzer, a count of total number 

of corners or protrusions was obtained on the boundary of the 

T 

: 

particle linage^ pro j ec t ed on the display screen (Pliite 3.6). The 
linage analyzer provides a total tangent count for both convex and 
concave features on the grain boundary. . Tangent count was taken 
a measure of grain angularity. The procedure involved in the 
stimation is as follows: 



f> i i * v •• .'vnalyais of each sand sample was carried out. 
Sufficient number of particles (25-30), from each siz<5 
fi action were placed on a glass slide which af Let enough 
tapping (to ensure that the particle rest..*; on the plane 
of its maximum projected area) was put on the movable 
stage under the microscope. 

Each grain on the display screen is boxed in the ’’tangent 
rectangular" (by movable frame facility) and using the 
count mode of the image analyzer, the number of tangents 
on the lower and upper grain boundary is 1 ecoi d< d. 

The stage was then rotated through 90 degrees and the 
tangent counts were recorded for the other two side.e. 

Sum of tangent counts obtain in steps 3 and 4, provide a 
tangent count value for each particle. 

The average value of tangents for each sieve fraction was 
then calculated as 

Tj=( I Tj)/n (3.5) 

where T i is the total number of tangent 
count for particle i , and 
n is the total number of pax Liclos 


examined . 

After computing Tj for each sieve fraction, an average 
value of tangent count for the sand under study was 
determined from the relationship 


T = ( It j P j/n)/100 (3.6) 

where Pj is the percent (by weight) of the 
particle in the corresponding sieve 


frac t ion . 
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The average value of tangent count thus; obtained war; directly 
unnd for comparing different sanda . However, tor comparing 
particle angularity for different aieve fraction of a given a and , 
the average tangent count for each fraction wan normalised with 
the average circumference of the grain.'! in that fraction. 

It in important to note that both roumlnoss index and 

• * \* 

tangent count arc? obtained by observing the grain images in two 
diiiumu iorni . A procedure similar to the one adopted hero for 
tangent count has also been used for averaging the values* of 
particle sphericity and shape factor. 
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CHAPTER IV 

EXPERIMENTAL DATA 

4.1 General 

Experimental investigations have been carried oat to 
characterize all the sand samples involved in the present study 
in terms of their mineralogi cal , physical and morphological 
characteristics. The samples subjected to the compressibility 
tests have been analyzed for their particle size and shape to 
enable an understanding of their role in the compressibility and 
also the extent of crushing the particles in the system undergo. 
In this Chapter, the experimental results are presented. 
Interpretation of these results along with the discussion of the 
inter relationshi p between various parameters involved during 
compressibility and crushing of the sands is dealt in the 
succeeding chapter. 

4.E Mineralogi cal Composition 

X-ray diffraction analysis has been carried out for the 
different sands used in the present study. Various mineral 
constituents have been identified on the basis of their 
characteristic, peaks. Semi-quantitative estimation of the 
mineral constituents has been made which has been confirmed to be 
fairly within the range on the basis of microscopic examination 
as well with the acid dissolution method in the case of carbonate 
f ract ion . 

Typical diffraction pattern for the Ganga sand is presented 
in Fig. 4.1. Quartz is the dominant mineral constituent (60 to 
657. ) associated with feldspar (20-257). In addition, presence of 
mica (8-107) is confirmed to be present from the occurrence of 
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-ray diffraction pattern of Ganga sand . 
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the peaks around 10 A and 5 A as well as others marked in the 
pattern. Minor amount of carbonate in the form of calcite is 
also present. The peak around 12.4° (20) could either be of 
kaolinite or chlorite or both. From the presence of peaks around 
36.1° and 38.2° (20) , it has been confirmed that kaolinite is 
present in a very minor amount. On the other hand, the peak at 
6.1° (20) corresponding to 14.2 A confirms minor amount of 

chlorite also. Thus chlorite and kaolinite are both present upto 
about 37. within the sand. 

The X-ray diffraction pattern for Kalpi sand (Fig. 4.2) 
reveals distinct mineralogi cal composition as compared to Ganga 
sand. This sand is relatively rich in calcite with the carbonate 
fraction estimated to be around 187.. Quartz and feldspar are in 

equal proportion (around 40% each). In this sand, the feldspar 

I 

however is of two varieties (alkali and potash types) indicating 

|h: 

a distinctly different source for .the sand as compared to that of 
Ganga. The occurrence of mica is in minor quantity (upto 27) in 

this case. 

The diffraction pattern for calcareous sand (Fig. 4.3) has 
confirmed the presence of carbonate upto 957. of its composition 
with quartz as the associated mineral (upto 57). X-ray analysis 
has very clearly established the presence of two carbonate 
minerals namely aragonite and calcite. As already indicated, the 
association of aragonite reflects contribution of carbonate from 
marine organisms. Electron microscopic examination of carbonate 
particles in this sand has revealed a distinct mo r pho log i cal 
pattern for the grains typical of a sea shell surface. 

I The standard sand used is a commercially prepared one. It 
is a pure quartz sand without any other mineral constituents as 


- Quartz 

- Feldspar 



Fig. 4.2 X-ray diffraction pattern of Kalpi sand* 
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confirmed from its X-ray diffraction pattern (Fig. 4.4). 

X— ray analyses have also been carried out for estimating 
the mine ral og i cal composition of other sand samples namely the 
Ganga sand after removal of mica, Ganga with enhanced mica, the 
sand obtained by mixing Ganga, calcareous and Standard (referred 
to as coarse Ganga), artificial calcareous sand, prepared by 
mixing Ganga, Kalpi, and Standard sands and the two fractions of 
Kalpi sands with different grain size distribution (designated as 
Kalpi 1 and Kalpi 2). The respective X-ray diffraction patterns 
for these sands are presented in Figs 4.5 - 4.10. Estimated 
mineral composition for each of these sands is summarized in 
Table 4.1 and 4.2. The estimation has been fairly within the 
range as has been checked using the acid .dissolution technique as 
in the case of carbonate fractions. : 

4.3 Physical Characteristics 

The grain size distributions for four main sands 
(Ganga, Kal pi , Calcareous and Standard) are shown in Fig. 4. 11, 
whereas Fig. 4.12 shows the gradation curves for Ganga with 
enhanced mica, Ganga with mica removed. Coarse Ganga, artificial 
calcareous, Kalpi 1 and Kalpi 2 sands. The grain size 
distribution for Kalpi sand (Fig. 4. 11) is distinctly different 
Compared to the other sands. It may be indicated that Kalpi sand 
las its source from the Yamuna with a relatively shorter distance 
)f transport. 

The relation between the relative density versus height of 
all using raining technique for Ganga, Kalpi, Calcareous and 
standard Sands is shown in Fig. 4.13. For a comparison data for 
he Chattahoochee river sand (Vesic 1963) is also represented on 
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Fig. 4.4 X-ray diffraction pattern of Standard sand. 
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Fig. 4.7 X-ray diffraction pattern of Coarse Ganga sand 
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Fig. 4.10 X-ray diffraction pattern of Kalpi 2 sand. 
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Table 4.1 Mineralogi cal Composition of .the Test Samples 


SI. Sand Quartz Feldspar Mica Carbonate Chlorite 

No. Sample ■/. '/. 7. 7. and Kaoli- 

n i te , V. 


1 

Ganga 60- 

-65 

20-25 

8-10 

2-3 

2~3 

2 

Ganga with 30- 

enhanced 

-40 

6 

50-60 

4 

3 


mica 






3 

Ganga with 
mica removed 

65 

30 

1-2 

3-4 

- 

4 

Coarse Ganga 26- 

-27 

2-3 

71 


- 

5 

Kal pi 

40 

40 

1-2 

18 

- 

6 

Kalpi 1 

45 

30 

- 

25 

- 

7 

Kalpi 2 

50 

18-20 

1-2 

30 

- 

8 

Cal careous 

5 

— 

— 

95 

— 





( aragon i te 

> and 

cal ci te ) 

9 

Artificial 

calcareous 

50 

20 

5 

25 

- 

•10 

Standard 

100 

- 

- 

- 

- 

11 

Kalpi (after 
Acid Treatment) 

50 

48 

- 

2 

- 

12 

Kalpi 1 (After 
Acid Treatment) 

55 

40 

- 

5 

- 

13 

Kalpi 2 (After 
acid treatment) 

55 

40 

- 

5 

- 
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Table 4.2: Estimates of Carbonate and Noncarbonate Mineral 
Contents (on the basis of acid dissolution method) 


SI .No. 

Sand 

Carbonate , 
7. 

Other minerals 

7. 

1 

Gan ga 

5 

95 

2 

Coarse Ganga 

71 

29 

3 

Kal pi 

20 

80 

4 

Kalpi 1 

26 

74 

5 

Kalpi 2 

29 

71 

6 

Cal careous 

91 

9 

7 

Artif i cial 
Calcareous 

20 

80 


8 


Standard 


0 


■100 



ercent Fingf Percent Finer 
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Fig. A. 11 Particle size distribution. 
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this Figure. The relationships for the rest of the sands are 
represented on Fig. 4.14a. As indicated earlier, the density of 

the sand is a unique function of height of free fall of sand as 

long as the rate of flow remains constant and beyond about 20 cm 
of height of fall, the rate of increase of relative density 
drastically varies. Effect of rate of flow has also been 

reported by Sahu (1988) and his results are presented in Fig. 
4.14b. These relationships show that the particle gradation has a 
strong influence on the degree of compaction of sands. Various 
physical properties of sands is summarized in Table 4.3. 

4.4 Morphological Characteristics 

Image Analyzer System was< used to find out the three 

diameters the long (d L ), the intermediate (dj) and the short (dg) 
diameters respectively. Using these three values, ratios of the 
intermediate to the long diameter (elongation ratio) (djVd|_) and 
the short to the intermediate diameter (Flatness ratio) (dg/dj) 
are determined for individual particles pertaining to each of the 
different size fractions. The particle shape classification for 
the sands, before and after loading, was done based on the Zingg 
shape concepts (Figs. 4.15 to 4.35). The relation between the 
shape factor and the flatness ratio (dg/dj) for sands before 
and after loading is indicated in Figs. 4.35 to 4.56. In addition 
to the study of particle shape, data were also obtained for each 
sieve fraction in term of number of tangents which reflect the 
surface characteristics of the edges of the particles. 

4.5 Compressibility 

To evaluate the compression characteristics, tests on dry 
sands at loose and dense relative densities (desired relative 
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Table 4.3: 

Physical 

Properties of 

the Sands 




SI . 

Sand 

G 

d 6Q 

d 50 

d 10 

c u 

e max 

e min 

1 

Ganga 

2.691 

0.218 

0.18 

0.085 

2.565 

1 . 250 

0.577 

2 

Kal pi 

2.676 

1.300 

1.00 

0.270 

4.814 

0.910 

0.480 

3 

Cal careous 

2.700 

0.480 

0.43 

0.270 

1 .778 

1 . 050 

0.570 

4 

k 

3 

Standard 

2.650 

0.500 

0.475 

0.360 

1.389 

0.810 

0.490 

Ganga with 
mica 

2.730 

0.199 

0.200 

0.080 

2.488 

1.229 

0.557 

6 

Ganga without 2.670 
mica 

0.218 

0.180 

0.100 

2. 180 

1.110 

0.718 

7 

Coarse Ganga 

2.675 

0.460 

0.400 

0.180 

2.565 

0.773 

0.485 

8 

Artif i cial 
Cal careous 

2.669 

0.500 

4.800 

0.280 

1 . 780 

0.904 

0.516 

9 

Kal pi 1 

2.695 

1.050 

1.020 

0.950 

1.073 

0.932 

0.540 

10 

Kal pi 2 

2.710 

1.950 

1 . 700 

1 . 180 

1.613 

0.989 

0.576 
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Fig. 4.16 Shape classification using Zingg 
sand subjected to 22.04 kg /cm 2 . 
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Fig. A. 17 Shape classification using Zingg 
sand subjected to 176.32 kg /cm 2 . 
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Fig. 4.18 Shape classification using Zingg 
Ganga sand. 
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Fig. 4.20 Shape classification using Zingg diagram for Kalpi 
sand subjected to 22.04 kg/cm 2 . 
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Fig. 4.24 Shape classification using Zingg diagram for calcareous 
sand . 
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Fig. A. 25 Shape classification using Zingg diagram for calcareous 
sand subjected to 22.04 kg /cm 2 . 
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Fig. A. 26 Shape classification using Zingg diagram for calareous 
sand subjected to 176.32 kg / cm 2 . 
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Fig. 4.27 Shape classification using Zingg diagram for artificial 
calcareous. 
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Fig. 4.29 Shape classification using Zingg diagram for 
Ottawa sand , 





110 



Fig. A. 30 Shape classification using Zingg diagram for 
San Fernando =# 5 sand . 
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Fig. 4.31 Shape classification using Zingg diagram for 
San Fernando# 6 sand. , 
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Fig. A. 32 Shape classification using Zingg diagram for 
San Fernando 7 sand. 
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Fig. 4.33 Shape classification using Zingg diagram for 
Leighton Buzzard sand. 
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Fig. 4.34 Shape classification using Zingg diagram for 
Lag uni lias (94A-M19) sand. 
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Fig. 4.35 Shape classification using Zingg diagram tor 
Lagunillas (94A+9AB) sand. 
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Fig. 4.60 Particle size gradation before and after loading for 

Ganga sand . 



Fig. 4.61 Particle size gradation before and after loading for 

Calcareous sand 
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Fig. 4.63 Particle size gradation before and after loading for 


Standard sand. 
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Fig. 4.64 Particle size gradation before and after loading for 

Kalpi 1 sand. 
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Fig. 4.65 Particle size gracjation before and after loading for 

Kalpi 2 sand. 
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Fig. 4.66 Particle size gradation before and after loading for 
Artificial Calcareous sand. 
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Fig. 4.67 Particle size gradation before and after loading for 
Coarse Ganga sand. 
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Fig. 4.68 Particle size gradation before and after loading 
for Ganga mica removed and Ganga enhanced mica 





CHAPTER V 


INTERPRETATION AND DISCUSSION OF TEST DATA 

Ihe engineer ing behaviour of a sand is controlled by 
its grain morphology and mineral composition. While sporadic 
attempts have been made to characterize sand in terms of its 
size and shape attributes, there is need for a quantitative 
classification of the same in terms of size, angularity 
(roundneoB ) , sphericity and shape factor . Such a 

characterization is essential for continuously tracking the 
grain modification at different stress levels. The 

mineralogy of the sand system plays an important role in the 
grain modification. As detailed in the previous chapter, 

series of experiments have been conducted to accomplish not 
only the quantification of the grain morphology but also to 
gain an understanding of various parameters operative during 
the sand compressibility and interrelationship between the 
same . In the present chapter, the test data obtained are 
interpreted in the light of various operating parameters and 
the details are presented on morphological classification, 
role of grain modification and crushing on compressibility 
and constrained modulus. The role of mineralogical 
composition in the behaviour of the sands under stress has 
also been brought out. 

5.1 GRAIN MORPHOLOGY 

In this section results of studies related to particle 
sphericity, shape factor and angularity are discussed on the 
basis of which the particle morphology of a variety of 
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natural and artificial sandy was characterized. 
Interrelationship between various parameters used for 
description o£ ssui face characteristics is’ established. 
Typical distributions of surface characteristics for 
different systems from rounded to highly angular nature of 
particles are presented. 

5.1.1 Sphericity and Shape Factor 

Particle shape can be designated as oblate, equant 
(spherical ) , prolate (toller, or rod like) and b laded (flaky) 
depending the location of the data plots in quadrants I, II, 
III, IV respectively on the zingg diagram (Fig. 2.1). The 
data pertaining to sphericity and shape factor for each of 
the size fractions of the sand used in the present study have 
been characterized for the particle shape using the Zingg 
diagram (Fig. 4.11 - 4.56). In case of four sands studied in 
detail , about 831 of particles for Standard sand lie within 
quadrant II, Indicating spherical nature of the grains and 


about 

15% 

of 

the 

part icl es 

lie 

along 

the 

bo uridary 

of 

quadrants 

I I 

and 

III. 

In sharp contrast 

to 

the 

Standard 

sand , 

Calcareous 

and 

Gang a 

Sands 

have 

only 

2 2 

to 

25% 


particles in quadrant II and roughly 50* of the particles for 
these sands lie in quadrants I and IV, depicting the plate 
like and flaky nature of the grains. In case of Kalpi sand 
about 37% of the grains lie within quadrant II with an equal 
percentage of grains located in quadrant III indicating 
elongated nature of the grains. 
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nont; of the other natural sands examined (Table 5.1) 
show toughly 1U -5U-*. gj a ins in quadrant 1 1 and 30-10% in 
quadrant 111 indicating abundance of bulky spherical and 
elongated grains. Only 10-20% of the grains are oblate and 
5% are flaky. In addition to Zings diagram with superposed 
sphericity value-,, the shape factor (SP) for all the sands 
was also calculated (Pigs. 1.36 - 1.56). Shape factor data 
gives much better separation of different grain shapes 
belonging to various quadrants in Zingg diagram. For 
example, in case of Standard and San Feranando sands, Zingg 
diagrams (Figs, 4.20 and 1.30} would appear to suggest as if 
they have similar distribution of grain shapes. However, 
plots on the shape factor diagrams (Figs. 1.19 and 1.51) 
clearly bring out She abundance of plate and rod like 
part id oh in San Fernando sand as compared to Standard sand. 
Hence for a complete characterization of the shape of sand 
grains, shape factor diagrams are essential in addition to 
the Zingg diagram. The percentage of particles in each 
quadrant for different sands is presented in Table 5.1. 

5.1.2 Particle Angularity 

In addition to the study of particle sphericity and 
shape factor , to estimate particle angularity data were also 
obtained for each sieve fraction in terms of number of 
tangents which reflect the surface characteristics of the 
edges of the particles. As explained earlier, the image 
analyzer was used to examine the two dimensional projection 
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! 

Table 5.1 : Percentage of Particles in Each Quadrant 


Percentage in different quadrants 
Sands 



i 

I I 

III 

IV 

Remarks 

Standard 

2.47 

82 .72 

14.81 

0.00 

Plainly spherical 
grains 

Kalpi 

17.31 

37.50 

35.58 

9.61 

Rod like and 
bulky grains 

Calcareous 

29.79 

22.34 

25.53 

22.2 4 

Oblate and 

Ganfia 

28.05 

2 5.61 

26.83 

19.51 

Flakey grains 

Kalpi 1 

15. 96 

42.55 

32.98 

8.51 

Bulky and rod 

Kalpi 2 

2 5.3 1 

40.26 

18.07 

15 .66 

like grain 

Art i f 1 c i a 1 

18.49 

30.16 

36.97 

5.88 

Calcareous 

Coarse 

20 . 88 

32 . 29 

30.21 

16.67 

Oblate and Flakey 

Ganga 

Laguni lias 
94A - 94B 

14.10 

51 ,28 

30.77 

3.85 

Rod like and 
bulky grains 

Laguni 1 1 an 

2 0.00 

40.00 

32.50 

7.50 

II ft 

HA ~ Ml 9 

|, 

Leighton 

Buzzard 

11.39 

40.50 

46.84 

1.27 

II II 

SanFernando 
! No. -5 

13.04 

46.39 

39.13 

1.44 

Bulky and rod 
like grains 

SanFernando 

18.42 

38.16 

38.16 

5.26 

If II 

No. -6 

SanFernando 

15.79 

50.00 

28.95 

5.2 6 

II 17 

No. -7 

| Ottawa 

16.25 

41 . 25 

36.25 

6.25 

Spherical bulky 
grains 
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of the particles. Along with the tangent count, particle 
roundness index (R) (Powers, 1953; Youd, 1973), was also 
evaluated for each fraction and weighted average of both 
these measurements for a sand sample was calculated. On the 
basis of projected surface area data for particles obtained 
from the image analyzer, the circumference of an equivalent 
spherical grain was calculated, and the wave length ( A ) of 
surface protrusions (depicted by the number of tangents) was 
evaluated. The amplitude (h) of the protrusion was estimated 
on the basin of the measured dimension d L and the diameter, 
dp of the inscribed equivalent circle. The values of A and 
h were taken as 

Circumference of the inscribed circle 

X = (5.1) 

(No. of Tangents - 1) 

h = (d L - d p )/2 (5.2) 

The ratio h/A was taken as the measure of angularity. 

Since the average number of tangents representing a 
given sand fraction would vary with the particle size, it is 
necessary to normalize the tangent counts with the possible 
average circumference of particles for a meaningful 
comparison between the different sand fractions. However for 
comparison between different sands, average value of the 
number of tangents for each sand (results averaged over all 
the sizes and their proportions) is a good indicator of 
particle angularity. Thus the particle angularity (PA) can 
also be expressed in terms of average number of tangents (T) 
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in <*- k \ • v " 11 !si»*vr» fraction and the average mean circumference 
of an <*<iu i va 1 i ri l c i i c 1 e (having the same surface area as that 

of the particle) an 



vh er e 


(5.3) 


dp i « the equi val ent average mean diameter of the 
• circle having the same surface area as that 

of the particle. 

K e n u 1 to of particle angularity (PA) calculated on the 
above basis when plotted against the respective particle 
nir.es (Fig. 5 . la) revealed an inverse relationship. This is 
in agreement. with the observation by Twenhofel. (1950). 
However , the exact nature of such a relationship for any sand 
would depend on its grain morphology and the distance of 
sediment transport. The weighted average values of particle 
sphericity ( fj) ) » shape factor (SF), number of tangents and h/A 
for the sands investigated are presented in Table 5.2. 

5.1.3 Shape Factor Versus Sphericity 

On the basis of available .data a relationship between 
shape factor (SF) and sphericity ( ¥ ) is indicated (Fig. 
5.1b), Sands containing flaky grains show sphericity value of 
less than 0.68 and SF is 0.54, as compared to sands (most 
alluvial sands ) with bulky spherical grains having l|l >0.74 
and SF of about .65. For sands containing elongated to bulky 


grains , 


sphericity value varies between 0.68 to 0.74. 


This 
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Tabl 8 5,2 Ueight ed Average of Shape Factor, Sphericity, Number 
of Tangent, h/A and Roundness Index 


Sands 

Average 

Shape 

Factor 

Average 
Spher i - 
city 

Average 

Tangent 

Average 

h/A 

Average 

Roundness 

Standard 

0 .766 

0.810 

4 

0.0465 

0.78 

-0.80 

Kalpi 

0 .596 

0 .698 

12 

0 . 2511 

0.10 

-0 . 20 

Calcareous 

0.550 

0.670 

6 

0 .1218 

0.50 

-0 . 60 

Ganga 

0.539 

0 .678 

11 

0 .2740 

0.15 

-0 . 25 

Kalpi 1 

0 .616 

0 .710 

10 

0.3010 

0.10 

-0 . 20 

Kalpi 2 

0 . 578 

0.708 

13 

0.3770 

0.10 

-0 . 20 

Coarse 

Ganga 

0.553 

0.674 

6 

0 .1673 

0 . 40 

-0 . 50 

Artif icial 
Calcareous 

0 .676 

0.706 

8 

0.2281 

0.30 

-0 . 40 

Laguni 1 1 as 

94A - Ml 9 

0.625 

0 .693 

6 

0.1410 

0 .40 

-0 . 50 

Laguni 1 las 

94A - 94B 

0 .624 

0 .715 

6 

0 .1418 

0 . 40 

-0.50 

Leighton 

Buzzard 

0 .623 

0 . 689 

7 

0 .1788 

0 . 30 

-0.40 

San Fernando 
No. -5 

0.663 

0.722 

9 

0.2374 

0 .26 

-0.35 

San Fernando 
No. -6 

0 . 632 

0 .712 

9 

0.2112 

0.28 

-0 .38 

San Fernando 
No. -7 

0 .648 

0 .731 

8 

0.2190 

0.30 

-0 . 40 

Ottawa 

0 . 630 

0 . 709 

5 

0 . 0820 

0.50 

-0 . 60 
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data, a i eliable method is developed to compute the Powers’ 
roundnesn index which can be used to classify the sand 
according to their angularity (Table 2.2). 

5.1.5 Distribution of Tangents 

Distribution of tangents representing protrusions on 
the grain surface was studied in great detail. Data were 
obtained for different sieve .sizes and a weighted average 
distribution curve was evaluated for each sand. Figures 5.2 
to 5.5 show the results for Standard, Calcareous, Ganga and 
Kalpi nand respectively. 

In case of Standard sand (Fig. 5.2), more than 70% of 
the grains indicate 4 tangents (which is taken as a measure 
of perfect roundness in this study) . This in conjunction 
with the earlier observation that more than 80% of grains are 
spherical (Table 5.1) would suggest that Standard sand 
consists of smooth, rounded and spherical particles. In case 
of Calcareous sand , 47% of grains have smooth surface with 4 
and 5 tangents (Fig. 5.3). This finding along with the 
earlier observation (Table 5.1) would suggest that this sand 
consists of subrounded, smooth, flaky and oblate particles, 
which is a good description of the shells comprising this 
materials. In case of Ganga sand only around 12% of the 
particles are having smooth surfaces (Fig. 5.4), the 
remaining having number of tangents varying from 6 to 17; 
some particles have as high as 24 to 33 tangents. These 
results along with the findings related to sphericity and 
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shape factor would suggest that Ganga sand has angular, 

oblate and flaky grains . 

Kalpi sand as indicated in Chapter 3, is distinctly 
different from others in more. than one way in respect of 
short distance of sediment transport from the source and also 
the formation of carbonate on the sand grains by the 
circulating ground water. This sand shows quite an unique 
distribution of tangents with a tbtal absence of any grains 
with four tangents (absence of smooth grains). Around 45 % 
grains show 6 to 12 tangents and 10% of grains reflect as 
high as 12 to 16 tangents (Fig. 5.5). Two fractions of this 
sand, called Kalpi 1 and Kalpi 2, were also investigated for 
angularity characteristics. The results presented in Figs. 
5.6 and 5.7 indicate that roughly 50% of the grains show 14 
to 18 tangents for Kalpi 2 sand. 

These angularity characteristics along with the earlier 
results on sphericity and shape factor would indicate that 
Kalpi sand grain® are highly angular, bulky and elongated in 
nature . 

Results for two artificially prepared sands, viz. 
coarse Ganga and artificial Calcareous are shown in Figs. 5.8 
and 5.9. In case of coarse Ganga the major component is 

K,- 

Calcareous * and the tangent distribution shown in Fig. 5.8 
agrees with the pattern shown for Calcareous sand (Fig. 5.3). 
Similarly* the tangent count distribution for artificial 
calcareous sand is influenced by the patterns of its 
components* viz . * Standard* Kalpi and Ganga sands. 


% of Particles 
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Fig. 5.6 Tangent count distribution for Kalpi 1 sand. 
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Fig. 5.7 Tangent count distribution for Kalpi 2 sand . 


8 12 
Tangent Count 


For sand with all fractions 



Tangent Count 

Fig. 5.8 Tangent count distribution for Coarse Ganga sand. 
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Tangent Count 

Fig. 5,9 Tangent count distribution for Artificial calcareous sand. 



0.212-0.150 I For sand with all fractions 



Tangent Count . 

Fig. 5.10 Tangent count distribution for Ottawa sand. 



0 . 212 - 0.150 For sand with all fractions 16H- < 0.075 



Tangent Count 

Fig. 5.11 Tangent count distribution for Leighton Buzzard sand 
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Tangent Count 

Fig. 5.13 Tangent count distribution tor Laguniltas (94A + 94B). 
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sapjUDd 10 °/ 0 


Tangent Count 

Fig. 5.14 Tangent count distribution for San Fernando mix 5 sand 


0.15 - 0.075 
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Tangent Count 

Fig. 5.15 Tangent count distribution for San Fernando mix 4k 6 sand. 
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1. Subrounded to rounded grains 

2 . Subrounded to subangular 

3. Subangular to angular 

4. Very angular 


(Fig. 5 .17a) 
(Fig. 5 . 17b) 
(Fig. 5 . 18a) 
(Fig. 5. 18b) 


These distributions are suggested as typical examples 
of sands with varying degrees of angularity. It is further 
suggested that the nature of tangent count distribution 
pattern may be used for the purpose of classification of 
sands on the basis of shape characteristics. Since the 
tangent count is inversely related to the roundness index (R) 
(Fig. 5 . 1 e ) , it is also possible to construct similar 
distribution patterns for different sands, in terms of 
roundness index ( R ) . 


5.1.7 Main Observations 

Study on particle morphology reported here leads to the 

following observations: 

On the basis of Zingg diagram and the shape factor 
versus sphericity relationship, sands can be classified with 

respect to the grain shape. 

The tangent count technique using the image analyzer 
provides an accurate and highly reliable method to quantify 
the particle shape characteristics. Typical tangent 
distributions for subrounded to rounded, subangular, angular 
and very angular grains have been developed and are 
recommended as a useful tool for classification of sands in 
respect of grain angularity. 




Fig. 5.17 Typical distribution of surface characteristics for subrounded to 
rounded and subangular to subrounded sands. 
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Fig. 5.18 Typical distribution of surface characteristics for angular and 
very angular sands- : 
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The tangent count has been employed to provide a 
systematic method for evaluation; of Powers ' roundness index. 
An inverse relationship between • average tangent count and 
roundness index (R) has been established. 

It has also been shown that while average tangent count 
is useful for comparing sands having grains of different 
angularity, the particle angularity 

average number of tangents in a sieve fraction 

(P A = ) 

average circumference of the particles in that fraction 

is a good measure of variation of angularity with particle 

size for a given sand. Particle angularity for alluvial 

sands is shown to decrease with increase in grain size. 

On the basis of these observations it is recommended 

that in addition to the usual classification tests for sands, 

the shape and angularity characteristics of sand grains 

should also be evaluated using the techniques and methodology 

as detailed in the present study. 

5.2 PARTICLE MODIFICATION AND CRUSHING DURING COMPRESSION 

5.2.1 General 

Sands unlike clays undergo changes in grain shape and 
size during compression. The degree of alteration of the 
grains depends upon the mineralogy, size, gradation, 
angularity and stress level applied. While sufficient data 
on particle crushing under high stresses are available in 
literature (Hardin, 1985), very limited information is 
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reported on the degree of particle' modification. In the 
present study, experiments were planned to attain a better 
understanding of the factors controlling particle 
modification during compression. A wide range of sands with 
different mineralogi cal composition, particle shape, size and 
gradation were studied. In this study, the overall sand 
particle breakage under stress has been separated into two 
components, viz. grain modification and grain crushing 
(changes in grain shape and size). Grain modification 
envisages breakage of angular surface protrusions at stress 
levels far less than that required for grain crushing. Such 
a modification would result in smoothening of particle 
boundaries (reducing the number of tangents ) whereas 
subsequent grain crushing at higher stresses would again 
produce angular grain (increasing the number of tangents). 

5.2.2 Modifications of Size and Shape of Grains 

The potential for the grain modifications in terms of 
size and shape alterations including the crushing is best 
demonstrated through the study of changes in grain size 
distribution curves of a sand before and after testing. 
Hardin (1985) proposed the calculation of breakage potential 
(designated as Bp) for particle degradation on the basis of 
the area of grain size gradation curve for sizes greater than 
0.074 mm, (Fig. 5.19 a). This parameter reflects the total 
potential for change of grain shApe and size under stress for 
a given sand and includes both tfte potential modification and 



Unit area 



Particle size distribution 
100 r r\ v v v \ \ v — i 


074 mm .074 mm 

Particle Diameter (mm) 

Fig. 5, 19(a) Definition of breakage potential. 



1 log 
cycle 



.074 mm 


after 

loadinc 


original 


.074 mm 

Particle Diameter (mm) 


Fig. 5.19(b) Definition of total breakage. 
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crushing. The grain size curves for the sands studied 
(Fig. 5.20) and those reported in literature (Fig. 5. 21) have 
been used to estimate Bp. The values of Bp along 
with the average particle size and coefficient of uniformity 
values of the respective sands are. listed in Table 5.3. The 
relationship between Bp and the average grain size d ^ q has 
been developed ( Fig . 5.22) on the basis of these estimations. 
Data from literature have been used to extend ■ the range of 
d 5 Q values, to test the validity of the proposed 
relationship. The relationship, tested for sands with a wide 
range of gradation curves, indicates that B p is uniquely 
related to the average particle size (d 50 ) . It may be 
indicated her© that this relationship has not been tested for 
gap graded sands. In practice, for most of the sands, the 
value of Bp may therefore be estimated on the basis of d^Q 
(in mm) as obtained from the sieve analysis by using the 
expression : 

Bp = 1.02 log (d 50 ) + 1.15 (5.5) 

The above simplified expression enables an estimation of B p 
as a preliminary measure without undertaking the detailed 
area estimates of the grain size curves. 

5.2.3 Particle Breakage 

Hardin (1985) used a terra B x to define the difference 
in area under the gradation curves before and after loading 
(Fig. 5.19b) as a measure of total particle breakage 
comprising of changes in grain shape and/or shape and size, 
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Table 5.3 : Gradation, Size m n 

- . 2;r r - — — 


Sand n 

Kalpl 

Ka I p I 1 

Kalpi 2 

Gang a 

Ganga 
( enhanced 
mica) 


4.814 
l .073 
1.613 
2.565 
2 . 500 


1 - 00 
1 .02 
1 . 70 
0.18 
0.20 


DR B p B t Br Minerals (%) 

Quartz^Feldspar Mica Carbonati 
°- 25 1 ’ U 0.115 <0 

1 - 14 °- 18 0-158 <5 30 _ .. 

°' 25 °' 31 «-»0 BO 15-20 i -2 30 

°' 25 0-37 0.065 0.166 62-65 20-25 8 -10 2-5 
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Average Particle Size (dso) (mm) 

Fig. 5.22 Relationship between total potential crushing and average particle size-. 
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depending on the stress level. Vari.f 

Variations of values of B „. 
with stress level f nr .... T 

(or a variety of sands tested lt - 
, . . reseed in this study 

and for those reported by other workers a • , 

are floated in Flo 

5.23. The gradation curves before and aft 

„ . and aft6r loading for the 

sands investigated have alreadv v 

already been presented earlier (Figs 

4-59 - 4.67). The gradation curves for Ta 

3 f ° r Ta youra sand (Miura 

and Yamanouchi t 19 7 1 .) + 

'' 0tta “ a (Robert and De Sousa, 

58), Chattahoochee river sand (Vesic. ms , a „d Hoi sand 

COe Beer, before and after loading are given in Fifis . 

5.24 and 5 * 2 S * 


5.2.4 Variation of B T with Stress 

Plots Of B x versus stress (Fig 5 23) for- H - f f 

' ° f°r different 

sands have revealed the following interesting trends: 

1. Standard sand and „.i sand (Be Beer »«) composed of 
well-rounded quarts grains do not undergo any breakage 
0>T*°) at stress of 176 Kg/cm 2 in sharp contrast to 
Kalpi 2 sand (composed of highly angular, elongated 
grains coated with carbonate) which attains a value of B T 
around 0.31 for the same stress level. Thus the effect 
°( particle angularity, mineralogy, site and relative 
density on total breakage is clearly evidenced from the 
d.ita presented in Fig. 5.23. 


2- Sands consisting entirely of quartz, with either 
perfectly rounded, subrounded or subangular grains, 
exhibit a typical variation of Bj with confining stress. 
For example, at stress level of 984 Kg/cm 2 , the poorly 
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Fig. 5.24 Particle size distribution before and after load tor 

Mol sand De Beer (1963) & Tayoura sand Miura etal.(1973) 
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Fig. 5.25 Particle size distribution before and after load for 
Chattahoochi river sand (Vesic, 1968) and Ottawa 
sand (Robert & De Souza, 1958). 
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era ' l0d 0ttaua oand < C “ * 1.5. d 50 = 0.6 ™0 indicated 

relatively higher B T (=0.274) as compared to the value of 

By (-.14) for a graded Ottawa sand sample (Cu = 2.75, dgg 

= .475 nun) . ' 

Values of B T are also dependent on the grain size of 
sand. For example, at the same stress level of 984 Kg/cm 2 , 
Dol sand gives B x = 0.124 compared to Bj = 0.274 for 
uniform Ottawa sand . While Cu values for both uniform 
Ottawa (1.44) and Mol (1.5) sands are practically same, 
t. h e 0 1 1 awa sand is medium to coarse ( d ^ q = 0.6 mm) and 
Mol sand is fin® (d 50 = . 195 ) . ■ 

The layout' a sand has subrounded to ' subangular quartz 
grains and shows B x = 0.035 for relative density of 0.24 
compared to B T = 0.14 at relative density of 0.90 for 
the same stress level of 176 Kg/cm 2 . The Ganga sample 
with mica removed (composed mainly of quartz and 

feldspars) shows behaviour similar to the Tayoura sand, 

2 

e.g. at stress level of 176 Kg/cm , the B T value of 
Ganga with mica removed (at relative density of .40) fits 
very well in between the B T values for Tayoura sand at 
relative density of 0.24 and 0.90. These two sands have 
approximately the same d$ 0 ( - - 2 mm). The effect of 
initial relative density on B T is clearly exhibited by 
the loose and dense samples of Tayoura sand. 

3. For the same mineralogy and degree of particle angularity 
( in 


case of crushed granite as used by Lee 


and 
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Fat hoomand , 1967), the value of B T increases with d 50 for 
the n&mti confining stress (points A, B, C, D and E in 

Fig. 5.23). Results for Napa basalt (Karachi et al . , 
1969) also depict similar behaviour. Identical behaviour 
io also illustrated by the data for Kalpi, Kalpi 1 and 
Kalpi 2 samples with d 5Q values of 1.0, 1.02 and 1.7 mm 
respectively. However, the disproportionately higher 
response observed in the case of Kalpi 2 sand can be 
attributed to the uniform and coarse nature of the grains, 
as well as the abundant carbonate coating on the same. 
Influence of d^Q on B-p is further illustrated by the B*p 
vermin stress behaviour of micaceous Ganga sand having 
the name degree of particle angularity as exhibited by 
Kalpi sand with 20% carbonate coating.' The smaller 
values of B^ at the same stress level in case of Ganga 
sand can only be attributed to the fact that d 50 (around 
0.18 mm) of Ganga sand is less than that for Kalpi sand 
( d r, q a 1.0 mm) . 

4. For the same d 50 (0.2 mm) and comparable relative 
density, Ganga sand having 8 - 10% mica has indicated 
higher B-p (0.0625) as compared to Tayoura sand and Ganga 
sand (with mica removed) at the same stress level of 176 

Kg/Cm 2 . 

5. In case of Ganga sand, at the same stress level (176.32 
kg/ cm 2 the value of B-p increases with the increase in 

mica content (points Gq, G and 


in Fig . 5.23). 
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Incidentally all three samples have same d 50 and similar 

gradation curves. 

6. Calcareous sand, composed of 9 5?; of carbonate material, 
with nubt ounded to subangular grains exhibited variations 
of with stress corresponding to those between 

micaceous angular Ganga sand and highly angular Kalpi 
sand with 2 0% carbonate coating. .It may be recalled that 
the Bp values for Ganga, Cacareous and Kalpi sands show a 
similar trend (0.39, 0.77 and 1.11 respectively). 

Particle breakage response of; different sands can be best 
expressed in terms of relative breakage B r = (Bp/Bp , the 
details of which are presented in the subsequent pages. 

It can thus be concluded tjtiat for the same mineralogy 
and particle angu 1 ar i ty , increases with dg 0 and decreases 
with C u . Further for the same d 50 increase in mica and 
carbonate contents results in the increase of B T . 

Except for Ganga sand, B T for all the sands tested in 

this study increases with grain angularity (h/ A ) at low 

2 

strews levels. Data at 22.04 Kg/cm is shown in Fig. 5.26a. 
The Ganga sand behaviour is perhaps due to the fact that the 
grain size (d 50 ) = 0.18 mm) is very fine, whereas for all 
other aands, the d^Q values vary from 0.4 to 1.7 mm (Fig. 
5.26a). Had Ganga grain size been coarser ( = 1.0 mm) its 
plot would have followed the trend depicted by other sands. 
Variation of B T with average number of tangent for these 
sands would give a similar trend. 
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Fig. 5.26(a) Relationship between B T and grain angularity (h/A) 

for different sands at RD of 0.25 and stress 22.04 kg/cm 2 . 



193 


5.2.5 Variation of B T with d 5Q 

Depend "nc« of B ? on d 50 f or different aanda is 

illustrated in Fig. 5.26b. Here the effect of particle 

angularity , mineralogy and stress level are brought out. At 

the same stress level, samples with same C u , grain 

mineralogy, and angularity (in case of particles for crushed 

granite) show consistent increase in B T with d 5Q . Similar 

trendn are shown by well-rounded quartz sands (Ottawa and Mol 

sands). Angular Kalpi and Napa basalt samples also exhibit 

an increase in B T with d 50 but the rate of increase is 

greater compared to the rounded and subrounded to subangular 

grains . lor loose samples of- micaceous Ganga sand and 

Tayoura sand with similar values of d 5Q , the value of B T for 

Ganga sand at 176 Kg /cm 2 is about 1.9 times greater than that 

of Tayoura sand at the same stress. The effect of initial 

relative density on B T in case of Tayoura sand (loose = 

0.24, dense = 0.9) is less than 20%.- However, the effect of 

stress level for the same particle angularity and relative 

density is much more pronounced. B T for loose Tayoura sand 

2 

sample increases from 0.035 at 176 Kg/cm to 0.095 at 450 
Kg/ cm 2 . It may be noted that as the stress increases by a 
factor of about 2.6, the corresponding increase in Bj is of 
the same order (2.7 times). From Fig. 5.26b, it may be 
concluded that as in case of Tayoura sand, for the same 
mineralogy, particle size (d 5 0 ) , and grain angularity, stress 
level controls the magnitude of particle breakage. Also for 
the same d$ 0 and at comparable relative density, sands of 
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iq. 5.26(b) Relationship between average particle size (dso) and total breakage (Bf) 
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similar mineralogy show increase in Bj with particle 
angularity when Mol and Tayoura sands in dense state are 

compared . 

5.2.6 Particle Modification 

The earlier discussion on particle breakage in terms of 
®T did not distinguish between two different phenomena 
leading to particle degradation. They are: 

1. Grain modification at low stress levels, when the 
grain boundary protrusions • are knocked off without 
any grain splitting. ■ This leads to reduction in 
particle angularity, 

2. Grain splitting at high stresses which would again 

produce angular grains .. 

There is liar d 1 y any data available on .grain modification. 
The basic questions to be answered are: 

1. A what stress level does the grain modification 
start for different sands , 

2. Upto what stress level does grain modification 

continue? 

Standard Ganga, Kalpi and Calcareous sands were 

investigated for grain modification. Since loose samples 

would experience more grain degradation, samples of these 

sands at relative density of 0.25 were subjected to 

compression under stress levels of 11.02, 22.04, 88.16 and 

176.32 Kg/cm 2 . At the end of each stress level the samples 

were taken out and were screened for grain size distribution. 

2 

Except for the samples tested at 11.02 Kg/cm 


studies for 
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surface modification were undertaken for all others using 

image analyzer. 

5. 2. 6.1 Changes in Grain Shape 

Based on flatness (dg/dj) and elongation (dj/dj^) ratios 

the changes in grain shape due to particle degradation under 

stress are best illustrated by the Zingg diagram with 

sphericity and shape factor superimposed. The results of 

experiments for Standard sand, Ganga Sand, Calcareous Sand 

and Kalpi Sand are presented in Fig. 5.27 ( a, b, c and d) . 

In the case of Standard sand where no particle breakage was 

observed even upto stress level of 176.32 kg/cm 2 , sieve 

analysis at all stress levels indicates no changes in grain 

size distribution. The measurement of particle dimensions 

(dg, dj and dj^) revealed no alteration in either flatness or 

elongation ratios. The initially rounded grains remained so 

2 

even after loading upto 176.32 Kg/cm . However, in the case 
of Ganga Sand (Fig. 5.27b), the initial average point 1, 

representing the original particle shape changes to position 
2 after loading to 22.04 Kg/cm 2 . Point 2 in quadrant II 
indicates an increase in particle shape factor and 
sphericity. This is due to the modification of the grain 
surface protrusion which reduces the particle angularity. 

41 

Further loading to 176.32 Kg/cm 2 leads to particle rupture 
for this sand resulting again in the production of grains 
with increased surface protrusions as indicated by point 3 
located in quadrant I once again. Position of points marked 



(a) Standard 


Stress 
° O.OOkg/cm 2 
A 22.04 kg /cm 2 
° 176.32kg/cm 2 


(b) Ganga 


© > Average points 
Bj 


.8. .9 .SF 




(c) Calcareous 


(d) Kalpi 





Fig. 5.27 Relationship beetween flatness and elongation ratios showing 
the particle degradation under stress for Standard, Ganga, 
Calcareous and Kalpi sands. 







198 


1,2,3 for Calcareous and Kalpi sands illustrate similar 
changes in particle shape for stress levels of 22.04 and 
176.32 Kg/cm 2 respectively. While in case of Ganga sand the 
final point 3 moves towards quadrant I, the final grain shape 
in case of Calcareous and Kalpi tends towards elongated 
grains represented by quadrant III. 

5. 2. 6. 2 Changes in Grain Angularity 

Based on tangent count, the distribution of tangents at 
stress levels of 22.04 and 176.32 Kg/cm 2 as well as original 
samples for Ganga , Calcareous and Kalpi sands are presented 
in Figs. 5.28, 5.29 and 5.30 respectively. Here again at 
stress levels of 22.04 Kg/ cm 2 , the number of tangents 
decreases as the surface protrusions are modified and the 

grains experience reduction in surface angularity. Further 

2 

loading to 176.32 Kg/cm leads to grain splitting which in 

turn leads to increased surface protrusions and increased 

tangent count. This evidence also supports the fact that at 

2 

a stress level of 22.04 Kg/cm , the grains had experienced 
surface modifications (as evidence from the lesser number of 
tangent counts) which resulted in lower particle angularity 
and at stress level of 176.32 Kg/cm crushing of the grains 
once again resulted in the increase in their angularity as 
observed from the increase in tangent counts. The Standard 
sand shows no change in tangent distributions. The tangent 
count data were also used to interpret the changes in 
percentage of particles having edges with one tangent count 
only. All tangent counts were, carried out on the four 
projected edges of the grain. All grain edges (1 to 4) with 





8 


Tangents 


Fig. 5.29 Distribution of tangents at <j'= 0.0, 22.04 and 
176.32 kg /cm 2 far Calcareous sand. 





201 


one tangent only are designated as smooth in Figs. 5.31 to 
5.33. The data for smooth edges for original samples and 

after subjected to stress levels of 0.0, 22.04 and 176.32 

Kg/cw for total sample and also for the different seive 

fractions in case of Ganga, Calcareous and Kalpi are 
presented in Figs. 5.31, 5.32 and 5.33 respectively. Here 
again at 22.04 Kg/cm , the number of particles with smooth 
edges increases (angularity decreases) and when grains are 
fractured at 176.32 Kg/cm^' the number of particles with 
smooth edges is once again reduced. These observations are 
complimentary to the data on tangent distribution count 
(Figs. 5.28, 5.29 and 5.30). 

Since h/ A 1 a a direct measure of grain angularity, 
values of average tangent have boon plotted against h/A (Fig. 
5.34). The data for Ganga, Calcareous and Kalpi sands at 
stress levels of 22.04 and 176.32 Kg/cm are also shown in 
Fig. 5.34. For a variety of original sand samples, a very 

good correlation exists between the number of tangents and 

2 

h/A. However, at a stress level of 22.04 Kg/cm as the number 

of tangents decreases, the grain angularity represented by h/A 

2 

also decreases and on grain crushing at 176.32 Kg/cm the 

data point tends towards higher initial tangent count and 

h/A . For example 1,2,3 represent loading stages for Ganga 

2 

under stress level of 0.0, 22.04 and 176.32 Kg/cm 
respectively. For three sands, viz., Ganga, Calcareous and 

Kalpi, the h/A versus tangents relationship shifts from 

1 

curve A to B and to C (close to A) as stress level increases 


of Particles 






»I 1 — 1 _ J I I i 1 1 L 

0 12 3 4 

Smooth Edges 

Fig. 5,31 The smooth edge distribution at d' = 0.0 , 22.04 and 
176,32 kg/cm 2 for Gcmga sand. 


/o of Particles % of Particles % of Particles 

















Original 
22,04 kg /cm 2 



Lq Leighton Buzzard 
O Ottawa 

Lagunillas (94A 94B) 
Lagunillas (94A-19M) 
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from 0.0 to 22.04 and 176.32 Kg/cm 2 . 

Thus the present investicAt ;«««, +u ■ 

cstigations on the grain shape and 

angularity have clearly brought out the fact that prior to 

the giain crushing at very high stresses, there exists a 

strum, level (^crushing stress) at which the grains undergo 

surface modifications. These surface modifications, as 

detailed in succeeding pages, exert a significant influence 

on the variation of constrained modulus with confining 

stress . Grain modification at a certain stress level 

( << crush i ng stress) under cyclic triaxial undrained loading 

leading to change from dilatant to compressive response of 

the sands has also been reported by Vaid and Chem (1985) who 

have suggested the dependence of this stress level on the 

initial relative density for angular mine tailing sand. 

Recognizing the grain modification at lower stresses, a 

detailed study was undertaken for all the sands tested under 

2 

stress levels of 11.02, 22.04 and 88.16 Kg/cm for 

establishing the threshold critical stress at which the grain 
modification begins . Since grain degradation is expected to 
be maximum under loose packing, all the tests were carried 


out on samples with initial relative density of 0.25 exccept 
for Standard sand for which the relative density used was 

i 

0.20. The results of this investigation are presented in 
Fig. 5.35 where the relative particle breakage B r (calculated 


as a ratio of B-j- and Bp) is plotted against the measure of 

2 

particle angularity h/A for data at 11.02 Kg/cm , 22.04 


Kg/cm' 


Kg/cm 2 stress levels 


11.02 Kg/ cm 


evidenced in Calcareous and Kalpi 2 


significant changes were 




Fig. 5.35 Relationship between the relative breakage and h/A for variety of 
sands . 
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a audit while Kalpi sand showed no grain modification. Major 
increase in B r takes place under stress level of 22.04 
Kfi/cm . lor stress increase from 22.04 to 88.16 Kg/Cm 2 , most 
sands except Kalpi 2 show minimal increase in B r . The 
behaviour of Kalpi 2 at 88.16 Kg/Cm 2 may be due to the early 
onset of grain crushing since this sand, as indicated 
earlier, has uniform grading, coarser grains and higher 
carbonate content. 

These results thus indicate that significant grain 

modification can occur even at a stress level less than 22.04 

2 

Kg/cm ' . Ait described in a subsequent section, the variation 
of modulus with stress has been used in the present study to 
determine the critical stress level at which the behaviour of 
the sand changes. This stress level is related to the stage 
at which grain modification also commences. A tentative 
relationship has also been proposed between this critical 
stress causing grain modification and relative density 
(Sec t i on 5.3.2). 

5.2.7 Particle Crushing 

As indicated by Marsal (1967), the phenomenon of 
particle fragmentation under stress, both during compression 
and shear, a very significant influence on the shear strength 
and compressibility of granular materials. Hardin (1985) has 
proposed a parameter called relative breakage (B r ) , to 
quantify the degree of particle; fragmentation or crushing, 
defined (Fig. 5.19) as the ratio of total particle breakage 
( B<£ ) and total particle potential breakage (Bp). In the 
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»ludy, all co,r, ssiM t.«t. ware oni -dimensional. 
Data from literature used for comparison also pertained to 

the cowin' fifiaion stage only. 

In the previous section, various aspects related to 
particle modification at low-stress levels were discussed. 
The overall variation of B r with stress is discussed in the 

following aections in some detail . 

5. 2. 7.1 Variation of B r with Stress 

In i i g . 5.36, experimental data for the sands investi- 
gated in thin study as well as for Tayoura sand (tliura and 
Yamanouchi, 1973), Ottawa Sand (Robert and De Souza, 1958), 
Chattahoochee river sand (Vosic 1 978), Mol sand (De Beer, 
1963) and Crushed Granite (Lee and Farhoomand, 1967), are 
presented. The same data alongwith Napa Basalts (Karachi et 
al , 1 969 ) has been shown on a semi-log graph (Fig. 5.37) in 

order to better illustrate the trends over a wide range of 
stresses. In the following discussions, the rounded to 
subrounded quartz sands of low to very low compressibility 
will be broadly referred to as strong sands. The angular to 
very angular, compressibi le to very compressible sands 
containing mica and carbonate will be called as weak sands. 
Angular rockfill falls in between these two categories. This 
categorization is more for the purpose of simplified 

presentation rather than a strict classification. From Figs. 
5.36 , 5.37 , 5.38 and 5.39 the following observation can be 


made : 



Ganga 
Calcareous 
Kalpi 
Kalpi i 
Kalpi 2 

Tayoura Loose 
Tayoura Dense 
Mol Sand 


w ^anga enhanced mica 
-f Ganga mica removed 
® Crushed Granite 
* Ottawa 

Chattahoochee 


CT (kg/cm 2 ) 

Fig. 5.36 Relationship between relative breakage and vertical 

effective stress. 





Calcareous 
Ganga 
Kalpl 
Kalpl i 

Kalpl 2 
Tayoura loose 
Ottawa 

Tayoura dense 
Ganga dense 
Kalpl dense 


Model Prediction 
Loose condition 
for Ganga sand^/ 


CT (kg/cm 2 ) 

Fig. 5.38 Variation of relative breakage at low stress and model 

prediction . 
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Fig. 5.39 Variation of relative breakage at high stresses. 
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1 . In cants of well rounded strong sands (e.g. Mol sand and 
ii*.UH3.ud sand) upto a stress level of 176 . 32 to 200 
Kt/'-'in , practically no grain degradation takes place 

b . 36 , 5.37), Ottawa sand (Robert and De Sourza, 

1 9bB) which is leatJ mounded than Mol and Standard sands 
shown B f - 0.01 at 100 Kg/cm 2 and 0.03 at 200 Kg/Cm 2 . The 
value of B r at stress Irvel of 500 Kg/Cm 2 is comparable 
with of weak sand at 100 Kg/Cm 2 . At stress levels greater 
than 500 Kg/cm 2 , B r increases at faster rate (0 . 0005 
Kg/cm ) • 1 he increased rate of change of B r for strong 

Hands in about 1/2 of that for weak sands. 

2. For v**«k Hands (such aa Ganga and Kalpi), values of B r 
b«tu««n 0.07 and 0.12 are obtained at a stress level of 
1 00 Kg/ cm 2 and roughly 2/3 of the grain modification has 
already taken place at a stress level of 20 Kg/cm . This 
fact wan used earlier to confirm the grain modif icat ion 
phenomenon at low stresses. At a stress level greater 
than 100 Kg/cm 2 , the values of B r increases at a faster 
rate (0.0012 Kg/cm 2 ) . 

3. The trend a of the subangular to subrounded sands (Tayoura 
sand) are intermediate between the same for weak and 

utrong aaniJii. 

4. Angular rockfill behaves almost similar to the weak sands 
at low stress levels (see Napa Basalt Fig. 5.37). At 20 
Kg/cm 2 the a lightly lower value of B r is due to the fact 
that this rock fill is stronger than the micaceous Ganga 
Hand and weakly cemented Kalpi sand. At stresses greater 
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' ,0 hi “ h Of rocMill , iT „ greater 

v.U.,.. of B r a. compand to o.al sands . Napa basalt 

t of. al , 1971), has d 5 Q ranging from 2-10 nun and 

uo " 1 '' higher B r at higher stresses. It may 

b * in t.r.ati„ e 4 ° “hserve (Table 5.4) that at the same 

0 1 r *» n a both B„ and Bm inrrMoo W r+. K • . ^ . . 

P D T in( -roase with increasing particle 

” izr; (d 50)- whereas the relative breakage B r . remains 

practically same for a given material. In case of Kalpi 

sand, B r value for Kalpi 2 is higher than that of Kalpi. 

This in because of two factors: 

i. Kalpi is well-graded (cu=4 .8) compared to Kalpi 2 (Cu 

* 3 .«). 


2 . Kalpi 2 consists of coarser fraction where the small 
grains have been cemented together to form the larger 

aggregate grains . 

From this limited information it may be possible to 
suggest that, except for the case where aggregated grains are 
involved, both Bp and B-p for a material of similar 
grad.it ion, angularity and mineralogy,, increase with particle 
size whereas B r at a given stress remains essentially the 
name. These observations are clearly supported by the trends 
shown in Figs 5.23 and 5.37. In Fig. 5.23, the Bj versus 
stress relationships are distinctive for each d 5Q in case of 
Kalpi, Napa basalt and crushed granite whereas the same data 
when normalized with respect t o • Bp shows almost identical 
behaviour for both particle sizes (Fig. 5.37). From Fig. 5.37 
it can a 1 ii o be seen that loose materials give higher Bj. for 
the b am e stress as compared to dense materials (Ganga, Kalpi 
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Table 5.4 


Values of den R u . n 

50' B p • and B r for 

at Effective Stress of 80 Kg/cm 2 


Different Sand 


Sands 

Particle 

shape 

d 5Q 

mm 

B P 

B *jr* 

B r 

Stress 

level 


~ 







Napa- 
Basalt 1 

Angular 

10.0 

2.16 

0.27 

0 .125 

80 

Kg/cm 2 

Napa- 
Basalt 2 

Angular 

2 . 3 

1.52 

0.19 

0 .125 

II 

II 

Crushed 1 
Granite 

Subrounded 
to Angular 

10.0 

2 . 16 

0 .18 

0.083 

W 

If 

Crushed 2 
Granite 

Subrounded 
to Angular 

2.0 

1.50 

0.14 

0 .093 

n 

If 

Kalpi 2 

Very Angular 

1 . 7 

1 .40 

0.14 

0 . 100 

n 

If 

Kalpi 

Very Angular 

1.0 

1 .10 

0.07 

0 .064 

if 

If 
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and Tayoura sands). The effect of particle gradation on B r 

has al t oady been pointed out in case of Kalpi and Kalpi 2 

sands i.e. for sands of same mineralogy and angularity, a 

poor 1 y ~ ad ed sand gives higher B r . compared to a well-graded 

sand. The effect of mica content on B r at same stress level 

is well bi ought out by Ganga sand with mica removed, Ganga 

sand and Ganga sand with enhanced mica shown as 1,2,3 

respectively in Fig. 5.37. Incidentally all the three 

samples of Gang a sand had identical gradation, d 50 and 

comparable relative density. The value of 'B r . increases as 

mica content increases (Fig. 5.40a). Test data point for 

mica-rich Chattachoochee river sand (Vesic, 1968) is also 

shown in Fig. 5.40a. The higher value of B r for 

Chattahoochee sand is due to the fact that it was tested at 
2 ? 

647 Kg/cm compared to 176.32 Kg/cm in case of Ganga sand. 
Also the d^Q of Chattahoochee river sand is 0.38 mm compared 

to 0.18 for Ganga sand. 

2 

The effect of mineralogy on the value of B r at 176.32 Kg/cm 
has also bo examined in terms of quartz or quartz plus 
feldspar content of sands . These relationships are indicated 
in Fig. 5.40 b for a variety of sands. An inverse linear 
relationship (which may be fortuitous) is evidenced between 
B r and quartz content. Of course, a relationship of this 
nature is expected. The relationship between B r and quartz 
plus feldspar is shown in Fig. 5.40 b. It may be noted that 
the Ganga Sand data point (Fig. 5.40b) lies closer to the B r 
versus quartz plus feldspar content curve than the B r 


versus 



% 




Fig. 5.40 Variation of Br with mineral content of sands. 
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quartz content relationshio r:+- 

may thus be noted that the 

B r Vernon ,««r.« content relationship is for non . 

, F s tor n °n micaceous 

sand a only. 

Baaed on the data for r a i« 

Calcareous, and Kalpi sands, the effect 

of carbonate material on partin* 

Particle crushing has been evaluated 

CFi *- 5 ' 40 c >- Ganga sand (with 

C Xh mica removed) has 4-5% 

carbonate content and the , • 4 . 

h« the data fits veil with Kalpi sand 

Calcareous . A d i *,+ ; + 

custmct increase of B wi+K • 

B r w ith increasing 

percentage of carhonat#* i. , 

oonate content is evidenced. 

These results clearly reflect „ 

y reflect the markedly different particle 

crushing behaviour of weak 8 a „H 0 

weak sands, as compared to strong 

quartz sand. The effect , 

effect of mineralogy, gradation and 

£l> .IT t» Jl €J I* ■© 1 1 1, Z ft l'i m #% t . a , 

has been brou « h t out and tentative 

relationships , based on available data v, , 

<Avai iaoi e data, have been suggested 

(Fig. 5 . 40a , b, c) . The : • 

’ j lhe Variati ° n of Particle crushing with 

stress is in general shovn to follow the pattern suggested by 
Terzaghi and Peck (1943) and De Beer (1963). However, the 
weak sand. exhibit a distinctly different relationship as 

compared to strong quartz sands for stress levels less than 

100 Kg/cm I. 2 . 

1 a caao of strong quartz sands (Fig. 5.39) (well- 
rounded Standard and Mol sand), the variation of B r with 

stress occurs in four phases : 

I. For stresses upto 200 Kg/cm 2 , there is no 

degradation of particles and B r is almost 
(close to the zero value). 
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II. For stresses between 200 to 400 Kg/cm 2 , B r 

increases relatively faster than that in phase I. 

In thia Phase grain crushing has started. 

Ill* A i ci jp i cl v a f 0 o£ incr a nf r _ _ r . 

wi. xuli ease or b r occurs for stresses 

between 400 to 1000 Kg /cm 2 . 

IV. For stress greater than a limiting stress of 1000 

2 

Kg/cm ' rate of increase of B r is relatively less 

than that in Phase III. 

For weak sands (Ganga, Kalpi, Calcareous and also crushed 
Granite and Napa basalt), the change of Bj.. with stress 
(Fig. 5.3?) may also be described as: 

i) For stresses less than 20 Kg/cm 2 , the increase of 
& r is considerable and particles undergo 
mod 1 f i cat ion . 

O 

ii) For stresses between 20 1;o 100 Kg/cm ' the 

modification continues but B r increases at a lower 

rate compared to the rate in phase ( i ) . 

2 

iii) For stresses greater than 100 Kg/cm , significant 
particle crushing begins and a significant 
enhancement in the rate of B r is evidenced upto a 
limiting stress which is likely tp be smaller than 

1000 Kg/ cm 2 . 

iv) For stresses greater than the limiting value, 
increase of B r would be at a lower rate than that 
in phase (iii). A measure of the magnitude of 
limiting stress upto which B r increases at higher 
rate and also the decrease in the rate of B^ at 
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A 


hifillAr than i ’ . , • 

g 1 tnan the limiting stress, are best 

illustrated bv thA + c * t 

uy me data for Napa basalt as 

i ml i cat od in Fig. 5.41. Here instead of ( j\ , a 
stress function C\, called the breakage 
effective stress (Hardin, 1985), has been used. 

tf'b = G'o U + 9( Tq/ g' 0 ) 3 ] (5.6) 

where n\ + 2(J\ 

Go = 

3 

, G'l - O ' 3 

T 0 = 

2 

It will be seen that in case of Karachi’s data, 
the enhanced rate of increase of B r is upto a 
limiting stress ((J b ) of 250 Kg/ cm^ as compared to 
G b = 1150 Kg/cm for strong Hoi sand. Further 
the increase of B r at a much smaller rate after 
the limiting stress are clearly noticed from the 
data of Karachi (1975) and Leslie (1975). 

The behaviour of subangular to subrounded sands 
confotmed to the range intermediate between that of weak and 


strong sands described above. 

5.2.8 Hardin's flodel for Particle Crushing 

The hyperbolic relationship (Hardin, 1985) between 
relative crushing ratio (B r ) and stress ratio of breakage 
effective stress a' b and breakage reference effective stress, 
G' r ’ is expressed as 

ca’ h /a- r ) n b 

pt a ' ■ — , C 5 . 7 J 

" iH0\/0' r ) n b 

wher# i & called the breakage number. 




Fig. 5.41 Relationship between relative particle breakage and breakage effective stress 
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Baaed on the analysis of available data on crushing for 
variety of sands, Hardin proposed a relationship between (j' c 

and n b as 

(j r = 800 P a (n b - 0.3) (5.8) 

C c has been shown to be related to particle shape number 
(n s ), crushing hardness (h), and initial void ratio ( e ^ ) . 
Using the values of (J / r * at different e ^ , Hardin proposed that 

( 1 + ei ) 

— = H (constant) (5.9) 

800 P a 

From the values of H, and the definition of . particle shape 
number ( n {} ) the particle crushing hardness (h) is estimated 

using the expression 

, 2 
h 

H ~ (5.10) 

n s 

The computed value of h is approximately equal to the scratch 
hardness (Ploh’s Scale). The value of h and n s for a variety 
of sands are reported by Hardin. In order to establish a 
stress path independent relationship between B r and effective 
stress, the breakage effective stress, Cf b , was defined by 
Hardin as 

cr' b = cr'oCi+n-r'o/cr'o) 3 ) (5.6) 


where 



is a first 
Cf\ + 2 



stress invariant, 

>nd 


2 
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Accord i ng to Hardin, a ' b normalized the relative breakage 
data from tests with varying values of c\/ (j'3 - R and to 
illustrate the above relationship crushed granite data of Lee 
and Farhoomand ( 1967) was used. The relationship of B with 
tf'b and °'l for crushed granite and Napa basalt are indicated 
in Figs . 5.42(a, b) and 5.43(a, b) . The data of crushed 
granite show that scatter of data point for different 
R( Cf'i/ 0'$) values is of the same order in Fig. 5.42(a) and 
5.42(b). All the data points for Napa basalt (Hardin 1985) 
have been plotted in terms of B r versus 0' b /P a (Fig. 5.43a) 
and B r versus 0 ( j /P a (Fig. 5.43b). Hence again there is no 
unique normalization with respect to R in either plot (spot R 
values indicated). These comparisons tend to suggest that 
both 0 ^ and 0 b produce more or less identical normalization 
of teat data with varying values of R-G'i/g'z* It is thus 
evident that Hardin’s particle crushing law (equation 5.7) 
could as well be stated in terms of g\ instead of 0^ . The 
advantages of such a formulation are obvious . 

Incidentally, it is interesting to examine the (j' h / 0 / 1 

ratio for a wide range of R= 0 1 / O 1 2 oc K- 0 3 / CF 


1 + 2K 
3 

1 + 2K 


1-K 


) [1+9 ( 3 / 2 ( ■ 


)) 3 ] 


1 + 2K 


1-K 


or 


0 1 C- 


0 b 


•) (1 + 30. 4(- 


) 3 ) 


1 + 2K 


1 + 2K 
1-K 


(1 + 30. 4(- 


> 3 ) 


3 


1 + 2K 


(5.11) 






ig. 5.43(a) Relationship between Op/Pa and relative 
rea age (Based on data from Hardin, 1985) 
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Fig. 5.43(b) Relationship betvyeen / Fti and relative 

breakage (Based on data from Hardin, 1985) 
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For values of K ranging from p.l (R=10) to 1.0 (R=1.0), 

O' b/ a 'l VAluftS are Plotted (Fig. 5.44). The relationship 
suggests a b / o\ = 5.5 for K = 0 . 1 (R=10) and 0,' h J o\ = 1.0 for 

K=0 .405 (R=2.06) and K=1 . 0 (R=10). . 

In case of Napa basalt, K Q under one-dimensional 
compression is 0.34 and the average value of K Q under plane 
strain is between 0.22 to 0.26. From Fig. 5.44, it can be 
seen that for compression under K Q and plane strain 

condition, the values of (j' b / q ' are 1.5 and 2.3 

respectively. Under isotropic compression the value of 
is equal to one. Under condition of one-dimensional 
compression of sands, where K Q under normally consolidated 
state ranges between 0.33 to 0.55 (Mayne and Kulhawy, 1982), 

/ f 

the values of 0 b / fl would vary between 1.65 to 0.95 

respectively. Therefore, as a first approximation, the value 
of Cf (j\ could be taken as unity for estimating the 
crushing behaviour of sands under one-dimensional compression 
using the guidelines for h,n a and o' c as proposed by Hardin 
(1985) . 

However, as shown in Figs. 5.4 2 and 5.43 the model 

parameters can as well be worked out in terms of O' ^ rather 

than 0 b in which case no approximations are necessary. The 
main difficulty in 0 b parameter is the necessity of obtaining 
reliable values of K 0 in one-dimensional and plane strain 
conditions of compression. In case of triaxial anisotropic 
compression 0^ can be readily estimated. 

Hardin’s particle crushing law (equation 5.7) assumes a 
hyperbolic relationship between B r and stress (0 b or 0^). 
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However, the proposed relationship has , two material 
parameters! (0 r and n b ) which, if varying over a wide range, 
ma y control the shape of the B r versus stress relationship as 
evident from Fig 5.45. The values of a' r and n b have been 
chosen for the extreme cases of weak and strong sands. It is 
clear that the form of B r versus stress relationship varies 

I 

with Cf r and values. The characteristic form as proposed 
by Hardin is valid for all sands having c' v between 300 to 
1000 kg/cm ( h- 6.3 to 6.5). In particular, this is true for 
angular weak sands. Available data for a variety of sands is 

also shown 1 n Fig. 5.45. 

Although at the onset, it appears tempting to conclude 
that within the scatter range the experimental data 3 eems to 
fit the computed relationship based on equation 5.7, it would 
be instructive to observe the rate of change of B r with 
stress for the weak as well as the strong sands and then 
compare these trends with the experimental observations made 
earl i or . 

A 

In case of weak sands ( 0 ^ = 300 - 1000 Kg/cm ), the 
model indicates initial higher rate of change of B r upto a 
limiting stress after which the rate decreases significantly. 
The limiting stress increases with increasing 0 r . However, 
in case of strong sands, a low initial rate of change of B r 
followed by an increased rate at high stresses is noticed. 
The experimental data for Napa basalt (Karachi et al 1969) 
and Mol sand (De Beer 1963) contradict this predicted trend 
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(Fig. b.ib). For all materials plotted here, either a' h was 
available in literature (Napa basalt) or was computed on the 


basis of available information 


Tq/ 0^0 anc * 0 o foe Hoi and 


Ottawa sands. The data for Tayoura sand is for isotropic 

consolidation and hence 0^ = (j'^. 

In case of Napa basalt, the data points plot well on 
the predicted curve upto 250 Kg/cm 2 after which the data 
transcends towards the predicted curve for a stronger sand. 
The comparison between data and the model prediction is 

illustrated in Fig. 5.46. The model would thus predict Bj, 

? 

well upto 2 5 0 Kg/cm for Napa baslat data which was used to 
calibrate the model. At higher stresses, the model however, 
would considerably over predict B r values. A similar trend 
is observed in case of Hoi sand (Fig. 5.45). The 
experimental evidence reviewed in Figs. 5.36, 5.37, 5.38, 

5.41 and 5.45 also suggest that a value of B r of around 0.25 

2 

will be attained at a stress levels of the order 200Kg/cm 

2 

and 900 Kg/cm' in case of weak sands and strong sands 


respectively. Therefore , it would seem that within 100 

7 

Kg/cm maximum B r values that could be attained for weak 
sands are of the order 0.15 and for strong sands B r =0. As 
already indicated values of B r between 0.25 to 0.3 may only 
be obtained at much higher stresses. In practice in the 


range of stresses encountered, the B r value is thus likely to 
be less than 0.5. The model prediction would estimate 
tf b =500 - 1500 Kg/cm 2 (weak to strong) for B r =0.5 and (J / h > 
10,000 Kg/cm 2 for B r =1.0 for all types of sands (Fig. 5.47). 
Clearly such high stresses needed to get values of B r > 0.5 
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are not 


likely to be encountered in routine 


engineering 


prac tire. 

it would also seem that in case of weak sands, grain 
crushing is of major concern in the study of strength and 
compressibility behaviour in the range of stresses 
encountered in case of deep foundation or during cone 
penetration testing. However, in case of strong sands, like 
Hoi and standard, grain crushing becomes important for 
stresses greater than 300 Kg/cm 2 . For subrounded to 
aubangu 1 at sands which fall in between weak and strong sands, 

there is likely to be significant grain crushing at stress 

2 

level upto 300 Kg/ cm . 

Based on these results and comparisons with prediction, 
the following comments on the Hardin’s model are pertinent: 

1. In case of prediction of relative breakage under stress 
based on the hyperbolic law, the shape of the predicted curve 
seems to be highly sensitive to the magnitude of reference 

/ 

breakage stress ( G r ) and n b , both of which depend strongly 

on the nature of sand (Fig. 5.45). 

2, The continuous increase of B r upto unity with stress, 

following the hyperbolic law, is not reflected by the 
experimental data, for both weak and strong sands (Figs. 
5 . 45 , 5.46). The available trends suggest flattening of the 

Bj. versus stress curves around B r between 0.2 to 0.3, and 
thus the predicted rates of Increase of B r beyond the 
corresponding limiting stresses for weak and strong sands are 
much higher than the actual experimental values (e.g. Napa 
basalt and Mol sand). 
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3 . For the data on crushed granite and Napa basalt used by 
Hardin to calibrate his model, it would appear that <j'i is as 
good a stress path normalizing function as Cf b (Figs. 5.42 and 
5.43). Thus it may be suggested that the model could be 
worked out in terms of which is much easier to obtain 
directly. With this modification, equation 5.7 may then be 
vriten in the form 

c cr'i / a'r) n b 

i + ( tf'i/ cr'r) n b 


B. 


(5.7) 


For the weak angular sands used in this study, the relevant 
values Of 0 r n b obtained from Hardin’s guidelines are 

ff'r * 511 - 800 
n b = 0.939 - 1.230 • 

For these sands n a = 25 and h varies from 6.5 to 7.2. 

It is interesting to compare the range of B r values at stress 

2 

level of 176 Kg/cm for these sands using Eqn . 5.7 with 

different Or and n b values. The values of B r ' works out to be 
around 0.27 and for 0*r = 511 and n b =0.939 whereas it is 

around 0.13 for Or = 80000 and n b = 1.230. This range of 

B r (0.13 to 0.27) compares very well (Fig. 5.37) with the 

experimental range of B r (0.13 to 0.23). In this range of 
stresses, the modified Hardin’s model clearly provides a 
reasonably good prediction of particle crushing for weak 
angular sands. 

4. The relative breakage parameter (B r) , proposed by Hardin 
is a much superior index for overall crushing than some of 
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the indices reported by others. For example, one such index 
designated as the crushing coefficient (C c ) expressed as "the 
ratio of percentage of particles after being subjected to 

stress finer than D 10 of the original sand, and percentage of 

particles of the original sand finer than D 10 of the original 
sand" ( Datta et al , 1980), developed for Calcareous sand 

takes into account only the D 10 fraction. In contrast to 

this, B r takes into consideration the changes in size for the 
entire range of grains. Also in its definition, the particle 
size is considered in terms of B p (Fig. 5.22) and the total 
breakage, Bj takes into account particle size and gradation 
shape and mineralogy (Fig. 5.46a). Thus the index C r in some 
canes way provide min loading results as it considers changes, 
in one particular size only. Values of C c for Kalpi 1, Kalpi 
2, Ganga and Gang a with enhanced mica (Fig. 5.48) reveal that 
the C c for Kalpi 1 is higher than Kalpi' 2, whereas the B r 
value for Kalpi 2 is roughly 1.5' times higher than than that 
of Kalpi 1 . Kalpi 2 has larger grains which are aggregates 
of small particles cemented together by carbonate and 
therefore at same stress a higher particle breakage B x is 
expected which is in fact obtained. The value of B x for 
Kalpi 2 in roughly 1.7 times higher than that of Kalpi 1. In 
case of Ganga sand samples similar trends can be observed. 
Again C c for Ganga works out higher than for that same sand 
with enhanced mica content (655: mica as compared to 10 < in 
original Ganga). The B r value however follows the expected 


trend . 
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Thus it can be inferred that the Hardin's B r index is a 
superior indicator of crushing of grains and therefore should 
be Preferred over all other indices that consider changes of 

one particular size only. 

5. The shape number (n g ) used in the computations of 
reference breakage stress may be correlated with the Powers’ 
roundness index (Fig. 5.49). This emperical correlation 
between n g and roundness index would enable a better 
classification of grain shape which correlates with number of 
tangents on grain boundary as discussed earlier (Table 5.5). 
Thus the four parts shape description of Hardin can further 
be refined. 

5.3 COMPRESSIBILITY AND CONSTRAINED MODULUS OF SANDS 

5.3.1 Compressibility 

5. 3. 1.1 General 

Stress-strain behaviour of sands, as studied in the 
oedomet «r test ( Hendron , 1963 ; Roberts, 1964 ) has been 
considered to be of three stages (Lambe and Whitman, 1969) as 

below: 

a) Locking behaviour during which sand gets "stiff er and 

stiffer a3 the level of stress increases" upto 150 
Kg/cm 2 and "finally, a stage is reached in which 
already dense arrays (o£ particles) are being squeezed 
more tightly together as contact points crush, thus 
following a little more sliding”. 
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0.3 


Rou 

. 5.49 Relationship 
roundness . 
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Tafel e 5.5 ; 

Value 

of Shape 

Number (n s ) t 

Roundness 


Index 

(R) and 

Tangent 

(T) for 


Different Particle Shapes 



-» «» — ^ 

. , ^ ^ 



Par tide 
Shape 


Value 
of n 3 

Value of 
roundness 

Average 

tangent 

Angular 


25 

0.21 

11 

Subangular 

20 

0.30 

8 

Subrounded 

17 

0.41 

6 

Rounded 


15 

0.84 

4 
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b) Yielding behaviour as a "result • of fracturing of 
individual sand particles which permits large relative 
motions between particles" - at stresses equal to and 
greater than 150 Kg/cm 2 . Considerable particle 
degradation takes place during this stage. 

c) Fracturing and crushing, leading to tighter packing of 
particles and increase of number of particle contacts, 
thus the sand once again becomes stiff er and stiffer as 
the s tress increases still further”. 

These distinct stages of compression were observed for 
Ottawa sand which is medium to coarse and well rounded quartz 
sand. However, actual stress level at which compression 
behaviour of sand changes from locking to yielding depends on 
many factors such as mineralogy, particle size, shape, 
roughness, gradation and relative density. The stress level 
at which crushing and fracturing of particle becomes 
significant is designated as the "critical stress”. This 
stress is smallest when the particle size is large, the soil 
is loose and has uniform gradation, the particles are angular 
and the strength of the Individual mineral particle is low 
(Lambe and Whitman, 1969). The compressibility of sands 
depends on a number of material-related and stress-related 
factors (Clayton et al , 1985). Relative density of sands 
alone is not sufficient to characterize the performance of 
different granular soilds. For example, the compressibility 
of two sands having same mineral composition and relative 
density la significantly influenced by particle shape 

T 

(Holubec and D’Appolonia, 1973 Kapoor , 1985). 


In general 
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tl.o compressibility is high when the sand is loose, the 
grains are ™ e «l,rr, the coefficient of uniformity (C u) is low. 
the average particle site is small, inter-particle 
cementation is ueak and the ageing effects are insignificant. 
Stress history, stress level and stress path have significant 
influence on compressibility of sands CYudhbir and Rahim, 
1987). This then portrays the broad picture of deformation 
bchtiviour of sands over a wide range of stresses. Much needs 
to be done to quantify the effect of these varied factors. 

5.3. 1.2 Test Results 

The compression tost data on oedometer for the sands 

investigated under both loose and dense states has already 

been presented in Chapter 4. In order to highlight the major 

factors influencing compressibility, test results of loose 

samples only will be considered here. The stress versus 

2 

strain relationships upto a stress level of 176.32 Kg/cm for 
rounded quartz Standard sand, angular Kalpi sand with 
carbonate coating, subrounded to subangular Calcareous sand 
and angular micaceous Ganga sand are presented in Fig. 5.50. 
Results for two fractions of Kalpi sand i.e. Kalpi 1 and 

T 

Kalpi 2, and two variations of Ganga sand i.e. Ganga with 
enhanced mica and with mica removed, are also shown in this 
figure. The stress-strain behaviour of these sands, in the 
low stress ranges is illustrated in Figs. 5.51 and 5.52a. 
At a given stress level upto 10 Kg/cm , the axial strain 
increases from low compressible well-rounded quartz Standard 



i 


8 



° Standard 
□ Calcareous 
x Kalpi 2 
• Kalpi 
A Ganga 

Ganga (enhanced 
mica) 

+- Kalpi 1 
■ Gangafmica 
removed ) 


Fig. 5.50 Stress- strain relationship for the loose sands at 
high stress level. 
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l ° VC ' ry hl * h «.»pro.«ibl. mica- enhanced fraction of 
Gai.K.'i a.md 60’. mica). Calcareous. Kalpi 2, Ganea with 

mica removed , Kalpi. Kalpi 1 and Ganea sands arranged in the 

orcl or of increasing axial strain. 

In case of Kalpi 1 , Kalpi 2, and Calcareous sands 

(Fig. 5.52a), their compressibility undergoes a significant 

change under stress level greater than 25-30 Kg/cm 2 . Kalpi 

sand and its fraction Kalpi 1 have identical behaviour upto 

a stress level of 40 Kg/cm 2 (Fig. 5.50) after which Kalpi 1 

exhibits higher compressibility compared to the Kalpi sand. 

Kalpi 1 being coarse-grained and poorly graded (C u =l . 07 ) 

shown higher compression an compared to a well-graded 

(C u =4 . 81 ) Kalpi sample at the same stress level. This marked 

difference in the stress-strain behaviour of Kalpi and Kalpi 

1 is related to the higher degree of grain modification under 

stress as indicated earliet- (Fig. 5.35). The Kalpi 2 

fraction, consisting of large angular grains, shows a better 

2 ■ 

interlocking upto 25 Kg/cnT after which due to very high 
degree of grain modification (Fig. 5.26), it experiences much 
greater strain than Kalpi and Kalpi 1 at the same stress 
level. The Calcareous sand exhibited even better 

interlocking behaviour than Kalpi 2 sand upto stress level of 
25 Kg/cm 2 after which both these materials exhibited 
identical stress-strtain response. The Ganga— fraction with 
mica removed has slender angular quartz and feldspar grains 
and upto a stress level of 18 Kg/cm it has indicated lower 
strain than Kalpi sand. However, at higher stress levels 
this Ganga- f r act ion experiences much greater ax:ial strain 



248 


and upto a stress level of 18 Kg/cm 2 it has indicated lower 
strain than Kalpi sand. However, at higher stress levels 
this Ganga- £ r act i on experiences much greater axial strain 
than Kalpi, with identical stress-strain behaviour as that of 
Kalpi 1. Here again slender, angular quartz and feldspar 
grains of poorly graded sand have experienced much more grain 
modification than the well-graded Kalpi sand. The stress- 
strain behaviour depicted by Standard sand, Ganga sand and 
Ganga with mica enhanced has consistently shown the 

increasing strain with stress right upto stress level of 

2 

176.32 Kg/cm , However, the rate of variation of strain with 

stress, progressively decreases as stress level increases 

throughout the loading in case of Standard sand. However, a 

progressive flattering of the slope of stress-strain curve 

2 

has been noticed for Ganga sand at a stress level of 4 Kg/cm 
and for the Ganga sands fraction with mica enhanced at a 

stress level of 10 Kg/cm 2 . 

The distinct differences in the stress-strain response 
of sands under high stress conditions as brought out in the 
present study can be attributed to variations in grain shape, 
size, gradation and mineralogy. The effect of grain 
angularity and gradation on the behaviour of Standard, 
Calcareous, Kalpi and Ganga sand at comparable density under 
stress was independently brought out by Sabu (1988) in a 
series of well conducted plate load tests on these sands, the 
results for which are presented in Fig. 5.52b. These tests 
represent the behaviour of the sands at extremely low 
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contrast to practically non-dilatant ' Standard sand composed 
of nearly non-interlocking well-rounded quartz grains. Ganga 
sand being most comprressible undergoes larger settlement for 
the same stress and does not exhibit general shear failure as 
exhibited by Kalpi and Calcareous sands. This is in 
conformity with the compressibility considerations reported 
by Vesic (1963). Inset of region A in Fig. 5.52b clearly 
shows that the slope of the stress settlement curves is in 
the order of compressibility of these sands (Fig. 5. 51). 

Yudhbir and Rahim (1987) in a review of the available 
information on sand compressibility have observed that the 
compress ibi lily is high in general when the sand is loose, 
the grains are angular, C u is low, the average particle size 
is small , the particle surface is rough, the grain mineral 
strength is small, interparticle comentation is weak and 
ageing effects are insignificant. Except for the effect of 
ageing, the results obtained in the present study (Figs. 5. 50, 
5.51 and 5.52) clearly highlight the effect of grain shape, 
size, roughness, gradation and mineralogy. However, stress- 
strain curve of sands having same mineral composition is the 
end result of a complex interaction of particle shape, size, 
angularity and gradation. This is well illustrated in case 
of Kalpi, Kalpi 1 and Kalpi 2 sands. Ganga sand and its 
variations (one with mica removed and other with enhanced 
mica) clearly bring out the effect of mica content on sand 
compressibility. Terzaghi (1925), Moore (1971) and Yudhbir 
(1978) have brought out the influence of mica content on the 

Mitchell ( 1976) lias observed that 


compressibility of sands . 
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ad ci 1 1 1 o i i of mica s 1 n, t\ 1 f i c an 1 1 v inrrpaoaa +. u ^ . , „ , . , 

increases the compressibility 

Of poorly graded Hands only. The three Ganga samples have 

prat J rally identical particle size and gradation. The 

Standard sand exhibits the characteristics of a sand made of 

smooth well-rounded quartz grains: 

The sti ess-strain relationships for the sands investi- 
gated in this study are presented in Fig. 5.53 along with the 
results for a variety of sands reported by earlier workers. 
Hoi sand, Standard sand, Ottawa sand, Sandy point sand, and 
Dune sand etc. depict stress-strain-behaviour of quartz sands 
of low compressibility having rounded to well-rounded grains. 
The Ottawa sand, Sandy point sand and Mol sand results are 
for dense- state whereas the data for Dune sand and Standard 
3and is for loose state. Results for angular ground quartz 
show clearly the effect of d^g , the higher compression 
recorded for ground quartz with 0.42-0.84 mm size (curves 6, 
Fig. 5.53) grains is no doubt due to higher grain crushing. 
Kalpi , Kalpi 1, and Ganga sand with mica removed are 

predominantly angular sands with dominant quartz and 

feldspar, with Kalpi samples having 18-24* carbonate (as 
coating on quartz and feldspar grains) and Ganga fraction 
having 2 - 5 *. mica. Fine-grained Ganga fraction gives higher 
compression than coarse-grained Kalpi, Kalpi 1 and Kalpi 2 
sands . Ganga sand, and Calcareous sand, in the high stress 
range give more or less identical compression. Extremely 
high compression was observed for mica-enhanced Ganga sand 
fraction and crushed mica (Oberoi 1983). 
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On the basis of results shown in Fig. 5.53, it may be 
concluded that typical compression response for three 
categories of sands can be clearly characterized. Sands with 
well-rounded to subrounded quartz grains behave as materials 
0 f low compressibility. Angular quartz and feldspar-rich 
sands with or without carbonate cementation depict behaviour 
of moderately compressible materials while the micaceous (10- 
15%) carbonate sands (~ 95%)’ behave as highly 

compressible materials. The flaky mica sample gives the 
highest compression compared to all the sands shown in Fig. 
5.53. The grouping of sands of' low, moderate and high 

compressibility i s indicated in Table 5.6. Even though 

) 

factors such as grain size, shape, angularity, gradation, 
relative density and mineralogy se6ra to govern the compressi- 
bility of sands in a complex manner, the grain angularity and 
mineralogy (for the same relative density) appear to be the 
significant controlling factors. 

5.3.2 Constrained Modulus 

Janbu (1963, 1985) proposed a parabolic law governing 

d O' 

variation of constrained tangent modulus, M ( = ) with 


vertical stress in the form 




(5.12) 


where 


P is atmospheric pressure, 

B. 

m 0 is a modulus number, and 


m l 


is the exponent indicating rate of 
’’change of M with 0 1 , taken as 0.5 


by Janbu. 
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Tab 1 a 5 - 6 


Classification of Sands According 

Corapr esit Ibi 1 1 ty 


to Their 


-* • ■ ~ “ ■* - “■ «* - ~ * 








Sand a 

Shown 1 n 

Mo. 

Assigned 

Grain 

i Size 


Shape 



Flfi. 5.53 


d 10 

d 50 

d 60 

DR 

Remarks 

Mo 1 Sand 
(De Beer 1961) 

1 

— 

0.2 

~~ 

rounded 

d ens e 


Ottawa 

(Robert 1964 ) 

2 

0.45 

0.6 

0,61 

rounded 

dense 


Sandy Point Sand 
(Lamb a 6 Whitman 
1969 ) 

3 

— 

* — 

— 

rounded 

dens e 

low 

compre- 

ssibility 

Standard 

4 

0.3 6 

q. 48 

0 . 50 

well 
round ed 

loose 


Dune Sand 
(Obnrol 19 8.3} 

7 

0.15 

0.3 

0 .30 

round ed 

loose 



Ground quart* 5 

(0.15 to . 4 8 mm 
Lamfee & Whitman) 


either dense or 
rounded to subrounded 
or both 


Ground '0«artx 

6 

— 

— 

— 




(.12 - o . ni ) 







Moderat e 

Kalpi 

Kalp! 1 

8 

0.27 

1.0 

1.3 0 

very angular 

loose 

compress- 

ibility 

9 

0.95 

1 .0 

1.05 

very angular 

loose 


Gangs (mica removed) 

10 

0.10 

0.48 

0 .22 

angular 

loose 


Gang a 

11 

0.09 

0.18 

0.22 

angular 

loose 

High compress ; 
ibility 

Calcar town 

12 

0.27 

0.43 

0 .48 

sub rounded 
to subangular 

loose 


Gar *ga (enhanced mica) 

13 

0.08 

0.17 

, 0 . 20 

very angular 

loose 

Very high com* 
ss ibility 

Micaceous sand 
(Oberol , 1 983 ) 

14 

0.14 

0.7 

0.90 

very angular 

loose 
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Yudhbir and Rahim (1987) have examined the stress range 

for which this law may be aplicable for sands with varying 

relative density and particle angularity. Belloti et al 

\ 

(1985), on the basis of a set of chamber test results for two 

sands , proposed that constrained modulus (M) is related to the 

relative density and applied stress as 

1 

ra = m c P a CO'/Pa) 1 "! exp (m 2 .D R )' (5.13) 

where 

* 0 , E*a' anc * m l have the same meaning as 

defined in eqn. 5.12 and 
m 2 is an experimental coefficient indicating 
the rate of change of H with relative 
density (D r -j. 

Using Eqn. 5.12, the variations of modulus H with stress and 
relative density have been investigated for the sands chosen 
in this study. Evaluation of experimental coefficients m Q , 


and m 2 is illustrated in Figs. 5.54 to 5.61. Coefficient 
m Q is the value of H at CT = 1.0 Kg/cm 2 for I) R =0-0, m 2 ia 
obtained as natural logrithm of the ratio of the values of H 
at Cf' = 1.0 Kfi/cm 2 for D R = 1.0 and the value of m 0 . The 


coefficient is given by 

1 M 

m l n 

2.303 m 0 


(5.14) 


where 

M is the modulus at stress level of 10 
Kg/cm 2 for D R = 0.0 













CT , kg /chn* 

Fig. 5.60 Evaluation of m 0 , m 1 and m 2 coefficients for 
Kalpi sand. 



(T , kg/cm z 


Fig. 5.61 Evaluation of mg,m-| and m 2 coefficients for 
Ganga sand. 
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The values of M for relative densities of 0.0 and 1.0 
are obtained by extrapolation of test results as in Figs. 
5.54 to 5.57. Following this procedure, the values of m 0 , m-^ 

and m 2 have been worked out fop a variety of sands (Figs. 

5.58 to 5.61) as shown in Table 5.7. 

I 

5. 3. 2.1 Analysis of Experimental .Results 1 

Test results in terms of M 'versus stress for Standard, 

Calcareous, Kalpi , Kalpi 2 and Ganga sands are presented in 

Figs. 5.62, 5.63, 5.64, 5.65 and 5.66 respectively. The 

results for loose Calcareous sand are also shown in Fig. 5.64 

? 

to bring out the fact that upto stress level of 10 Kg/cm , 
Calcareous and Kalpi Sands show almost similar behaviour. 
The detail in this stress range is indicated in the inset of 

the figure. However, for stresses greater than 10 Kg/cm , 

modulus versus stress relationship for Calcareous sand is 
markedly different from that of Kalpi and this is mainly due 
to difference in mineralogy, particle shape and the fdct that 
grain modification in case of Calcareous sand begins at a 
lower stress level. The results presented are for two 
relative density values. Also superimposed on these curves 
are the respective model predictions on the basis of the 
relationships as in Table 5.7. It has been observed that 
Standard sand shows more or less a parabolic relationship 
between modulus and stress, upto quite a high stress level. 
The agreement between model prediction and experimental data 
upto 80 - 100 Kg/cm 2 only suggests that experimental 

coef f icients m c , m-j. and m 2 have been well determined. In 
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Table 5.7 : Model Parameters for Modulus of Variety of Sands 


Sands 

m 0 

itijL 

m 2 

M = 

IT 

i 0 cr mi e m2 


Ganga 

94 

0 .573 

0 . 533 

M 

=i 

9 4 00.573 

e 0.533 

Calcareous 

180 

0 . 378 

0.773 

M 

= 

180 G 0 378 

e 0.773 

Kalpi 

130 

0 .466 

0 . 900 

M 

= 

130 a 0 - 466 

00.900 

Standard 

220 

0 .560 

1.410 

M 

= 

220 a 0 - 560 

1.410 

e 

Ganga 

(enhanced Mica) 

34 

0.573 

0 .593 

M 

= 

34 00.573 

e 0.593 

Ganga 
(with mica 
removed) 

120 

0.371 

0 . 900 

M 


120 a 0 - 371 

e 0.900 

Ticino 

230 

0.370 

1 .470 

M 

ss 

23 0 a 0 - 370 

e 1.470 

Hokksand 

208 

0 .425 

1.360 

M 

= 

208 a 0,425 

e 1.360 

Ottawa 

185 

0.750 

1 . 373 

M 

= 

185 0 0 - 750 

el * 373 

Mine Tailing 

128 

0.371 

0 . 900 

M 

= 

i28a 0 - 371 

e 0.900 



, kg /cm 



effective confining stress (CF') for Standard sand 



, kg/cm 


Kalpi 2 

° RD 25% 
A RD 65% 


0 40 80 120 160 200 

CT', kg /cm 2 

Fig. 5.65 Relationship between constrained modulus (M) and effective 
confining stress (O’ 7 ) for Kalpi 2 sand 




kg /cm" 



Fig. 5.66 Relationship between modulus (M) and effective con 
fining stress (CT') for Ganga sand. 
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eano of angular sands of medium to high compressibility, the 
picture is quite different. The values of modulus increase 
wi t.h stress upto a certain stress level (critical stress) 
and then deviate from the values predicted on the basis of 
parabolic increase of modulus with stress- The trends for 
* stress values greater than critical stress vary differently 
for different sands. The value of critical stress also seems 
to depend on the relative density, particle size, shape, 
angularity and mineralogy of a sand. While in the case of 
Standard sand, the range for .critical stress is 80-100 
Kg/ cm 2 , the coresponding ranges for Calcareous, Kalpi, Kalpi- 


?, and Ganga 

n and s 

are 9-13, 15- 

28, 5-13 

and 

7-12 

re«p ee t i v e 1 y . 

Exe ep t 

for the case of 

SI: and ard 

sand , 

whch 

connial & of well-rounded quartz grains 

, the lower 

value 

s of 

c r i t i c a X s tres s 

are for 

loose and the 

higher values are 

for 


medium dense states respectively. In case of angular 
compressible sands, this result is consistent with the fact 
that they indicate cons iderabl e grain modifications at low 
.•stress ranges (Fig. 5.35). Harr and Grives (1980) suggested 
that the product of particle contacts and porosity is 
constant, and they have recommended a value of 3 for the same. 
In caMo of loose angular compressible sands, higher porosity 
would suggest smaller number of particle contacts with 
consequent higher contact stress concentrations leading to 
enhanced grain modification and crushing. This would also 
predict lower critical stress for loose sand as compared to 


dense sand. 
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In the case of Standard sand, it is expected that 
critical stress values for loose and dense states show a 
similar trend for large differences in relative density. The 
results in Fig. 5.62 are for a very close range of relative 
density values. In this case .any changes at low stress 
levels below (200 Kg/cm ) would- be mainly due to particle 
adjustments. In the case of angular,, medium to high 

compressibility sands, the deviation of experimental trends 
from the parabolic law is suggested to be due to the changes 
in the nature of sand with the onset of grain modification. 
At this low stress level, as indicated earlier, no grain 
crushing would have taken place and all observed changes in 
gradation are attributed to grain modification. This 

procedure to locate the point of deviation of experimental 
modulus versus stress trends from the predicted parabolic 
law, and the earlier direct evidence of grain modification 
(shape and angularity) at low stresses (Fig. 5.35) have been 
used for defining a critical vertical effective stress level 
( O ' c ) at which significant grain modification commences. 
This grain modification leads to significant changes in the 
sand compressibility (as shown earlier) and therefore the 
modulus versus stress relationships are observed to be at 
variance with the parabolic law proposed by Janbu and 
Bellotti et al . Yudhbir and Rahim (1987 ) have shown that m 0 
and mj values upto Cf^ are different from the subsequent 
behaviour of sand at higher stresses. The laws proposed by 
Janbu and Bellotti et al are based on test data obtained at 
stress levels less than 8-10 Kg/cm 2 . Results for sand tested 
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in thin stress range are shown in the Fig. 5.67. In this 
fslronn tango, except for loose, very highly micaceous or very 
loose carbonate sands, no grain modification is expected. 
Thus it may be stated that the modulus predictive laws (Eqns. 
5.12 and 5.13) tacitly assume no change in the nature of 
sand during compression. The value of 0 indicates a sort 
of threshold level upto which only the experimental constants 
m o' Itl l anc ^ ^2 ’ as evaluated for low stress level test data, 
are applicable. On the basis of limited data obtained in 

this study and the results available in literature, a 
tentative relationship between critical stress and relative, 
density is proposed (Fig. 5.68).. Although trends seen .in 
Fig. 5 . 6 » are consistent with the expected response due to 
differences in mineralogy and particle angularity, extensive 
data are required to clearly define these relationships for 
sands of low, medium and high compressibility.' 

5. 3. 2. 2 Effect of Grain Mineralogy on the Relationship 
Between Modulus and Stress 

Effect of mineral. content on grain modification 
crushing, and compressibility has been stressed in the 
earlier sections. The relationship between modulus and 
stress for three fractions of Ganga sand with varying 
percentages of mica content is indicated in Fig. 5.70. All 
the three samples have practically same d^g and comparable 
gradation curves. It will be seen that as mica. content 
i nor ears es , the constrained modulus at a given confining 
stress den eases. Variation of modulus number m Q (Eqn. 5.13) 



• Standard (RD=20%) 
O Calcareous (RD = 257o) 
O Kalpi 

A Ganga » 



(T' } kg/cm^ 


Fig, 5.67 Comparison of M vs CT' relationship at low stress level 
for Ganga, Kalpi, Calcareous and Standard sands. 
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and critical stress (Fig. 5.67) with mineral content is 

i ^ ^ 0 ^ " a t e < j in Fig. 5.69. A d d i t i o n a 1 c o n I; r olio d t o s t d a t a 

are needed to establish quantitative relationships between 

constrained modulus and mineral composition for a variety of 

. * 

sands. 

5.3.3 Applications of Test Result to The Estimation of Cone- 
Penetration Resistance in Sands 

Based on the results of grain shape, angularity, 

modification, crushing, and compressibility studies discussed 
earlier and the shear test data on Standard, Kalpi and Ganga 

sands reported by Phool Chand (1988), it is possible to 

predict. the cone penetration resistance, t q r (measure of 
ultimate bearing capacity), as a function of relative density 
and confining effective stress for sands with varying 
compressibilities. In this section, effect of grain 

angularity on angle of shearing resistance, reported by Phool 
Chand (1988) and other workers will be first reviewed and 
then utilised to predict q c on the basis of Janbu’s theory 
(196), 1985). The effect of sand compressibility on the mode 

of failure will be incorporated in the form of angle P (Fig. 
5.71) to compute the appropriate values of bearing capacity 
factor N q of a sand based on its angle of shearing 
res i s t attc e . 

5. 3. 3.1 Calculation of q c for Sands with Varying ' Degree of 

Compressibilities 

Janbu and Senneset (1974) suggested a value of (-15) 
for P for sand' of low compresibili ty which exhibit general 
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shear failure mechanism around the cone. Robertson et al . 

(1981) (Fig. 5.71-Inset) on the basis of all available data 

for test chamber studies, concluded that the values for 

sand of low compressibility (at a given relative density and 

confining stress) vary from -15° to -32°. For the purpose of 

computing q c values, it was assumed that, as a first 

approximation, (3 = -15° to -32° be taken to represent failure 

in case of rounded quartz sands of low compressibility and 

|3=+15° be taken to represent the failure mode in case of 

sands of high compressibility. For sands of moderate 

compressibility [3 = 0.0 was adopted. Ideally in case of 

compressible sand, variable [3 values as suggested by Vesic 

( 1963 ) (Fig. 5.71) would be more appropriate. Since the 

objective of thi.s exercise is to develop a comparative feel 

for the overall trends of q c versus effective stress for 

sands with low, medium and high compressibility, a constant 

value of [3 was chosen. Values of q c were computed for each 

of the representative sands (Standard, Kalpi and Ganga) on 

the ban i s of <p' versus relative density variation at 

different. stress levels as shown in Figs. 5.72 and 5.73 

(Phool Chand 1988). The variation of cf)' with relative 

density for three sands in the stress ranges investigated is 

1 

presented in Fig. 5.74. For a' given relative density and 
confining effective stress, angular sand gives higher (p as 
compared to well-rounded sand. This is consist ant with 
results compiled by Uinterkorn and Fang (1986). 
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Fig. 5. 73 Relation between peak friction angle and 
Relative Density for Gdnga and Kalpi at 
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Following Bel lot ti et al (1985), q c may be expressed in 
terms of effective stress and relative density in the form 

q c ~ c o 0 Cl exp (C 2 .D r ) (5.15) 

where 

C Q » C 1 and c 2 are experimental coefficients. 

From the computed values for a range of stress and 
relative density values, the coefficients C Q , C 1 and C 2 were 
evaluated for each of the sands in the same manner as adopted 

for the computation of m 0 , ir^ and m 2 as detailed earlier. 

The C () , and C 2 values proposed for sands of low, medium 

and high compressibility are 1 given in Table 5.8. In 

addition, the corresponding values of m Q , m-^ and m 2 

characterizing the modulus for these sands are also presented 
in this Table. The classification of sands in these three 
categories based on their compressibility is already 
presented (Table 5.6, Fig. 5. 53). 

5. 3. 3. 2 Comparison Between Predicted and Measured Modulus 
versus Stress Relationships for Sands 
On the basis of predictive laws (Table 5.7), values of 
modulus at various stress values for a wide range of relative 
densities were computed for standard, Calcareous, Kalpi and 
Ganga sands. The predicted and measured results for these 
four sands are presented in Figs. 5.75 to 5.78. As can be 
noticed, there is a good agreement between the predicted and 
the observed values. In Fig. 5. 79, the behaviour of different 
sands at comparable relative density value is illustrated. 
Test. data fr,jf>m Fhool Chand (1988) are also included. Using 
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Table 5.8 : Suggested Model Parameter for Constrained 


1 

Modulus 

CM) and 

Cone Penetration 

Resistance (q c 

Type of Sands 

Measure 

of M Kg/cm 2 

Measure 

of q c 

Kg/cm 2 

mo 

1*1 

m 2 

Co 

C 1 

C 2 

Low Compress- 
ibility Sands 

220 

0 . 37 

1.41 

21 

0.51 

2.96 

Moderat e 

Compressibility 

Sands 

152 

0 . 42 

0 . 84 

15 

0.64 

2 . 61 

High 

94 

0 .57 

, 0 . 53 

10 

0.87 

1.99 


Compressibility 

Sands 



Modulus (M) , kg /cm 


2600 



R D = 1 00 % 

__ 

Predicted S' 

o 

Experimental s' 


2400 
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Fig. 5.7 ffect of relative densit 
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Fig. 5.79 Effect of RD and (T^on modulus for Ganga, Kalpi, standard 
and calcareous sands. 
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m o * m l anc * m 2 values as recommended in Table 5.8, a 
generalized comparative relationships of modulus and stress 
for sands of diferent compressibilities are presented in Fig. 
5 . 80 . Results are compared for relative densities of 0.40 
and 0.80 for all the three categories of sands. 

Similarly using C c , and C 2 values (Table 5.8), a 
generalized comparative relationship of q c and stress for 
three categories of sands are shown in Fig. 5.81. The 
relationship as reported by Robertson et al (1983) on the 
basis of sand compressibility is also shown as insert in this 
Figure. As can be seen, the patterns as predicted from the 
present study are in agreement with trends reported by 
Robertson and coworkers. It may be indicated that the 
patterns for the trends in Fig. 5.80 and 5.81 in terms of 
modulus and q c strikingly similar as influenced by sand 
compressibility. 

The values of m Q , m^ and m 2 as well as C Q , and C 2 
have been correlated with particle angularity. As indicated 
earlier, low, moderate and high compressibility of sands is 
directly related to particle angularity (for the same mineral 
composition). Particle angularity was related to Powers’ 
roundness index (Table 2.2). Relationship between 
coefficients m Q , , m 2 , C Q , and C 2 and Powers’ roundness 
index have been presented in Fig. 5.82. Although additional 
data points would be needed to evolve the exact nature of 
these relationships, the overall trends are distinct. 
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5. 3, 3. 3 Relationship Between 
Resistance 


l c a nd Angl, 


of 


Shearing 


Broms (1986) has reviewed t-t- 

10 available data on various 

relationships between q r . and &>' 

1 i3atlds - The influence of 

effective overburden pressure , 

U<! sand compressibility on 

such a relationship has been indict- , 

1 cat ed . I„ Fie . 5 . 82Cb)i the 

predicted trends for sands of l 0w 

and high compressibility 
from the present study are present-^ 

along with the reported 

values by Meyerhof ( 1976). The D ,- i- 

1 evicted trends are on the 

basis of the following relationships- 
tan Cp' = .64(cr') _ ' d ^ e -56DR 

„1 . 9 9 DR 




io ( a 1 ) ■ 87 


tan (p' ~ 0.5 ( cr ' ) ’ 07 


2 i ( a ' ) • 51 


. 83DR 


2 . 9 6 DR 


For sands of 
high compressi- 
bility 


For sands of 
low compressi- 
bility 


(5.16) 

(5.17) 

(5.18) 

(5.19) 

Relationships (5.16 and 5.18) ar* t- 

based on drained direct 

shear test data (Phool Chand 1988) ti 

' ille results shown in Fig. 

5.82(b) clearly indicate the infi,, 

luence of effective over- 
burden stress and sand compr ess i k ; i • , 

sibUxty. The relationship 

proposed by Meyerhof appears to h* , 

• 6 a PPlicable for sands of 

medium to low compressibility at i rtt , 

row efective stresses. As 

such there cannot be a unique ret-.-*-: 

•a-tionship between q c & 0 

irrespective of sand compr ess ib’i 1 ; t- , r 

lty and effective overbuden 

pressure . 


5. 3. 3. 4 Relationship Between M and n 

q c 

The modulus and q„ values as ^ . . 

c a Predicted, by using m Q , mj , 

mo, C n , Ci, and Co recommended for thr- Q „ 

£ ° 1 L tnree categories of sands 

have been compared over a wide ranoa ,, .. , . . 

aee °f relative densities. 
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.similar plots have also been reported by Robertson et al 
(1983). The relationships for sands for low, moderate and 
high compressibility are shown in Figs. 5.83 to 5.85 together 
with the variation of a (=M/q c ) versus relative density at 
stress level of 1 and 10 Kg/cm^. XJithin this stress range, 
the variation of (X values for relative density values of 0.0 
and 1.0 is shown in Table 5.9. While q c versus stress and M 
versus stress relationships at a given relative density are 
markedly different for sands of different compressibility, 
the values of (X for all the three categories of sands are 
within a narrow range for a given relative density. The 
average curves fitted through the data of a values for all 
the three categories of sands ; are presented in Fig. 5.86. 
The relationship between 0( and relative density (DR) may be 
represented as 

CX =10CT / -(- 22 ) e ~ ( 1 . 6DR) (5.20) 

t 

The above expression is valid for normally consolidated sands 

2 

upto a stress level of 10 Kg/cm J and the (X values are within 
the range reported in literature (Mitchell and Gardner 1975). 
Robertson et al . (1983) suggest the values for 0( in the 
range 3-11. The above equation thus appears to define a 
somewhat unique relationship between (X and relative density 
for normally consolidated sands. 

Thus in this section, the following generalizations 
have been arrived at: 

1 . Sands, depending on their mineralogy, particle size, 
shape, angularity and gradation can be grouped into three 
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Low Compressibility 
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300 


ft Values for Sands with Ranging 
ility 
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10 Kg/cm 2 

DR = 0 

DR = 1.0 
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4.8 
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broad categories (Table 5.6) as (i) ' sands of low 
t.osBpi' esa i toi 1 i t y , (ii) sands of moderate compressibility 
ami (iii) sands of high compressibility. 

2 . Even though factors such as grain size, shape, 
angularity, gradation, relative density and mineralogy 
seem to govern the compressibility of sands in a complex 
manner, the particle angularity and mineral content (for 
same relative density) appear to be the two most 
important factors. 

3. Stress-strain behaviour of compressible sands is markedly 
different in low stress ( 0 ' < 10 Kg/cm 2 ) and high stress 

( d ' > 30 Kg/cm 2 ) ranges. 

4. Constrained modulus (h) is shown to be related to 
effective stress and relative density through the 

expression 

H = m D a 1 e ‘ 

This is cons i a t ant with the model proposed by Bellotti et 
a l (1985). The values of m Q , m^ and m 2 are a function of 
sand compressibility (Table 5.8). 

5. The coefficients m Q , m-^ and m 2 as evaluated are valid for 
being the critical stress at which significant 
grain modification begins. Based on limited data, 

relationship between critical stress and relative 

! 

density for variety of sands has been proposed. The 
values of m 0 , m-^ and m 2 will be different for the 

dif event stress levels exceeding (j ' c , when grain crushing 


«: out me nee a . 



hroflrJ categories (Table 5. 6> as (i) sands of lov 
compressibility, (ii) sands of 'moderate compress ib il itrv 

5 

am! (iii) sands of high compressibility. 

Even though factors such as grain size, shape, 
angularity, gradation, relative density and mineralogy 
seem to govern the compressibility of sands in a complex 
manner, the particle angularity and mineral content (for 
same relative density) appear to be the two mos ~t 
important factors. 

Stress' strain behaviour of compressible sands is markedly 
different, in low stress ( CT 1 < 1U Kg/cm^) and high stress 
( Cj >10 K g/ cm^'J ranges. 

Constrained modulus (H ) is shown to be related to 
effective stress and relative density through. the 

express ion 

n - » c a-'"'! , ‘“! dr 

Thin i ; i cona Istant with the model proposed by Bellotti et 
a] (1 9R!> ) . The values of m () , her and. m 2 are a function of 

n.iiid compr «>a al b il ity (Table 5.8). 

Tin* coefficients m 0 , m j and m 2 as evaluated are valid £o*r 
(J / <q'-,0 c '> ei »>e the critical stress at which significant 
grain mod if ica ti on begins. Based on limited data , 
t el at. ion jih ip between critical stress and relative 
d on city for variety of sands has been proposed. The 
values of m <} , mj and m 2 will be different for the 
fill «>i ant. .* ? t * i* ms levels exceeding . when grain crushing 
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(> . Baaed on the effect of grain angularity on (f)' and sand 
compressibility on mode of failure around static 
penetration cone, the relation of q c with the effective 
stress and relative density is of the following form 
q c « C c cr' Cl e C2DR 

Beilotti et al (1985) proposed a similar form. The 

values of C 0 , and C 2 are proposed to be function of 

sand compressibility (Table 5.8). The values of C Q , C-^, 

/ 2 

and C2 as evaluated are valid for Cf <10Kg/cm . 
Coefficients C Q , C , and C 2 like m G , m-^ , and in 2 have been 
correlated with Powers’ roundness index. 

7 . It has been brought out (Fig. 5.82(b)) that q c versus 4 >' 
relationship is governed by sand compressibility and 
effective overburden stress. There is thus no unique 
relationship between q c and <£> applicable to all sands as 
suggested in the literature. 

0. Even though modulus versus stress and q c versus stress 
t t o» a given relative density) are shown to be intimately 
gov <m n ed by the sand compressibility, their ratio (X = 
appears to be uniquely related primarily to relative 
density and to a lesser degree to the stress. It is 
proposed that, as a first estimate, (X may be evaluated 
1 1 urn t he expression 

a = 11) 0'"C-22) ,-!■< »* 

The value of modulus may thus, be estimated using q c as 

obi, lined from the cone test. 



CHAPTER VI 


CONCLUSIONS 

The pattern of deformation of sands under different 
stress conditions is controlled by several parameters such as 
the grain size, shape, angularity, state of packing and 
mineralogy. While earlier studies have emphasized on the grain 
size and state of packing, importance of grain shape and 
angularity as well as the mineralogy has not been 
satisfactorily worked out. On the basis of the detailed 
investigations under taken in the present study to examine the 
interrelationships between various parameters operative during 
deformation of sands and quantitative evaluation of the same, 
the following aspects are established: 

6.1 The particle shape characterization for each of the sand 
fraction in the sand used in the present study have been 
carried out using the Zingg diagram and shape factor. On 
the basis of Zingg diagram and the shape factor versus 
sphericity relationship (Fig. 5. 1b), sands can be 
classified with respect to grain shape. 

6.2 The tangent count technique, using the image analyzer 
provides an accurate and highly reliable method to 
quantify particle surface characteristics of a variety of 
sands. Typical average tangent count distributions for 
subrounded to rounded, subangular , angular and very 
angular grains have been developed (Figs. 5.17 and 5.18) 
and are recommended as a useful tool for classification 
of sands in respect of grain angularity. The average 
tangent count has been employed to provide a systematic 



method for evaluation of Powers’ roundness index, R 
(Fig. 5.1c) which can be used to classify sand based on 
the grain angularity. 

6.3 Relationship between angularity and size for a particular 
alluvial sand can be demonstrated in terms of average 
tangent count per Unit circumferential length of 
particles for each sieve fraction. An inverse relation- 
ship betwen particle, angularity and size has been 
established in case of alluvial sands (Fig. 5.1a). 

6.4 It has been shown that while average tangent count is 
useful for comparing sands having grains of different 
angularity, the particle angularity defined as 

average number of tangents in a sieve 
fraction 


average circumference of the particles 
in that fraction 

is a good measure of variation of angularity with 
particle size for a given sand. Particle angularity for 
alluvial sands is shown to decrease with increase in 

grain size. 

On the basis of these conclusions, it is recommended that in 
addition to the usual classification tests for sands 
(e max ,e mi n» e, grain size distribution etc.), the shape and 
angularity characteristics of sand grains should also be 
evaluated using the technique demonstrated in this study. 

6.4 For sands with a wide range of particle gradation 
characteristics a relationship between total potential 
for particle breakage <B p) (Fig. 5. 22) and d 50 (mm) has 
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been established ass 

Bp a 1.02 lOg djQ + 1.15 

6.6 At low stress levels, (<22 kg/cm 2 ) the total particle 
breakage (Bj), increases with grain angularity (h/A) for 
a wide variety of sands (Fig. 5.26a). 

6.7 For same minerology and particle angularity, total 
particle breakage (Bj) (Fig. 5.19) of a sand at a given 
relative density increases with dgg (Fig. 5.26b) over a 
wide range of stress. 

6.8 Based on compression tests at different stress levels 

(11.02, 22.04, 88.16 kg/cm ) along with subsequent 

particle shape analysis (Fig. 5. 27), tangent count 
distribution (Figs. 5.28 to 5.30), smooth edge 

distribution (Figs. 5.31 to 5.33), particle angularity 
(Fig. 5.34) and particle breakage analysis (Fig. 5.35), 
it has been established that there exists a critical 
stress level, 0^ at which grains of alluvial sands 

undergo modification. This critical stress is much less 
than that required for actual grain crushing and its 
value depends on mineral composition, size, shape, 
angularity, gradation and relative density of the sand. 

On the basis of limited test data available, a tentative 
relationship between critical stress 0 C and relative 
density for a variety of sands has been proposed (Fig. 
5.68). 

6.9 The grain modification continues upto a stress level of 
about 80 Kg/cm 2 after which actual particle breakage 
(crushing) starts. The variation of particle breakage 
(expressed in terms of B r ) with stress is in general 
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agreement with the pattern suggested by Terzaghi and Peck 
(1943) and De Beer (1963). However, the B r versus o' c 

relationship depends on the nature of the sand. 

6.10 For weak sands (Ganga, Kalpi, Calcareous, Crushed Granite 

and Napa Basalt etc.) the change of B r with stress may be 

considered in four phases (Fig. 5.37). 

P 

i) For stress <20 Kg/cm , the rate of increase of B r 
with stress is high and particles undergo 

modification. 

ii) For stresses between 20 and 100 Kg/cm 2 , the 

modification continues but B r increases at a lower 

rate compared to that in phase (i). 

iii) For stresses greater than 100 Kg/cm 2 , significant 

particle crushing begins and B r experiences the 

highest rate of increase upto a limiting stress 

2 

which is likely to be smaller than 1000 Kg/cm (Fig. 

5.39) . 

iv) For stresses greater than the limiting stress, B r 
increases at a lower rate than that in phase (iii). 

6.11 The variation of B r with stress in case of strong sands 
(well rounded quartz. Standard and Mol Sand, etc.) also 
occurs in four phases (Fig. 5.37) 

I) For stresses upto 200 Kg/cm 2 , the increase of B r is only 
marginal, with some modification/crushing in case of 
subrounded grains. 

11) For stresses between 200 to 400 Kg/cm , B r increases 
relatively faster than that in phase 1. Grain crushing 

has started in this phase. 



Ill) The highest rate of increase of B r occurs for stresses 
between 400 to 1000 Kg/cm 2 (Fig. 5. 39). 

IV) For stresses greater than a limiting stress (1000 
Kg /cm ) rate of increase of B r is less than that in 

phase III. 

6.12 The behaviour of subangular to subrounded sands (as in 
the case of Toyoura Sand) is intermediate to that of weak 

c" 

and strong sands. 

6.13 Relative breakage B r increases with mica content 
(Fig. 5. 40a) . 

6.14 There exist a relationship between B r and quarts content 

(also quartz plus feldspar) for different sands 

(Fig .5.40b) . 

6.15 Based on the data mainly for Calcareous and Kalpi sands, 
the effect of carbonate on particle crushing has been 
evaluated (Fig. 5. 40c). 

6.16 The relative breakage parameter B r is a much superior 
index of particle crushing (modification) than some of 
the other indices available in the literature. 

Based on grain size, shape, angularity, mineralogy, gradation 
and relative density, natural sands experience particle 
degradation to varying degrees under stress. Studies related 
to particle degradation for weak sands under stress should 
distinguish between grain modification which starts at a 
critical stress level, and grain breakage (crushing) which 
commences at a stress level of about 100 Kg/cm 2 . The 
results of changes in particle gradation under stress should 
be expressed in terms of Hardin's B r parameter. 



309 


6.17 Based on the test results and compar isons with the 
predictions from Hardin's Model, the following inferences 

are made* 

a) The hyperbolic law to predict relative breakage under 
stress is not unique for all types of sands. The 
shape of the predicted curve seems to be highly 
sensitive to the magnitude of the reference breakage 
stress, a' r , (Fig. 5. 44a). 


b) Normalization 

of 

the 

test results 

under 

different 

stress paths 

can 

be 

achieved using 


and a' b 

(Fig. 5.42) .• 

The 

model parameters 

can as 

well be 


worked out in terms of cf^ rather than which 

requires reliable values of K Q for the sand. The 
ratio (j' b / CT^ varies between 0.9 to 1.5 over the 
range of K Q values for normally consolidated sands 
(Fig. 5. 44) . 

c) The range of B r values attainable in practice is 
likely to be limited to 0.5 (Figs. 5.36, 5.37, 5.38, 

5.41, 5.45) even when, the stress level is as high as 

1000 Kg/cm 2 . 

6.18 The shape number (n s ) of, Hardin's classification for the 
computation of reference breakage stress, is correlated 
with Powers' roundness index (Fig. 5. 47). This enables a 
better classification of grain shape (angularity) which 
correlates with average number of tangents on grain 
surface (Table 5.5). Hardin's four part shape 

(angularity) classification can thus be refined. 
Parameters for Hardin's hyperbolic law may as well be worked 
out in terms of <j\ rather than 0 b . The value of Hardin's 
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* h *P« number (n f ), for a sand may be obtained from the 

relationship between n s and Powers* roundness index (R) 

<Fifl.3.49). 

6. 19 Stress-strain behaviour of compressible sands is 

markedly different in low stress ( &'< 10 Kg/cm 2 ) and 

high stress ( O' > 30 Kg/cm 2 ) ranges (Figs. 5.50, 5.51, 

3.52) . 

6.20 Of the factors such as grajn size, shape, angularity, 
gradation and mineralogy that govern the compress- 
ibility of sands (for a given relative density) in a 
complex manner, the particle angularity and mineralogy 
appear to be the two most important ones (Fig. 5. 53). 

6.21 On the basis of mineral ogy ,' par t i c 1 e size, shape, angu- 
larity and gradation, sands can be grouped (Table 5.6) 
into three broad categories in terms of their 
compressibility into 

(i) sands of low compressibility, 

(ii) sands of moderate compressibility, and 

(iii) sands of high compressibility. 

6.22 The relationship between the constrained modulus (M) and 

effective stress as proposed by Janbu (1963, 85) and 

Bel Dili et al (1985) are valid in the stress range less 

than the critical stress G' c . At stresses beyond this 

2 

critical value (and less than 100 kg/cm ), the weak 
sands undergo significant changes due to grain 
modification and the relationship between M and 
effective stress as obtained from experimental data in 
the present study significantly differs from those that 
predicted earlier (Figs. 5. 62 to 5.66). 
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6.23 At stresses less than the critical stress, M is proposed 
to be related to effective stress and relative density 
by the following relationship 

M = m 0 *""1 exp(mg Dp) 

The values of m 0 , m-j and mg for sands of varying 
compressibility are also estimated (Table 5.7). 

6.24 The effect of compressibility on M - (j * relationship for 
sands has been highlighted (Fig. 5.80). 

6.25 Based on effect of grain angularity on angle of shearing 

resistance and sand compressibility on mode of failure 
around static penetration cone, q c has been computed at 
different values of ^and relative density. The depen- 
dence of relationship between q c and (j * on sand 

compressibility has been brought out (Fig ,5.81a) . These 
trends agree with those presented by Campenella et al . 
(Fig .5.81 in set ) . 

6.26 A relationship between q c , relative density and 

effective stress, is worked out in the present study, is 
similar to the one reported by Belotti et al (1985) and 
is of the form 

q c = C Q (T /C 1 exp (Cg D R ) 

The values of C D , and Cg are proposed (Table 5.7) for 


sands 

of varying 

compressibility. As <p' decreases 

at 

high 

stresses. 

the proposed relationship and 

the 

coefficients C Q , 

C. 1 and C 2 are applicable strictly 

for 


2 

stress levels less than 10 Kg/cm . 

6.27 Even though the relation of CT ' with M as well as with 
q c (at a given relative density) are shown to be 
intimately governed by the sand compressibility, the 



ratio, 0 “ M/q c however' appears to be uniquely related 
to relative density and effective' stress. The range of 
values of (X (3— 11) is in agreement with the reported 
values of Compenella et al. 

6.28 It has been brought out (Fig. 5. 82(b) ) that q c versus 
<p' relationship is gov erned by sand compressibi 1 ity and 
effective overburden stress, and hence no unique 
relationship between q c and Cp' is applicable to all sands 
as suggested in the literature. 

6.29 It is proposed that as a first estimate, CX may be 
evaluated from the following relationship! 

(X = W-°* 22 exp (-1.6 DR) 

The value of M thus can be estimated on the basis of q c 
obtined from the cone test. 

The compressibility of sands is strongly controlled among 
other factors by mineralogy and grain angularity. At low 
stresses, the stress-strain behaviour of weak sands is mainly 
governed by grain interlocking. At higher stress 030 
Kg/cm ) , the stress-strain behaviour is significantly 
different from that at low stresses due to the grain 
modification and crushing. This distinction should be 
clearly recognized. 

The power relationship of modulus with stress and 
relative density tacitly assumes that the sand remains 
unchanged (without any change in its stress-strain behaviour) 
during loading to high stresses. The fact that the stress- 
strain behaviour of the sand undergoes changes due to 
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particle modification and crushing at high stresses, these 
power laws and the suggested coefficients m Q , and m 2 are 
strictly applicable for stress levels less than critical 
stress. Both q c versus o' and M versus 0 ' relationships are 
strongly controlled by the compressibility of the sand. 
Therefore any relationship developed for a particular sand is 
only ’ applicable for sands of similar compressibility 
characteristics. Correlation developed for relatively 
incompressible rounded to subrounded artificial quartz sands 
is therefore not to be applicable to natural sands which are 
as a rule subangular to angular in character and therefore 
compressible in behaviour. 

The expression cx=10CJ , ~^*^ exp ( — 1.6 DR) may be employed to 
estimate the value of M on the basis of cone penetration 
resistance value (q c ) commonly obtained as part of soil 

investigation program. 
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